PLoS ONE
Home Tyrosyl-DNA phosphodiesterases are involved in mutagenic events at a ribonucleotide embedded into DNA in human cells
Tyrosyl-DNA phosphodiesterases are involved in mutagenic events at a ribonucleotide embedded into DNA in human cells
Tyrosyl-DNA phosphodiesterases are involved in mutagenic events at a ribonucleotide embedded into DNA in human cells

Competing Interests: The authors have declared that no competing interests exist.

Article Type: research-article Article History
Abstract

Ribonucleoside triphosphates are often incorporated into genomic DNA during DNA replication. The accumulation of unrepaired ribonucleotides is associated with genomic instability, which is mediated by DNA topoisomerase 1 (Top1) processing of embedded ribonucleotides. The cleavage initiated by Top1 at the site of a ribonucleotide leads to the formation of a Top1-DNA cleavage complex (Top1cc), occasionally resulting in a DNA double-strand break (DSB). In humans, tyrosyl-DNA phosphodiesterases (TDPs) are essential repair enzymes that resolve the trapped Top1cc followed by downstream repair factors. However, there is limited cellular evidence of the involvement of TDPs in the processing of incorporated ribonucleotides in mammals. We assessed the role of TDPs in mutagenesis induced by a single ribonucleotide embedded into DNA. A supF shuttle vector site-specifically containing a single riboguanosine (rG) was introduced into the human lymphoblastoid TK6 cell line and its TDP1-, TDP2-, and TDP1/TDP2-deficient derivatives. TDP1 and TDP2 insufficiency remarkably decreased the mutant frequency caused by an embedded rG. The ratio of large deletion mutations induced by rG was also substantially lower in TDP1/TDP2-deficient cells than wild-type cells. Furthermore, the disruption of TDPs reduced the length of rG-mediated large deletion mutations. The recovery ratio of the propagated plasmid was also increased in TDP1/TDP2-deficient cells after the transfection of the shuttle vector containing rG. The results suggest that TDPs-mediated ribonucleotide processing cascade leads to unfavorable consequences, whereas in the absence of these repair factors, a more error-free processing pathway might function to suppress the ribonucleotide-induced mutagenesis. Furthermore, base substitution mutations at sites outside the position of rG were detected in the supF gene via a TDPs-independent mechanism. Overall, we provide new insights into the mechanism of mutagenesis induced by an embedded ribonucleotide in mammalian cells, which may lead to the fatal phenotype in the ribonucleotide excision repair deficiency.

Takeishi,Kogashi,Odagiri,Sasanuma,Takeda,Yasui,Honma,Suzuki,Kamiya,Sugasawa,Ura,Sassa,and Sobol: Tyrosyl-DNA phosphodiesterases are involved in mutagenic events at a ribonucleotide embedded into DNA in human cells

Introduction

In eukaryotic cells, ribonucleoside triphosphates are occasionally incorporated by DNA polymerases (Pols) during DNA replication [1]. It is estimated that more than 1,000,000 ribonucleotide molecules are embedded into genomic DNA during single-round replication in mammals [2]. The transient incorporation of ribonucleotides by Pols is beneficial for efficient DNA repair. For example, genomic ribonucleotides serve as a signal for strand discrimination during DNA mismatch repair [3,4]. Ribonucleotide incorporation during DNA double-strand break (DSB) repair also promotes efficient DNA ligation reactions [5]. However, dysfunction of the appropriate repair pathway for embedded ribonucleotides causes the accumulation of ribonucleotides, leading to genomic instability and resulting in tumorigenesis and disease [68]. Therefore, genomic ribonucleotide incorporation is a double-edged sword for genome integrity.

Embedded ribonucleotides are primarily repaired by RNase H2-initiated ribonucleotide excision repair (RER). During this process, RNase H2 excises the 5′-side of the incorporated ribonucleotide into DNA [9,10]. Then, Pols δ and ε with PCNA displace the DNA strand with a ribonucleotide, followed by the excision of the flapped DNA by flap endonuclease FEN1 and nick sealing by DNA ligases [11]. RNase H2 enzyme inactivation leads to severe cellular abnormalities, including persistent DNA damage response and p53-mediated apoptosis, which results in embryonic lethality in mice [2,12]. This fatal phenotype can be caused by the alternate false-processing of a ribonucleotide initiated by DNA topoisomerase 1 (Top1). The Top1 cleaves the 3′ of the embedded ribonucleotide and forms a covalent Top1-DNA cleavage complex (Top1cc) at the 3′-terminus [13,14]. The phosphate of Top1cc is susceptible to nucleophilic attack of the 2′-hydroxyl of the ribonucleotide sugar backbone, generating a 2′,3′-cyclic phosphate end and releasing Top1. Then, Top1 cleaves in the 5′-upstream region of the nick containing the 2′,3′-cyclic phosphate, forming a second Top1cc and DNA gap. In yeast, this Top1 processing of a ribonucleotide causes short deletion mutations at tandem repeat sequences [14]. However, such mutations have not been observed in the genomic DNA of RNase H2-deficient mouse-derived cells [7]. The trapped Top1cc also causes DNA replication collapse and DSBs formation [15]. Importantly, it has been reported that Top1-mediated processing of ribonucleotides induces DSBs in yeast and mammalian cells [16,17]. In such cases, the trapped Top1cc with a DSB must be repaired by the downstream repair factors.

In eukaryotes, trapped topoisomerase-DNA complexes are resolved by proteasomal degradation, followed by hydrolysis of the phosphotyrosyl bond between topoisomerases and the DNA backbone by tyrosyl DNA phosphodiesterases (TDPs). TDP1 hydrolyzes the phosphodiester bonds at the 3′-DNA terminus linked to tyrosyl moiety within the trapped Top1cc [18]. In addition, TDP1 has broad substrate specificity for various 3′-blocking lesions resulting from different types of DNA damage and chain-terminating anticancer nucleosides [1921]. Thus, the biological and clinical importance of the function of TDP1 has been well investigated [22]. TDP2 can resolve 5′-phosphotyrosyl bonds between Top2 and the DNA strand break terminus in Top2cc [23]. In addition to its canonical function, TDP2 also possesses weak cleavage activity for 3′-phosphotyrosyl bonds [24], suggesting the overlapping roles of TDP1 and TDP2 in Top1cc repair [25]. TDP1 and TDP2 are also involved in the repair of topoisomerase-induced DSBs. For example, TDP1 facilitates non-homologous end-joining (NHEJ) pathway via functional interactions with NHEJ factors Ku70/80 and XLF in human cells [26]. In addition, TDP2 is required for efficient NHEJ upon the formation of Top2-mediated DSBs [27].

Based on the biochemical and biological functions of TDP1 and TDP2, these enzymes appear to affect ribonucleotide-induced genomic instability through the processing of Top1-DNA complexes. Based on an in vitro study, TDP1 excises the 3′-blocking topoisomerase-derived lesions linked to a single embedded ribonucleotide [28]. However, the disruption of TDP1 in RNase H2-deficient human cells does not exhibit the obvious effect based on growth and survival rates [29]. There is limited evidence supporting the role of TDP1 and TDP2 in mutagenesis of an embedded ribonucleotide in mammalian cells.

We previously reported that a single ribonucleotide induced deletion mutations in the simian virus 40 (SV40)-based shuttle vector that can be propagated in human cells [30]. SV40 DNA replication is carried out by human replication factors including PCNA, Top1, and DNA polymerase δ [31]. Unlike the genome in which DNA replication timing is strictly regulated, the shuttle vector containing SV40 replication origin is rapidly propagated in the presence of SV40 large-T antigen in mammalian cells [32]; the copy number of the vector can reach 100,000 copies per cell within 48 h after transfection. Due to its high replication rate, it is expected that DNA lesions incorporated into the vector could escape from the canonical repair, thereby giving rise to mutation frequencies in the range of 10−3 [33,34]. Using the same experimental technique that quantifies the mutagenic properties of an embedded ribonucleotide in cells, we assessed the effect of TDPs deficiency on the mutagenic potential of a single riboguanosine (rG) incorporated into DNA. The disruption of TDP1 and TDP2 markedly suppressed the mutagenic potential of an embedded rG. Our results provide in cellulo evidence supporting that TDPs mediate mutagenic processing of a ribonucleotide incorporated into DNA rather than protect the DNA against it in mammalian cells, which might contribute to the global genomic instability in RER deficiency. A possible mechanism for the involvement of DNA repair pathways in ribonucleotide-induced mutagenesis is discussed.

Materials and methods

Cell culture

The human TSCER2 cell line, derived from human lymphoblastoid TK6 [35], and its derivatives were cultured in RPMI-1640 medium (Nacalai Tesque) supplemented with 10% fetal bovine serum, 200 μg/mL sodium pyruvate, 100 U/mL penicillin, and 100 μg/mL streptomycin at 37°C in an atmosphere of 5% CO2 and 100% humidity. TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells were established as previously described [36]. The loss of TDP1 or TDP2 expression in cell lines was confirmed in parallel in a previous study [37].

Construction of a closed circular double-stranded plasmid DNA site-specifically containing a single ribonucleotide

A plasmid containing a single rG at position 160 of pMY189, the numbering used is that of the supF reporter gene in the plasmid with denoting the EcoRI site before supF, was constructed as previously described [38]. In detail, the 5′-phosphorylated 21-mer primer (5′-CGACTTCGAAGXTTCGAATCC-3′, where X represents dG or rG) was obtained from Tsukuba Oligo Service Co. Ltd. (Ibaraki, Japan). The primer (16 pmol) was annealed with single-stranded DNA (5 μg) of pMY189 that was isolated from JM109 Escherichia coli (E. coli) strain harboring the plasmid using a VCSM13 helper phage. Then, the annealed DNA was incubated with 2 units of T4 DNA polymerase (Sigma) and 250 units of T4 DNA ligase (Takara) at 37°C for 3 h in the reaction mixture in a final volume of 100 μl contained 50 mM Tris-HCl (pH 7.9), 50 mM NaCl, 7 mM MgCl2, 1 mM dithiothreitol, 1 mM ATP, 0.5 mM dNTPs, and 12 μg of T4 gene 32 protein (Nippon Gene). The reaction product was treated with 16 units of dam methyltransferase (New England Biolabs) to restore the bacterial methylation pattern. The closed circular double-stranded DNA with dG or rG was purified using ultracentrifugation in a cesium chloride-ethidium bromide gradient. Homogeneity and no degradation of the purified DNA was confirmed by agarose gel electrophoresis before the transfection experiment.

Forward mutation assay in the supF gene

The closed circular double-stranded DNA with or without rG (500 ng) was transfected into TSCER2 cells (2 × 106 cells) using NEPA21 electroporator (Nepa Gene Co. Ltd.) following the manufacturer’s instructions. After incubation for 48 h at 37°C, the replicated plasmid DNA was extracted from cells and treated with DpnI to digest unreplicated plasmids. The resultant DNA was introduced into the KS40/pOF105 indicator E. coli strain [39]. The transformed cells were plated onto Luria-Bertani (LB) agar selection plates containing nalidixic acid (50 μg/ml), streptomycin (100 μg/ml), ampicillin (150 μg/ml), chloramphenicol (30 μg/ml), 5-bromo-4-chloro-3-indolyl-β-D-galactoside (80 μg/ml), and isopropyl β-D-thiogalactopyranoside (23.8 μg/ml) to select the E. coli cells harboring a mutated supF gene. Then, the cells were plated onto LB titer plates containing ampicillin (150 μg/ml) and chloramphenicol (30 μg/ml) to determine the total number of transformants. The mutant frequency of supF was determined based on the numbers of white and faint blue colonies on selection plates and the number of colonies on titer plates. For sequencing analysis to determine the mutation spectrum, the supF gene was amplified from mutant colonies by PCR with primers 5′-AGTGCCACCTGACATCTA-3′ and 5′-CAGCAGATTACGCGCAGA-3′. The resulting fragment was treated with ExoSapIT (Thermo Fisher Scientific) to degrade excess primers and deoxynucleoside triphosphates and sequenced using BigDye terminator v3.1 (Thermo Fisher Scientific) on an ABI 3130 genetic analyzer (Thermo Fisher Scientific).

Determining the recovery ratio of propagated plasmids from human cells

An internal control vector pST189K was constructed by replacing the ampicillin resistance gene and M13 intergenic sequence of pZ189-StuI with the kanamycin resistance gene [40]. The synthetic closed circular DNA containing dG or rG (500 ng) and an equal amount of pST189K (500 ng) were co-transfected into TSCER2 cells (2 × 106 cells) by using NEPA21 electroporator (Nepa Gene Co. Ltd.). After 48 h of incubation, the replicated plasmids were extracted from cells, digested with DpnI, and then introduced into the KS40/pOF105 strain as described above. The transformed cells were plated onto LB plates containing ampicillin (150 μg/ml) or kanamycin (50 μg/ml) in the presence of chloramphenicol (30 μg/ml). The recovery ratio of replicated plasmids was indicated by the number of ampicillin-resistant colonies relative to the number of kanamycin-resistant colonies.

Statistical analysis

The statistical significance of the supF mutant frequency was evaluated using Dunnett’s multiple comparison test. The mutation spectra between the wild-type, TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells were statistically analyzed using Fisher’s exact test with multiple testing corrections using Holm’s method. Statistical significance was indicated by P-values of <0.05.

Results

Mutagenic events induced by a ribonucleotide embedded into a plasmid in human cells

To investigate the mutation induced by a ribonucleotide embedded into DNA, we constructed the shuttle vector site-specifically containing rG and introduced into human TSCER2 cells. The propagated plasmids were extracted from cells, and further introduced into E. coli KS40/pOF105 indicator strain to determine the supF mutant frequency and spectrum. Mutant frequencies were determined as the number of mutant colonies divided by the total number of colonies recovered (S1 Table). In the wild-type cells, the mutant frequency of the plasmid containing rG ((3.4 ± 0.49) × 10−3) was significantly higher than that of the control plasmid ((0.64 ± 0.044) × 10−3) (P < 0.01) (Fig 1A and 1B). Based on detailed DNA sequencing analysis in the supF gene, the most predominant mutations induced by rG were large deletion mutations ranging from 9 to 304 bp in length in the wild-type cells (Table 1 and S2 Table), which was consistent with previous observations [30]. The majority of large deletions were extended to the 5′-upstream region from the position of rG. The second most common mutations were “untargeted” mutations, i.e., base substitution mutations outside the site of rG, which was preferentially observed at guanine bases (Table 1 and S3 Table). In contrast, large deletion mutations were infrequently observed in the control plasmid.

Frequency of supF mutants in a shuttle vector containing (A) dG or (B) rG replicated in wild-type (WT), TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells. Mutant frequencies were calculated from the data of S1 Table. Data for control dG are expressed as the mean ± standard error (SE) of at least three independent experiments from two independent clones. Data for rG are expressed as the mean ± SE of three to six independent experiments from two independent clones. *Significant difference between the assessed and WT cells; P < 0.05.
Fig 1

Frequency of supF mutants in a shuttle vector containing (A) dG or (B) rG replicated in wild-type (WT), TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells. Mutant frequencies were calculated from the data of S1 Table. Data for control dG are expressed as the mean ± standard error (SE) of at least three independent experiments from two independent clones. Data for rG are expressed as the mean ± SE of three to six independent experiments from two independent clones. *Significant difference between the assessed and WT cells; P < 0.05.

Table 1
Mutation spectrum in supF induced by a single ribonucleotide in wild-type, TDP1-/-,TDP2-/-, and TDP1-/-/TDP2-/- cells.
0GaMutationWild-typeTDP1-/-TDP2-/-TDP1-/-/TDP2-/-
dGMutations at site of 0G
G:C → T:A6 (15)8 (26)4 (10)4 (12)
G:C → A:T0 (0)1 (3.0)1 (3.0)2 (6.0)
Large deletion3 (7.5)3 (10)5 (13)1 (3.0)
Base substitution, 1–2 bp deletion or insertion outside 0G31 (78)19 (61)28 (74)27 (79)
Total40 (100)31 (100)38 (100)34 (100)
rGMutations at site of 0G
G:C → T:A1 (1.7)0 (0)2 (3.0)3 (5.0)
G:C → C:G1 (1.7)2 (3.0)3 (4.5)0 (0)
G:C → A:T1 (1.7)0 (0)1 (1.5)0 (0)
Large deletion33 (55)21 (35)27 (41)15 (24)**
Base substitution, 1–2 bp deletion or insertion outside 0G23 (38)37 (62)29 (45)44 (71)
Other1 (2.0)0 (0)3 (5.0)0 (0)
Total60 (100)60 (100)65 (100)62 (100)

a 0G indicates the 160th base position in pMY189 where dG or rG was introduced.

b Numbers in parentheses represent the percentage of total number of mutants.

** Significant difference between the assessed cells and WT cells; P < 0.01.

Effect of the disruption of TDP1 and TDP2 on the mutagenic potential of an embedded ribonucleotide

We examined the impact of TDP1 and TDP2 deficiency on the mutation induced by a ribonucleotide incorporated into DNA. To compare the supF mutant frequency between cell lines, it should be noted that the population doubling times of TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells (14 ± 1.6, 13 ± 0.10, and 14 ± 1.2 h, respectively) were similar to that of the wild-type cells (14 ± 0.66 h). Although the mutant frequency induced by rG was slightly decreased in the absence of TDP1 ((2.2 ± 0.28) × 10−3) compared with the wild-type cells ((3.4 ± 0.49) × 10−3), this difference was not statistically significant (Fig 1B). The mutant frequencies in supF gene in the control ((0.90 ± 0.11) × 10−3) and rG-containing plasmid ((3.5 ± 0.089) × 10−3) in TDP2-/- cells were comparable to those in the wild-type cells ((0.64 ± 0.044) × 10−3 and (3.4 ± 0.49) × 10−3, respectively) (Fig 1A and 1B). TDP2 possesses a weak 3′-tyrosyl DNA phosphodiesterase activity, which contributes to resolving Top1-mediated DNA lesions in the absence of TDP1 [24]. Thus, we suspected that the disruption of both TDPs leads to the marked modulation of the mutagenic potential of rG embedded into DNA. For results, the frequency of mutants induced by rG was significantly lower in TDP1-/-/TDP2-/- cells ((1.8 ± 0.26) × 10−3) than in the wild-type cells ((3.4 ± 0.49) × 10−3) (Fig 1B).

Deletion mutations caused by an embedded ribonucleotide are suppressed in the absence of TDP1 and TDP2

We also assessed whether the loss of TDPs modulated the mutation spectrum induced by a ribonucleotide in supF gene. Among the mutations observed, we focused on large deletion and untargeted mutations that are predominantly induced by rG. In TDP1-/- cells, the most predominant mutations induced by rG were large deletion mutations as observed in the wild-type cells (Table 1): the ratio of large deletion mutations in TDP1-/- cells (35%) was lower than that in the wild-type cells (55%), although the difference was not significant. Based on the mutant frequency and mutation spectra given in Fig 1B and Table 1, the normalized mutant frequency derived from rG-induced large deletion mutations was calculated to be 0.77 × 10−3 in TDP1-/- cells, which was approximately 2.4-fold lower than that of the wild-type cells (1.9 × 10−3) (Fig 2A). However, the frequency of untargeted mutations was comparable between the wild-type and TDP1-/- cells (1.3 × 10−3 and 1.4 × 10−3, respectively) (Fig 2B). Regardless of the status of TDP2, the mutation spectrum in the supF gene was not substantially different (Table 1). The normalized frequency of large deletion mutations was slightly decreased by the disruption of TDP2 (1.4 × 10−3) when compared with the wild-type cells (1.9 × 10−3) (Fig 2A).

The normalized mutant frequencies of (A) large deletion and (B) untargeted mutations detected in the supF gene in the plasmids replicated in the wild-type, TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells. The frequency values of each mutant were calculated based on the number of mutants represented in Table 1. A significant difference was determined based on the number of mutations shown in Table 1.
Fig 2

The normalized mutant frequencies of (A) large deletion and (B) untargeted mutations detected in the supF gene in the plasmids replicated in the wild-type, TDP1-/-, TDP2-/-, and TDP1-/-/TDP2-/- cells. The frequency values of each mutant were calculated based on the number of mutants represented in Table 1. A significant difference was determined based on the number of mutations shown in Table 1.

In TDP1-/-/TDP2-/- cells, the ratio of rG-mediated large deletion mutations (24%) was significantly lower than that in the wild-type cells (55%) (Table 1), according to the decreased mutant frequency. The frequency derived from large deletion mutations in TDP1-/-/TDP2-/- cells (0.43 × 10−3) was 4.3-fold lower than that in the wild-type cells (1.9 × 10−3) (Fig 2A). On the other hand, the frequency of untargeted mutations was comparable between the wild-type and TDP1-/-/TDP2-/- cells (1.3 × 10−3 and 1.3 × 10−3, respectively) (Fig 2B), suggesting that such untargeted mutations related to rG were not associated with the function of TDP1 or TDP2. Notably, there was no significant difference between the TDP1-/-/TDP2-/- and wild-type cells in the mutant frequency and mutation spectrum in the supF gene of the control plasmid (Fig 1A and Table 1).

Tyrosyl-DNA phosphodiesterases deficiency alters the length of large deletion mutations

To examine whether the loss of both TDP1 and TDP2 affects the magnitude of deletions induced by a ribonucleotide, we compared the length of rG-induced large deletion mutations between cells. The median length of large deletion mutations was determined to be 25 bp in TDP1-/-/TDP2-/- cells, which was remarkably shorter than that observed in the wild-type cells (145 bp; Fig 3A). These results raise a question regarding the loss of TDPs affecting the recovery of replicated plasmids from cells resulting from the reduced frequency and length of large deletion mutations. Thus, we compared the relative recovery ratio of propagated plasmids between wild-type and TDP1-/-/TDP2-/- cells. To determine the recovery ratio, the synthetic vector containing dG or rG was co-transfected with an equal amount of the internal control vector pST189K harboring kanamycin-resistant gene. After 48 h incubation, the replicated plasmids were extracted and introduced into the indicator E. coli strain, which was then plated onto LB plates. The relative recovery ratio was calculated from the numbers of ampicillin-resistant colonies relative to kanamycin-resistant colonies. The recovery ratio of replicated plasmids after introducing the shuttle vector containing rG (0.55 ± 0.087) was lower than that related to dG (1.8 ± 0.10) in the wild-type cells. Interestingly, we found that the recovery ratio after the transfection of the shuttle vector with rG increased with disruption of TDPs (0.86 ± 0.11) (Fig 3B), suggesting that TDPs-dependent processing of an embedded ribonucleotide is mutagenic and compromises the plasmid propagation.

Effect of tyrosyl-DNA phosphodiesterase deficiency on the length of large deletion mutations and the recovery ratio of propagated plasmids.
Fig 3

Effect of tyrosyl-DNA phosphodiesterase deficiency on the length of large deletion mutations and the recovery ratio of propagated plasmids.

(A) Box plot representing the length of large deletion mutations observed in the indicated genotypes. A significant difference is indicated by an asterisk (*P <0.05). (B) The synthetic closed-circular DNA with or without rG was co-transfected with pST189K into cells. After 48 h of incubation, the replicated plasmids were recovered and introduced into the KS40/pOF105 strain. The transformed E. coli cells were plated onto LB plates with ampicillin (150 μg/mL) or kanamycin (50 μg/mL) in the presence of chloramphenicol (30 μg/mL). The numbers of ampicillin-resistant colonies relative to those of kanamycin-resistant colonies are presented. Values presented are the mean ± SE from six independent cultures derived from two transfections. A significant difference is indicated by an asterisk (*P <0.05).

Discussion

The impact of ribonucleotide accumulation into the genome has been investigated in vitro and in vivo [1,41,42], and the false processing of incorporated ribonucleotides has been found to lead to genomic instability [14,16,17]. Despite extensive studies on the role of TDPs in genome stability, little is known about the involvement of TDPs in the processing of the ribonucleotide incorporated into DNA in mammalian cells. In principle, disruption of the canonical RER is considered to be required to assess the mutagenic processing initiated by Top1. RER-deficient mammalian cells, however, are viable only in the absence of the p53 tumor suppressor gene. The loss of p53 causes another problem whereby the background mutation rate is dramatically increased because of the roles of p53 in genome integrity [43], which could complicate the interpretation of results [44]. To specifically visualize the mutagenic events at an embedded ribonucleotide that is readily removed by RER, we employed the shuttle vector system that is rapidly propagated in cells, in which some extent of an embedded ribonucleotide might escape the canonical repair in RER-proficient TK6 cells. Based on our results, a single embedded rG induced predominantly large deletion mutations in which the deleted regions were extended to the 5′-direction from the site of rG in the supF gene (S2 Table). In cells, the mutagenic processing of an embedded ribonucleotide is initiated by the cleavage at the 3′-side of a ribonucleotide, generating a 3′-dirty DNA end containing 2′,3′-cyclic phosphate through the formation of Top1cc. Then, the second Top1cc is formed at 5′-upstream of the nick, generating the gap [45]. Because Top1cc is also associated with the formation of DSBs upon DNA replication in cells [16,17], the deletion mutations observed in this study can be the result of these Top1-mediated fault processing events.

TDP1 has an important role of tolerating DNA damage resulting from the formation of Top1cc [28]. When investigating the frequency of each type of mutations, the TDP1 deficiency caused a 2.4-fold decrease in the frequency of large deletion mutations (Fig 2A). This moderate effect observed in TDP1-/- cells is possibly explained by the contribution of TDP2 activity to resolve Top1-induced DNA damage in the absence of TDP1 [24]. In addition, the disruption of both TDP1 and TDP2 resulted in a significant reduction of the rG-induced large deletion mutations (Fig 1B and Table 1). Thus, our results suggest that TDPs mediate the mutagenic processing of a ribonucleotide incorporated into DNA.

Importantly, TDP1 and TDP2 deficiencies decreased the length of rG-mediated deletion mutations (Fig 3A). In detail, most rG-induced large deletion mutations (27/33) were longer than 50 bp, and the remaining (6/33) were smaller deletions (<50 bp) in the wild-type cells (S2 Table). A fraction of >50 bp deletion mutations (5/15) decreased in TDP1-/-/TDP2-/- cells, whereas smaller deletions ranging from 3 to 50 bp were still evident (10/15) (S2 Table). Consistent with the reduced length of the deletion mutations, the recovery ratio of propagated plasmids related to rG was enhanced in TDP1-/-/TDP2-/- cells (Fig 3B). These results imply that the distinct repair pathway can play a complementary role in the TDPs-dependent mutagenic pathway in cells. Recently, it was reported that human AP endonuclease 2 (APE2) catalyzes 3′-5′ nucleolytic resection at the Top1cc terminus and 2′,3′-cyclic phosphate-DNA end [29]. The exonuclease activity of APE2 is essential for the activation of ATR-Chk1 DNA damage-response [46]. The ATR-Chk1 pathway regulates DNA replication fork speed and prevents further DNA damage [47,48]. Therefore, the APE2-mediated processing cascade can be crucial for suppression of the ribonucleotide-induced mutagenesis. Obviously, the loss of APE2 or ATR results in a synthetic lethal phenotype in RNase H2 deficiency, but cells lacking both TDP1 and RNase H2 are viable [29,49], suggesting that APE2 plays a pivotal role in protecting the genome in the absence of the canonical RER. Consistent with these studies, the yeast homologue of APE2 suppresses Top1cc-mediated mutations in RER-deficient cells [50]. It is noted that a substantial number of rG-induced smaller deletions has been observed in the absence of TDPs. These mutations might be mediated by a distinct machinery such as XPF-ERCC1 that is involved in the repair of Top1-induced DNA damage via formation of DSB [51].

In contrast, it has been proposed that TDPs are involved in repairing stalled topoisomerase-DNA complexes in conjunction with the NHEJ pathway in mammalian cells [26,27]. Upon the formation of Top1-induced DSBs, TDP1 is phosphorylated by ATM and DNA-dependent protein kinase (DNA-PK) [52]. TDP1 also promotes DNA binding by the initial NHEJ components Ku70/80 and XLF and stimulates the activity of DNA-PK, facilitating the efficient NHEJ-mediated DSB repair [26]. Consistently, the depletion of TDP1 reduces NHEJ efficiency in cells [53]. An earlier study reported that the lack of the NHEJ pathway confers resistance to camptothecin, a selective Top1 inhibitor that stabilizes Top1cc, in chicken DT40 cells [15]. Thus, in a similar manner, the disruption of the TDPs-dependent NHEJ pathway might alleviate the mutagenic consequences of ribonucleotide processing through an alternate repair such as the APE2-dependent pathway. Notably, a question arises about how large deletion mutations occur as the consequence of the TDP-dependent pathway. Although its detail mechanism remains unclear, a recent study reported that similar large deletion mutations were induced by the trapped Top1cc in yeast [54]. Based on our results, most deletion mutations were extended to the 5′-upstream region from the site of an embedded ribonucleotide (S2 Table). Because the frequency of these mutations was remarkably decreased by the depletion of TDPs (Fig 2A), such mutations might have been induced via the formation of Top1cc 5′-upstream of the ribonucleotide. Based on in vitro studies, however, the second Top1cc is formed only a few bases upstream from the nick [14,28,45]. The second Top1-cleavage could possibly occur over a range of longer nucleotides in cellular circumstances, resulting in large deletion mutations [54]. Alternatively, the cleavage of 3′-blocking lesions by TDPs may be followed by further resection reaction in certain situations.

By taking together our data and those of other studies, we proposed a model of ribonucleotide-induced mutagenesis in mammalian cells (Fig 4). The formation of 2′,3′-cyclic phosphate by Top1-mediated cleavage of a ribonucleotide is followed by a second Top1cc formation, leaving a DNA gap between the 3′-Top1cc and 5′-hydroxyl termini. The Top1cc-induced DSB via replication fork collapse can be repaired by the TDPs-dependent NHEJ pathway. Then, deletion mutations might be formed as a result of the unfavorable DSB repair. Alternatively, the other distinct factor, APE2, is involved in the resection of a trapped Top1cc-derived adduct and the 2′,3′-cyclic phosphate-DNA end, which may suppress the mutagenic consequences of ribonucleotide processing via the activation of ATR-Chk1 pathway.

Predicted model of mutagenic events induced by a ribonucleotide incorporated into DNA.
Fig 4

Predicted model of mutagenic events induced by a ribonucleotide incorporated into DNA.

The incorporated ribonucleotide is primarily repaired by RNase H2-initiated RER in an error-free manner. Alternate processing of the embedded ribonucleotide is initiated by Topoisomerase 1 (Top1)-mediated cleavage at the 3′-side of the ribonucleotide to generate a transient Top1-DNA cleavage complex (Top1cc). Then, the nucleophilic attack of the 2′-hydroxyl group to the phosphotyrosyl bond can release Top1 with the formation of 2′,3′-cyclic phosphate (depicted as a blue triangle). This 3′-blocking lesion can be processed by possible pathways as described below. A second Top1cc formation leads to the release of a DNA fragment containing 2′,3′-cyclic phosphate and the formation of a DNA gap. During the excision of a second Top1cc by TDP1 and TDP2, DNA replication fork stalling is induced by the Top1cc adduct, resulting in the formation of DSBs. The following NHEJ repair pathway might result in the formation of deletion mutations. Alternatively, the DNA terminus containing 2′,3′-cyclic phosphate and a second Top1cc adduct is subjected to resection catalyzed by the exonuclease, AP endonuclease 2 [29], which might be performed in an more error-free manner via the activation of ATR-Chk1 pathway. Untargeted base substitution mutations are induced by a ribonucleotide embedded into DNA through the TDP-independent mechanism.

Interestingly, a substantial amount of untargeted mutations were detected in the supF gene derived from plasmids containing rG in all cell lines (Fig 2B and Table 1). It was unexpected to observe that such mutations were preferentially induced at guanine sequences (S3 Table). Although the mechanistic insight of these mutations is unclear, an unrepaired ribonucleotide in the template DNA might influence the fidelity of local DNA replication through unknown mechanism (Fig 4).

Overall, we provided evidence that mutagenesis induced by an embedded ribonucleotide is mediated by TDP1 and TDP2 in human cells. In mammals, loss of the canonical RER is associated with the risk of prostate adenocarcinoma, glioma, and colorectal tumor due to the large scale genomic instability [7,49,55]. The TDPs-dependent repair cascade may therefore contribute to the tumorigenesis in RER-deficient cells. If so, such alternate repair pathway can be a candidate target for cancer therapy. Further investigations are needed to thoroughly understand the impact of the processing pathway for the accumulated ribonucleotides in the mammalian genomic DNA.

Acknowledgements

We thank Dr. Tomonari Matsuda (Kyoto University) for providing us with pMY189 and Dr. Tatsuo Nunoshiba (International Christian University) for providing us with KS40/pOF105. We are grateful to Dr. Masami Yamada (National Defense Academy) for the helpful suggestions. We wish to thank Ms. Akiko Ukai (National Institute of Health Sciences) for technical support.

References

ASassa, MYasui, MHonma (2019) Current perspectives on mechanisms of ribonucleotide incorporation and processing in mammalian DNA. Genes Environ 41: 3 10.1186/s41021-019-0118-7

MAReijns, BRabe, RERigby, PMill, KRAstell, et al (2012) Enzymatic removal of ribonucleotides from DNA is essential for mammalian genome integrity and development. Cell 149: 10081022. 10.1016/j.cell.2012.04.011

MMGhodgaonkar, FLazzaro, MOlivera-Pimentel, MArtola-Boran, PCejka, et al (2013) Ribonucleotides misincorporated into DNA act as strand-discrimination signals in eukaryotic mismatch repair. Mol Cell 50: 323332. 10.1016/j.molcel.2013.03.019

SALujan, JSWilliams, ARClausen, ABClark, TAKunkel (2013) Ribonucleotides are signals for mismatch repair of leading-strand replication errors. Mol Cell 50: 437443. 10.1016/j.molcel.2013.03.017

JMPryor, MPConlin, JCarvajal-Garcia, MELuedeman, AJLuthman, et al (2018) Ribonucleotide incorporation enables repair of chromosome breaks by nonhomologous end joining. Science 361: 11261129. 10.1126/science.aat2477

BHiller, AHoppe, CHaase, CHiller, NSchubert, et al (2018) Ribonucleotide Excision Repair Is Essential to Prevent Squamous Cell Carcinoma of the Skin. Cancer Res 78: 59175926. 10.1158/0008-5472.CAN-18-1099

KAden, KBartsch, JDahl, MAMReijns, DEsser, et al (2019) Epithelial RNase H2 Maintains Genome Integrity and Prevents Intestinal Tumorigenesis in Mice. Gastroenterology 156: 145159 e119. 10.1053/j.gastro.2018.09.047

YJCrow, ALeitch, BEHayward, AGarner, RParmar, et al (2006) Mutations in genes encoding ribonuclease H2 subunits cause Aicardi-Goutieres syndrome and mimic congenital viral brain infection. Nat Genet 38: 910916. 10.1038/ng1842

PSEder, JAWalder (1991) Ribonuclease H from K562 human erythroleukemia cells. Purification, characterization, and substrate specificity. J Biol Chem 266: 64726479.

10 

HSJeong, PSBacklund, HCChen, AAKaravanov, RJCrouch (2004) RNase H2 of Saccharomyces cerevisiae is a complex of three proteins. Nucleic Acids Res 32: 407414. 10.1093/nar/gkh209

11 

JLSparks, HChon, SMCerritelli, TAKunkel, EJohansson, et al (2012) RNase H2-initiated ribonucleotide excision repair. Mol Cell 47: 980986. 10.1016/j.molcel.2012.06.035

12 

BHiller, MAchleitner, SGlage, RNaumann, RBehrendt, et al (2012) Mammalian RNase H2 removes ribonucleotides from DNA to maintain genome integrity. J Exp Med 209: 14191426. 10.1084/jem.20120876

13 

JSekiguchi, SShuman (1997) Site-specific ribonuclease activity of eukaryotic DNA topoisomerase I. Mol Cell 1: 8997. 10.1016/s1097-2765(00)80010-6

14 

NKim, SNHuang, JSWilliams, YCLi, ABClark, et al (2011) Mutagenic processing of ribonucleotides in DNA by yeast topoisomerase I. Science 332: 15611564. 10.1126/science.1205016

15 

NAdachi, SSo, HKoyama (2004) Loss of nonhomologous end joining confers camptothecin resistance in DT40 cells. Implications for the repair of topoisomerase I-mediated DNA damage. J Biol Chem 279: 3734337348. 10.1074/jbc.M313910200

16 

MZimmermann, OMurina, MAMReijns, AAgathanggelou, RChallis, et al (2018) CRISPR screens identify genomic ribonucleotides as a source of PARP-trapping lesions. Nature 559: 285289. 10.1038/s41586-018-0291-z

17 

SNHuang, JSWilliams, MEArana, TAKunkel, YPommier (2017) Topoisomerase I-mediated cleavage at unrepaired ribonucleotides generates DNA double-strand breaks. EMBO J 36: 361373. 10.15252/embj.201592426

18 

SWYang, ABBurginJr., BNHuizenga, CARobertson, KCYao, et al (1996) A eukaryotic enzyme that can disjoin dead-end covalent complexes between DNA and type I topoisomerases. Proc Natl Acad Sci U S A 93: 1153411539. 10.1073/pnas.93.21.11534

19 

HInterthal, HJChen, JJChampoux (2005) Human Tdp1 cleaves a broad spectrum of substrates, including phosphoamide linkages. J Biol Chem 280: 3651836528. 10.1074/jbc.M508898200

20 

SYHuang, JMurai, IDalla Rosa, TSDexheimer, ANaumova, et al (2013) TDP1 repairs nuclear and mitochondrial DNA damage induced by chain-terminating anticancer and antiviral nucleoside analogs. Nucleic Acids Res 41: 77937803. 10.1093/nar/gkt483

21 

JMurai, SYHuang, BBDas, TSDexheimer, STakeda, et al (2012) Tyrosyl-DNA phosphodiesterase 1 (TDP1) repairs DNA damage induced by topoisomerases I and II and base alkylation in vertebrate cells. J Biol Chem 287: 1284812857. 10.1074/jbc.M111.333963

22 

YPommier, SYHuang, RGao, BBDas, JMurai, et al (2014) Tyrosyl-DNA-phosphodiesterases (TDP1 and TDP2). DNA Repair (Amst) 19: 114129. 10.1016/j.dnarep.2014.03.020

23 

FCortes Ledesma, SFEl Khamisy, MCZuma, KOsborn, KWCaldecott (2009) A human 5'-tyrosyl DNA phosphodiesterase that repairs topoisomerase-mediated DNA damage. Nature 461: 674678. 10.1038/nature08444

24 

ZZeng, ASharma, LJu, JMurai, LUmans, et al (2012) TDP2 promotes repair of topoisomerase I-mediated DNA damage in the absence of TDP1. Nucleic Acids Res 40: 83718380. 10.1093/nar/gks622

25 

MTsuda, KKitamasu, CKumagai, KSugiyama, TNakano, et al (2020) Tyrosyl-DNA phosphodiesterase 2 (TDP2) repairs topoisomerase 1 DNA-protein crosslinks and 3'-blocking lesions in the absence of tyrosyl-DNA phosphodiesterase 1 (TDP1). DNA Repair (Amst) 91–92: 102849 10.1016/j.dnarep.2020.102849

26 

JHeo, JLi, MSummerlin, AHays, SKatyal, et al (2015) TDP1 promotes assembly of non-homologous end joining protein complexes on DNA. DNA Repair (Amst) 30: 2837. 10.1016/j.dnarep.2015.03.003

27 

FGomez-Herreros, RRomero-Granados, ZZeng, AAlvarez-Quilon, CQuintero, et al (2013) TDP2-dependent non-homologous end-joining protects against topoisomerase II-induced DNA breaks and genome instability in cells and in vivo. PLoS Genet 9: e1003226 10.1371/journal.pgen.1003226

28 

JLSparks, PMBurgers (2015) Error-free and mutagenic processing of topoisomerase 1-provoked damage at genomic ribonucleotides. EMBO J 34: 12591269. 10.15252/embj.201490868

29 

AAlvarez-Quilon, JLWojtaszek, MCMathieu, TPatel, CDAppel, et al (2020) Endogenous DNA 3' Blocks Are Vulnerabilities for BRCA1 and BRCA2 Deficiency and Are Reversed by the APE2 Nuclease. Mol Cell 78: 11521165 e1158. 10.1016/j.molcel.2020.05.021

30 

ASassa, HTada, ATakeishi, KHarada, MSuzuki, et al (2019) Processing of a single ribonucleotide embedded into DNA by human nucleotide excision repair and DNA polymerase eta. Sci Rep 9: 13910 10.1038/s41598-019-50421-8

31 

EFanning, KZhao (2009) SV40 DNA replication: from the A gene to a nanomachine. Virology 384: 352359. 10.1016/j.virol.2008.11.038

32 

LCTsui, MLBreitman, LSiminovitch, MBuchwald (1982) Persistence of freely replicating SV40 recombinant molecules carrying a selectable marker in permissive simian cells. Cell 30: 499508. 10.1016/0092-8674(82)90247-1

33 

JBNeto, AGentil, RECabral, ASarasin (1992) Mutation spectrum of heat-induced abasic sites on a single-stranded shuttle vector replicated in mammalian cells. J Biol Chem 267: 1971819723.

34 

SSHeinzel, PJKrysan, MPCalos, RBDuBridge (1988) Use of simian virus 40 replication to amplify Epstein-Barr virus shuttle vectors in human cells. J Virol 62: 37383746. 10.1128/JVI.62.10.3738-3746.1988

35 

MHonma, MIzumi, MSakuraba, STadokoro, HSakamoto, et al (2003) Deletion, rearrangement, and gene conversion; genetic consequences of chromosomal double-strand breaks in human cells. Environ Mol Mutagen 42: 288298. 10.1002/em.10201

36 

NNHoa, TShimizu, ZWZhou, ZQWang, RADeshpande, et al (2016) Mre11 Is Essential for the Removal of Lethal Topoisomerase 2 Covalent Cleavage Complexes. Mol Cell 64: 580592. 10.1016/j.molcel.2016.10.011

37 

RAkagawa, HTTrinh, LKSaha, MTsuda, KHirota, et al (2020) UBC13-Mediated Ubiquitin Signaling Promotes Removal of Blocking Adducts from DNA Double-Strand Breaks. iScience 23: 101027 10.1016/j.isci.2020.101027

38 

TSuzuki, HHarashima, HKamiya (2010) Effects of base excision repair proteins on mutagenesis by 8-oxo-7,8-dihydroguanine (8-hydroxyguanine) paired with cytosine and adenine. DNA Repair (Amst) 9: 542550. 10.1016/j.dnarep.2010.02.004

39 

FObata, TNunoshiba, THashimoto-Gotoh, KYamamoto (1998) An improved system for selection of forward mutations in an Escherichia coli supF gene carried by plasmids. J Radiat Res 39: 263270. 10.1269/jrr.39.263

40 

KSatou, MHori, KKawai, HKasai, HHarashima, et al (2009) Involvement of specialized DNA polymerases in mutagenesis by 8-hydroxy-dGTP in human cells. DNA Repair (Amst) 8: 637642. 10.1016/j.dnarep.2008.12.009

41 

VKellner, BLuke (2020) Molecular and physiological consequences of faulty eukaryotic ribonucleotide excision repair. EMBO J 39: e102309 10.15252/embj.2019102309

42 

SMCerritelli, RJCrouch (2019) RNase H2-RED carpets the path to eukaryotic RNase H2 functions. DNA Repair (Amst) 84: 102736 10.1016/j.dnarep.2019.102736

43 

MHonma, MHayashi, TSofuni (1997) Cytotoxic and mutagenic responses to X-rays and chemical mutagens in normal and p53-mutated human lymphoblastoid cells. Mutat Res 374: 8998. 10.1016/s0027-5107(96)00223-0

44 

JWhitwell, RSmith, KJenner, HLyon, DWood, et al (2015) Relationships between p53 status, apoptosis and induction of micronuclei in different human and mouse cell lines in vitro: Implications for improving existing assays. Mutat Res Genet Toxicol Environ Mutagen 789–790: 727. 10.1016/j.mrgentox.2015.05.011

45 

SYHuang, SGhosh, YPommier (2015) Topoisomerase I alone is sufficient to produce short DNA deletions and can also reverse nicks at ribonucleotide sites. J Biol Chem 290: 1406814076. 10.1074/jbc.M115.653345

46 

MAHossain, YLin, SYan (2018) Single-Strand Break End Resection in Genome Integrity: Mechanism and Regulation by APE2. Int J Mol Sci 19 10.3390/ijms19082389

47 

XHan, FMayca Pozo, JNWisotsky, BWang, JWJacobberger, et al (2015) Phosphorylation of Minichromosome Maintenance 3 (MCM3) by Checkpoint Kinase 1 (Chk1) Negatively Regulates DNA Replication and Checkpoint Activation. J Biol Chem 290: 1237012378. 10.1074/jbc.M114.621532

48 

EPetermann, MWoodcock, THelleday (2010) Chk1 promotes replication fork progression by controlling replication initiation. Proc Natl Acad Sci U S A 107: 1609016095. 10.1073/pnas.1005031107

49 

CWang, GWang, XFeng, PShepherd, JZhang, et al (2019) Genome-wide CRISPR screens reveal synthetic lethality of RNASEH2 deficiency and ATR inhibition. Oncogene 38: 24512463. 10.1038/s41388-018-0606-4

50 

FLi, QWang, JHSeol, JChe, XLu, et al (2019) Apn2 resolves blocked 3' ends and suppresses Top1-induced mutagenesis at genomic rNMP sites. Nat Struct Mol Biol 26: 155163. 10.1038/s41594-019-0186-1

51 

YWZhang, MRegairaz, JASeiler, KKAgama, JHDoroshow, et al (2011) Poly(ADP-ribose) polymerase and XPF-ERCC1 participate in distinct pathways for the repair of topoisomerase I-induced DNA damage in mammalian cells. Nucleic Acids Res 39: 36073620. 10.1093/nar/gkq1304

52 

BBDas, SAntony, SGupta, TSDexheimer, CERedon, et al (2009) Optimal function of the DNA repair enzyme TDP1 requires its phosphorylation by ATM and/or DNA-PK. EMBO J 28: 36673680. 10.1038/emboj.2009.302

53 

JLi, MSummerlin, KCNitiss, JLNitiss, LAHanakahi (2017) TDP1 is required for efficient non-homologous end joining in human cells. DNA Repair (Amst) 60: 4049. 10.1016/j.dnarep.2017.10.003

54 

JECho, SJinks-Robertson (2019) Deletions associated with stabilization of the Top1 cleavage complex in yeast are products of the nonhomologous end-joining pathway. Proc Natl Acad Sci U S A 116: 2268322691. 10.1073/pnas.1914081116

55 

UBeyer, FBrand, HMartens, JWeder, AChristians, et al (2017) Rare ADAR and RNASEH2B variants and a type I interferon signature in glioma and prostate carcinoma risk and tumorigenesis. Acta Neuropathol 134: 905922. 10.1007/s00401-017-1774-y