Complex Breakpoints and Template Switching Associated with Non-canonical Termination of Homologous Recombination in Mammalian Cells

One-ended and two-ended breaks at nickase-broken replication forks

Replisome collision with a nicked parental DNA template can lead to the formation of a replication-associated double strand break (DSB). How this break is repaired has implications for cancer initiation, cancer therapy and therapeutic gene editing. Recent work shows that collision of a replisome with a nicked DNA template can give rise to either a single-ended (se) or a double-ended (de)DSB, with potentially divergent effects on repair pathway choice and genomic instability. Emerging evidence suggests that the biochemical environment of the broken mammalian replication fork may be specialized in such a way as to skew repair in favor of homologous recombination at the expense of non-homologous end joining.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. We used inverse PCR of non-B microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures integrated at a common ectopic chromosomal site to show that these non-B DNAs generate highly mutagenized eccDNAs by replication-dependent mechanisms. Mutagenesis occurs within the non-B DNAs and extends several kilobases bidirectionally into flanking and nonallelic DNA. Each non-B DNA exhibits a different pattern of mutagenesis, while sister clones containing the same non-B DNA also display distinct patterns of recombination, microhomology-mediated template switching and base substitutions. Mutations include mismatches, short duplications, long nontemplated insertions, large deletions and template switches to sister chromatids and nonallelic chromosomes. Drug-induced replication stress or the depletion of DNA repair factors Rad51, the COPS2 signalosome subunit or POLη change the pattern of template switching and alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA double strand breaks to account for the generation and circularization of mutagenized eccDNAs and the appearance of genomic homologous recombination deficiency (HRD) scars. These results may help to explain the appearance of tumor eccDNAS and their roles in neoantigen production, oncogenesis and resistance to chemotherapy.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. In tumors, highly transcribed eccDNAs have been implicated in oncogenesis, neoantigen production and resistance to chemotherapy. Here we show that unstable microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures generate eccDNAs when integrated at a common ectopic site in human cells. These non-B DNA prone microsatellites form eccDNAs by replication-dependent mechanisms. The microsatellite-based eccDNAs are highly mutagenized and display template switches to sister chromatids and to nonallelic chromosomal sites. High frequency mutagenesis occurs within the eccDNA microsatellites and extends bidirectionally for several kilobases into flanking DNA and nonallelic DNA. Mutations include mismatches, short duplications, longer nontemplated insertions and large deletions. Template switching leads to recurrent deletions and recombination domains within the eccDNAs. Template switching events are microhomology-mediated, but do not occur at all potential sites of complementarity. Each microsatellite exhibits a distinct pattern of recombination, microhomology choice and base substitution signature. Depletion of Rad51, the COPS2 signalosome subunit or POLη alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA breaks for the generation and circularization of mutagenized eccDNAs and genomic homologous recombination deficiency (HRD) scars.

Concatenation of Transgenic DNA: Random or Orchestrated?

Generation of transgenic organisms by pronuclear microinjection has become a routine procedure. However, while the process of DNA integration in the genome is well understood, we still do not know much about the recombination between transgene molecules that happens in the first moments after DNA injection. Most of the time, injected molecules are joined together in head-to-tail tandem repeats-the so-called concatemers. In this review, we focused on the possible concatenation mechanisms and how they could be studied with genetic reporters tracking individual copies in concatemers. We also discuss various features of concatemers, including palindromic junctions and repeat-induced gene silencing (RIGS). Finally, we speculate how cooperation of DNA repair pathways creates a multicopy concatenated insert.

Noncanonical outcomes of break-induced replication produce complex, extremely long-tract gene conversion events in yeast

Long-tract gene conversions (LTGC) can result from the repair of collapsed replication forks, and several mechanisms have been proposed to explain how the repair process produces this outcome. We studied LTGC events produced from repair collapsed forks at yeast fragile site FS2. Our analysis included chromosome sizing by contour-clamped homogeneous electric field electrophoresis, next-generation whole-genome sequencing, and Sanger sequencing across repair event junctions. We compared the sequence and structure of LTGC events in our cells to the expected qualities of LTGC events generated by proposed mechanisms. Our evidence indicates that some LTGC events arise from half-crossover during BIR, some LTGC events arise from gap repair, and some LTGC events can be explained by either gap repair or "late" template switch during BIR. Also based on our data, we propose that models of collapsed replication forks be revised to show not a one-end double-strand break (DSB), but rather a two-end DSB in which the ends are separated in time and subject to gap repair.

Mako: A Graph-based Pattern Growth Approach to Detect Complex Structural Variants

Complex structural variants (CSVs) are genomic alterations that have more than two breakpoints and are considered as the simultaneous occurrence of simple structural variants. However, detecting the compounded mutational signals of CSVs is challenging through a commonly used model-match strategy. As a result, there has been limited progress for CSV discovery compared with simple structural variants. Here, we systematically analyzed the multi-breakpoint connection feature of CSVs, and proposed Mako, utilizing a bottom-up guided model-free strategy, to detect CSVs from paired-end short-read sequencing. Specifically, we implemented a graph-based pattern growth approach, where the graph depicts potential breakpoint connections, and pattern growth enables CSV detection without pre-defined models. Comprehensive evaluations on both simulated and real datasets revealed that Mako outperformed other algorithms. Notably, validation rates of CSVs on real data based on experimental and computational validations as well as manual inspections are around 70%, where the medians of experimental and computational breakpoint shift are 13 bp and 26 bp, respectively. Moreover, the Mako CSV subgraph effectively characterized the breakpoint connections of a CSV event and uncovered a total of 15 CSV types, including two novel types of adjacent segment swap and tandem dispersed duplication. Further analysis of these CSVs also revealed the impact of sequence homology on the formation of CSVs. Mako is publicly available at https://github.com/xjtu-omics/Mako.

Thermodynamic analysis of DNA hybridization signatures near mitochondrial DNA deletion breakpoints

Broad evidence in the literature supports double-strand breaks (DSBs) as initiators of mitochondrial DNA (mtDNA) deletion mutations. While DNA misalignment during DSB repair is commonly proposed as the mechanism by which DSBs cause deletion mutations, details such as the specific DNA repair errors are still lacking. Here, we used DNA hybridization thermodynamics to infer the sequence lengths of mtDNA misalignments that are associated with mtDNA deletions. We gathered and analyzed 9,921 previously reported mtDNA deletion breakpoints in human, rhesus monkey, mouse, rat, and Caenorhabditis elegans. Our analysis shows that a large fraction of mtDNA breakpoint positions can be explained by the thermodynamics of short ≤ 5-nt misalignments. The significance of short DNA misalignments supports an important role for erroneous non-homologous and micro-homology-dependent DSB repair in mtDNA deletion formation. The consistency of the results of our analysis across species further suggests a shared mode of mtDNA deletion mutagenesis.

How Chaotic Is Genome Chaos?

Simple Summary

Cancer genomes can undergo major restructurings involving many chromosomal locations at key stages in tumor development. This restructuring process has been designated “genome chaos” by some authors. In order to examine how chaotic cancer genome restructuring may be, the cell and molecular processes for DNA restructuring are reviewed. Examination of the action of these processes in various cancers reveals a degree of specificity that indicates genome restructuring may be sufficiently reproducible to enable possible therapies that interrupt tumor progression to more lethal forms.

Abstract

Cancer genomes evolve in a punctuated manner during tumor evolution. Abrupt genome restructuring at key steps in this evolution has been called “genome chaos.” To answer whether widespread genome change is truly chaotic, this review (i) summarizes the limited number of cell and molecular systems that execute genome restructuring, (ii) describes the characteristic signatures of DNA changes that result from activity of those systems, and (iii) examines two cases where genome restructuring is determined to a significant degree by cell type or viral infection. The conclusion is that many restructured cancer genomes display sufficiently unchaotic signatures to identify the cellular systems responsible for major oncogenic transitions, thereby identifying possible targets for therapies to inhibit tumor progression to greater aggressiveness.

XRCC1 promotes replication restart, nascent fork degradation and mutagenic DNA repair in BRCA2-deficient cells

Homologous recombination/end joining (HR/HEJ)-deficient cancers with BRCA mutations utilize alternative DNA double-strand break repair pathways, particularly alternative non-homologous end joining or microhomology-mediated end joining (alt-EJ/MMEJ) during S and G2 cell cycle phases. Depletion of alt-EJ factors, including XRCC1, PARP1 and POLQ, is synthetically lethal with BRCA2 deficiency; yet, XRCC1 roles in HR-deficient cancers and replication stress are enigmatic. Here, we show that after replication stress, XRCC1 forms an active repair complex with POLQ and MRE11 that supports alt-EJ activity in vitro. BRCA2 limits XRCC1 recruitment and repair complex formation to suppress alt-EJ at stalled forks. Without BRCA2 fork protection, XRCC1 enables cells to complete DNA replication at the expense of increased genome instability by promoting MRE11-dependent fork resection and restart. High XRCC1 and MRE11 gene expression negatively impacts Kaplan-Meier survival curves and hazard ratios for HR-deficient breast cancer patients in The Cancer Genome Atlas. The additive effects of depleting both BRCA2 and XRCC1 indicate distinct pathways for replication restart. Our collective data show that XRCC1-mediated processing contributes to replication fork degradation, replication restart and chromosome aberrations in BRCA2-deficient cells, uncovering new roles of XRCC1 and microhomology-mediated repair mechanisms in HR-deficient cancers, with implications for chemotherapeutic strategies targeting POLQ and PARP activities.

Joint single cell DNA-seq and RNA-seq of gastric cancer cell lines reveals rules of in vitro evolution

Cancer cell lines are not homogeneous nor are they static in their genetic state and biological properties. Genetic, transcriptional and phenotypic diversity within cell lines contributes to the lack of experimental reproducibility frequently observed in tissue-culture-based studies. While cancer cell line heterogeneity has been generally recognized, there are no studies which quantify the number of clones that coexist within cell lines and their distinguishing characteristics. We used a single-cell DNA sequencing approach to characterize the cellular diversity within nine gastric cancer cell lines and integrated this information with single-cell RNA sequencing. Overall, we sequenced the genomes of 8824 cells, identifying between 2 and 12 clones per cell line. Using the transcriptomes of more than 28 000 single cells from the same cell lines, we independently corroborated 88% of the clonal structure determined from single cell DNA analysis. For one of these cell lines, we identified cell surface markers that distinguished two subpopulations and used flow cytometry to sort these two clones. We identified substantial proportions of replicating cells in each cell line, assigned these cells to subclones detected among the G0/G1 population and used the proportion of replicating cells per subclone as a surrogate of each subclone's growth rate.

DNA barcoding reveals that injected transgenes are predominantly processed by homologous recombination in mouse zygote

Mechanisms that ensure repair of double-strand DNA breaks (DSBs) are instrumental in the integration of foreign DNA into the genome of transgenic organisms. After pronuclear microinjection, exogenous DNA is usually found as a concatemer comprising multiple co-integrated transgene copies. Here, we investigated the contribution of various DSB repair pathways to the concatemer formation. We injected mouse zygotes with a pool of linear DNA molecules carrying unique barcodes at both ends and obtained 10 transgenic embryos with 1–300 transgene copies. Sequencing the barcodes allowed us to assign relative positions to the copies in concatemers and detect recombination events that occurred during integration. Cumulative analysis of approximately 1,000 integrated copies reveals that over 80% of them underwent recombination when their linear ends were processed by synthesis-dependent strand annealing (SDSA) or double-strand break repair (DSBR). We also observed evidence of double Holliday junction (dHJ) formation and crossing over during the concatemer formations. Sequencing indels at the junctions between copies shows that at least 10% of DNA molecules introduced into the zygotes are ligated by non-homologous end joining (NHEJ). Our barcoding approach, verified with Pacific Biosciences Single Molecule Real-Time (SMRT) long-range sequencing, documents high activity of homologous recombination after DNA microinjection.

Rad51 recruitment and exclusion of non-homologous end joining during homologous recombination at a Tus/Ter mammalian replication fork barrier

Classical non-homologous end joining (C-NHEJ) and homologous recombination (HR) compete to repair mammalian chromosomal double strand breaks (DSBs). However, C-NHEJ has no impact on HR induced by DNA nicking enzymes. In this case, the replication fork is thought to convert the DNA nick into a one-ended DSB, which lacks a readily available partner for C-NHEJ. Whether C-NHEJ competes with HR at a non-enzymatic mammalian replication fork barrier (RFB) remains unknown. We previously showed that conservative "short tract" gene conversion (STGC) induced by a chromosomal Tus/Ter RFB is a product of bidirectional replication fork stalling. This finding raises the possibility that Tus/Ter-induced STGC proceeds via a two-ended DSB intermediate. If so, Tus/Ter-induced STGC might be subject to competition by C-NHEJ. However, in contrast to the DSB response, where genetic ablation of C-NHEJ stimulates HR, we report here that Tus/Ter-induced HR is unaffected by deletion of either of two C-NHEJ genes, Xrcc4 or Ku70. These results show that Tus/Ter-induced HR does not entail the formation of a two-ended DSB to which C-NHEJ has competitive access. We found no evidence that the alternative end-joining factor, DNA polymerase θ, competes with Tus/Ter-induced HR. We used chromatin-immunoprecipitation to compare Rad51 recruitment to a Tus/Ter RFB and to a neighboring site-specific DSB. Rad51 accumulation at Tus/Ter was more intense and more sustained than at a DSB. In contrast to the DSB response, Rad51 accumulation at Tus/Ter was restricted to within a few hundred base pairs of the RFB. Taken together, these findings suggest that the major DNA structures that bind Rad51 at a Tus/Ter RFB are not conventional DSBs. We propose that Rad51 acts as an "early responder" at stalled forks, binding single stranded daughter strand gaps on the arrested lagging strand, and that Rad51-mediated fork remodeling generates HR intermediates that are incapable of Ku binding and therefore invisible to the C-NHEJ machinery.

Mechanism of tandem duplication formation in BRCA1-mutant cells

Small, approximately 10-kilobase microhomology-mediated tandem duplications are abundant in the genomes of BRCA1-linked but not BRCA2-linked breast cancer. Here we define the mechanism underlying this rearrangement signature. We show that, in primary mammalian cells, BRCA1, but not BRCA2, suppresses the formation of tandem duplications at a site-specific chromosomal replication fork barrier imposed by the binding of Tus proteins to an array of Ter sites. BRCA1 has no equivalent role at chromosomal double-stranded DNA breaks, indicating that tandem duplications form specifically at stalled forks. Tandem duplications in BRCA1 mutant cells arise by a replication restart-bypass mechanism terminated by end joining or by microhomology-mediated template switching, the latter forming complex tandem duplication breakpoints. Solitary DNA ends form directly at Tus-Ter, implicating misrepair of these lesions in tandem duplication formation. Furthermore, BRCA1 inactivation is strongly associated with ~10 kilobase tandem duplications in ovarian cancer. This tandem duplicator phenotype may be a general signature of BRCA1-deficient cancer.

FANCJ helicase controls the balance between short- and long-tract gene conversions between sister chromatids

The FANCJ DNA helicase is linked to hereditary breast and ovarian cancers as well as bone marrow failure disorder Fanconi anemia (FA). Although FANCJ has been implicated in the repair of DNA double-strand breaks (DSBs) by homologous recombination (HR), the molecular mechanism underlying the tumor suppressor functions of FANCJ remains obscure. Here, we demonstrate that FANCJ deficient human and hamster cells exhibit reduction in the overall gene conversions in response to a site-specific chromosomal DSB induced by I-SceI endonuclease. Strikingly, the gene conversion events were biased in favour of long-tract gene conversions in FANCJ depleted cells. The fine regulation of short- (STGC) and long-tract gene conversions (LTGC) by FANCJ was dependent on its interaction with BRCA1 tumor suppressor. Notably, helicase activity of FANCJ was essential for controlling the overall HR and in terminating the extended repair synthesis during sister chromatid recombination (SCR). Moreover, cells expressing FANCJ pathological mutants exhibited defective SCR with an increased frequency of LTGC. These data unravel the novel function of FANCJ helicase in regulating SCR and SCR associated gene amplification/duplications and imply that these functions of FANCJ are crucial for the genome maintenance and tumor suppression.

Endogenous sequence patterns predispose the repair modes of CRISPR/Cas9-induced DNA double-stranded breaks in Arabidopsis thaliana

The possibility to predict the outcome of targeted DNA double-stranded break (DSB) repair would be desirable for genome editing. Furthermore the consequences of mis-repair of potentially cell-lethal DSBs and the underlying pathways are not yet fully understood. Here we study the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-induced mutation spectra at three selected endogenous loci in Arabidopsis thaliana by deep sequencing of long amplicon libraries. Notably, we found sequence-dependent genomic features that affected the DNA repair outcome. Deletions of 1-bp to <1000-bp size and/or very short insertions, deletions >1 kbp (all due to NHEJ) and deletions combined with insertions between 5-bp to >100 bp [caused by a synthesis-dependent strand annealing (SDSA)-like mechanism] occurred most frequently at all three loci. The appearance of single-stranded annealing events depends on the presence and distance between repeats flanking the DSB. The frequency and size of insertions is increased if a sequence with high similarity to the target site was available in cis. Most deletions were linked to pre-existing microhomology. Deletion and/or insertion mutations were blunt-end ligated or via de novo generated microhomology. While most mutation types and, to some degree, their predictability are comparable with animal systems, the broad range of deletion mutations seems to be a peculiar feature of the plant A. thaliana.

FOXM1 expression is significantly associated with chemotherapy resistance and adverse prognosis in non-serous epithelial ovarian cancer patients

BACKGROUND

Epithelial ovarian cancer (EOC) is a spectrum of different diseases, which makes their treatment a challenge. Forkhead box M1 (FOXM1) is an oncogene aberrantly expressed in many solid cancers including serous EOC, but its role in non-serous EOCs remains undefined. We examined FOXM1 expression and its correlation to prognosis across the three major EOC subtypes, and its role in tumorigenesis and chemo-resistance in vitro.

METHODS

Gene signatures were generated by microarray for 14 clear-cell and 26 endometrioid EOCs, and 15 normal endometrium snap-frozen biopsies. Validation of FOXM1 expression was performed by RT-qPCR and immunohistochemistry in the same samples and additionally in 50 high-grade serous EOCs and in their most adequate normal controls (10 luminal fallopian tube and 20 ovarian surface epithelial brushings). Correlations of FOXM1 expression to clinic-pathological parameters and patients' prognosis were evaluated by Kaplan-Meier and Cox proportional-hazards analyses. OVCAR-3 and two novel deeply characterized EOC cell lines (EOC-CC1 and OSPC2, with clear-cell and serous subtype, respectively) were employed for in vitro studies. Effects of FOXM1 inhibition by transient siRNA transfection were evaluated on cell-proliferation, cell-cycle, colony formation, invasion, and response to conventional first- and second-line anticancer agents, and to the PARP-inhibitor olaparib. Gene signatures of FOXM1-silenced cell lines were generated by microarray and confirmed by RT-qPCR.

RESULTS

A significant FOXM1 mRNA up-regulation was found in EOCs compared to normal controls. FOXM1 protein overexpression significantly correlated to serous histology (p = 0.001) and advanced FIGO stage (p = 0.004). Multivariate analyses confirmed FOXM1 protein overexpression as an independent indicator of worse disease specific survival in non-serous EOCs, and of shorter time to progression in platinum-resistant cases. FOXM1 downregulation in EOC cell lines inhibited cell growth and clonogenicity, and promoted the cytotoxic effects of platinum compounds, doxorubicin hydrochloride and olaparib. Upon FOXM1 knock-down in EOC-CC1 and OSPC2 cells, microarray and RT-qPCR analyses revealed the deregulation of several common and other unique subtype-specific FOXM1 putative targets involved in cell cycle, metastasis, DNA repair and drug response.

CONCLUSIONS

FOXM1 is up-regulated in all three major EOCs subtypes, and is a prognostic biomarker and a potential combinatorial therapeutic target in platinum resistant disease, irrespective of tumor histology.

Genomic rearrangements induced by unscheduled DNA double strand breaks in somatic mammalian cells

DNA double-strand breaks (DSBs) are highly toxic lesions that can lead to profound genome rearrangements and/or cell death. They routinely occur in genomes due to endogenous or exogenous stresses. Efficient repair systems, canonical non-homologous end-joining and homologous recombination exist in the cell and not only ensure the maintenance of genome integrity but also, via specific programmed DNA double-strand breaks, permit its diversity and plasticity. However, these repair systems need to be tightly controlled because they can also generate genomic rearrangements. Thus, when DSB repair is not properly regulated, genome integrity is no longer guaranteed. In this review, we will focus on non-programmed genome rearrangements generated by DSB repair, in somatic cells. We first discuss genome rearrangements induced by homologous recombination and end-joining. We then discuss recently described rearrangement mechanisms, driven by microhomologies, that do not involve the joining of DNA ends but rather initiate DNA synthesis (microhomology-mediated break-induced replication, fork stalling and template switching and microhomology-mediated template switching). Finally, we discuss chromothripsis, which is the shattering of a localized region of the genome followed by erratic rejoining.

Correction: Complex Breakpoints and Template Switching Associated with Non-canonical Termination of Homologous Recombination in Mammalian Cells

[This corrects the article DOI: 10.1371/journal.pgen.1006410.].