Transcription-associated topoisomerase 2α (TOP2A) activity is a major effector of cytotoxicity induced by G-quadruplex ligands
Abstract
G-quadruplexes (G4) are non-canonical DNA structures found in the genome of most species including human. Small molecules stabilizing these structures, called G4 ligands, have been identified and, for some of them, shown to induce cytotoxic DNA double-strand breaks. Through the use of an unbiased genetic approach, we identify here topoisomerase 2α (TOP2A) as a major effector of cytotoxicity induced by two clastogenic G4 ligands, pyridostatin and CX-5461, the latter molecule currently undergoing phase I/II clinical trials in oncology. We show that both TOP2 activity and transcription account for DNA break production following G4 ligand treatments. In contrast, clastogenic activity of these G4 ligands is countered by topoisomerase 1 (TOP1), which limits co-transcriptional G4 formation, and by factors promoting transcriptional elongation. Altogether our results support that clastogenic G4 ligands act as DNA structure-driven TOP2 poisons at transcribed regions bearing G4 structures.
Introduction
In recent years, evidence has accumulated to indicate that transcription is a major source of genomic instability (Aguilera, 2002; Gaillard et al., 2013; Gómez-González and Aguilera, 2019). Transcription-dependent DNA double-stranded breaks (DSBs) are mainly associated with RNA-polymerase (RNA-Pol) arrests provoked by different non-exclusive factors including DNA torsional stress, inhibition of transcription elongation, and formation of secondary structures, such as G-quadruplexes (G4) and R-loops (Chedin and Benham, 2020; Crossley et al., 2019; García-Muse and Aguilera, 2019; Kotsantis et al., 2020; Miglietta et al., 2020). G4 are four-stranded secondary structures formed at guanine-rich tracts (Burge et al., 2006). Present throughout the human genome (Chambers et al., 2015; Hänsel-Hertsch et al., 2016; Hänsel-Hertsch et al., 2018), G4 have been associated with spontaneous DNA breaks, hotspots for chromosomal translocations and several human syndromes (Hänsel-Hertsch et al., 2017; Murat and Balasubramanian, 2014; Maizels, 2015; Puget et al., 2019). In proliferating cells, G4 act as replication fork barriers, provoking fork collapses, the activation of the DNA damage response and the induction of replication-dependent DSBs (Lerner and Sale, 2019). In addition, increasing evidence also indicates a significant impact of G4 structures on genomic stability through transcription-dependent processes (Kotsantis et al., 2020; Miglietta et al., 2020; De Magis et al., 2019; Rodriguez et al., 2012; van Wietmarschen et al., 2018; Kim, 2019).
G4 mapping in the human genome shows a significant enrichment of these structures within promoter and 5′ UTR regions of highly transcribed genes, and several in vitro and cellular studies show that the stabilization of G4 structures by small compounds, G4 ligands, generally represses transcription of genes containing G-rich tracts (Chambers et al., 2015; Hänsel-Hertsch et al., 2016; Hänsel-Hertsch et al., 2017). During transcription, while G4 structures located on template DNA could act as physical barriers blocking RNA-Pol II progression, the formation of G4 on the opposite strand could promote and stabilize secondary structures that block transcription elongation (Gómez-González and Aguilera, 2019; Kotsantis et al., 2020; Miglietta et al., 2020). Genome-wide analyses of G4 motifs in human cells indicate that these structures are highly correlated with RNA-Pol II pausing sites and R-loop-forming regions, two different factors promoting RNA-Pol II arrests and transcription-dependent DNA breaks (Puget et al., 2019; Chen et al., 2017; Eddy et al., 2011). Additionally, the formation of highly stable DNA secondary structures, such as G4, has been shown to promote the formation of DNA topoisomerase 2 (TOP2)-mediated DNA breaks (Szlachta et al., 2020). Interestingly, recent studies demonstrate a major contribution of TOP2 activity in the generation of DSBs in highly transcribed genes (Canela et al., 2019; Gittens et al., 2019; Gothe et al., 2019). Moreover, several studies demonstrated that DNA topoisomerase 1 (TOP1) and TOP2 recognize and preferentially cleave DNA at regions forming stable secondary structures (Froelich-Ammon et al., 1994; Jonstrup et al., 2008; Mills et al., 2018; West and Austin, 1999). In eukaryotic cells, TOP1 and TOP2 activities are required to resolve topological stresses resulting from DNA transactions (for a review, see Pommier et al., 2016). These enzymes relax topological constraints through the formation of transient single-stranded (TOP1) or double-stranded DNA breaks (TOP2), in which the enzymes are covalently linked to the DNA backbone. In humans, TOP2 activity is supported by two isoenzymes, TOP2α (TOP2A) and TOP2β (TOP2B), that are encoded by two different genes. TOP2A plays key roles in DNA replication and chromosome segregation, while TOP2B is mainly associated with transcription (Pommier et al., 2016; Nitiss, 2009a; Madabhushi, 2018). TOP2 are poisoned by small molecules that trap the transient TOP2-DNA complex, also known as ‘cleavage complex’ (TOP2cc), during the enzyme catalytic cycle (Deweese and Osheroff, 2009; Nitiss, 2009b; Pommier et al., 2010). The repair of TOP2cc requires a sequential process consisting in the removal of TOP2 protein from DNA through TOP2 proteolysis (Fan et al., 2008; Mao et al., 2001; Zhang et al., 2006) or nucleolytic degradation of the associated DNA (Aparicio et al., 2016; Nakamura et al., 2010) and the repair of the resulting DSB by non-homologous end joining (NHEJ) (Gómez-Herreros et al., 2013) or homologous recombination (HR), respectively (Ashour et al., 2015).
Bruno et al., 2020 and Olivieri et al., 2020 have recently reported that the induction of DNA breaks by two clastogenic G4 ligands, CX-5461 and pyridostatin (PDS), may result from a TOP2-poisoning-like mechanism, forming TOP2cc, but the respective contribution of TOP2A/TOP2B to the clastogenic activity of both G4 ligands is still unclear. Moreover, a recent report argues that cytotoxicity of CX-5461 may rely on irreversible inhibition of RNA-Pol I transcription initiation (Mars et al., 2020). Here, we report a differential contribution of TOP2A and TOP2B to DNA break production in response to G4 ligands PDS and CX-5461 and identify TOP2A as a major effector of cytotoxicity induced by the two G4 ligands. We also report that transcription plays an essential role in the production of TOP2-dependent DNA breaks induced by both G4 ligands. We show that G4 ligand-induced DSBs are countered by TOP1 that limits co-transcriptional G4 formation and by factors that promote transcriptional elongation. Taken together, our results support the concept that CX-5461 and PDS act as DNA structure-driven TOP2 poisons at G-rich transcribed genomic regions.
Results
TOP2A mutations confer resistance to clastogenic G4 ligands CX-5461 and pyridostatin
The chemical compound CX-5461, currently in phase I/II clinical trials for cancer treatments (Khot et al., 2019), was first described as an RNA-Pol I inhibitor (Haddach et al., 2012). Although the cytotoxic effects induced by this compound have been related to rDNA-transcription inhibition (Mars et al., 2020; Negi and Brown, 2015), CX-5461 has also been shown to be a potent G4 stabilizer and to provoke rapid induction of DSBs through a TOP2-dependent mechanism (Bruno et al., 2020; Olivieri et al., 2020; Xu et al., 2017). In order to more clearly define how CX-5461 mediates its cytotoxicity, we adopted an unbiased approach based on the selection and characterization of cells resistant to this drug. To do this, human near haploid HAP1 cells were randomly mutagenized with ethyl methane sulfonate (EMS) based on previous work (Forment et al., 2017) and clones resistant to a lethal CX-5461 concentration of 0.3 µM were isolated (CX-5461-resistant [CXR] clones). Resistance of seven of these clones to CX-5461 was further confirmed by cell survival assays showing IC50 values on CXR clones ranging from 0.22 to 0.34 µM CX-5461, corresponding to an average ninefold increase in the IC50 value compared to the 0.03 µM IC50 on wild-type HAP1 cells (WT) (see Figure 1A and Supplementary file 1). In addition, CXR clones did not show cross-resistance to the unrelated drug and efflux-pumps substrate nocodazole (Kasap et al., 2014), excluding a multidrug resistance (MDR) phenotype (Supplementary file 1).
Inspired by previous work (Wacker et al., 2012), we analyzed the selected CXR clones through a global RNA-sequencing approach (RNA-seq) to identify non- and mis-sense mutations in coding genes that could account for the observed resistance. Around eight genes per clone were found with non- or mis-sense mutations through this approach (Figure 1B and Supplementary file 2). Unbiased analysis of genes mutated in several clones revealed that each clone, except CXR #A2 and CXR #A6, carried a homozygous mutation in the TOP2A gene, encoding for the TOP2A protein (Figure 1B). Manual analysis of the sequencing data for the TOP2A gene confirmed these mutations and revealed that the CXR #A2 clone carried the S654N mutation, while the CXR #A6 clone carried a homozygous mutation of the first nucleotide of the last intron, resulting in intron retention and replacement of the last 42 TOP2A amino acids, carrying the nuclear localization signal (NLS), by 18 unrelated amino acids (Figure 1—figure supplement 1A–C). From these analyses, TOP2A emerged as the only gene with coding mutations in all resistant clones. Four clones had a mutation in the ATPase domain (F85I), while two had mutations in the DNA-binding region (S654N and L703I) (Figure 1C). Immunoblotting analysis revealed that none of the identified TOP2A mutations resulted in the loss of TOP2A protein, in agreement with its essential function in proliferating cells (Akimitsu et al., 2003; Carpenter and Porter, 2004). In addition, TOP2B expression level was unaffected in these clones (Figure 1D).
To test whether the catalytic activity of the TOP2A protein was altered in CXR clones, we determined the sensitivity of WT and CXR HAP1 clones to the TOP2 poison etoposide (ETP), a chemotherapeutic drug that acts by stabilizing TOP2cc and the cytotoxicity of which is therefore dependent on TOP2 activity (Walker and Nitiss, 2002). All CXR presented a strong resistance to ETP with resistance indexes ranging from 8- to 19-fold relative to control cells (Figure 2C and Supplementary file 1). These results support the idea that TOP2A point mutations in CXR clones reduce TOP2A activity. To confirm this, we adapted a heparin-based extraction protocol (de Campos-Nebel et al., 2016) coupled with immunoblotting to monitor the accumulation of TOP2Acc following ETP treatment. In this assay, TOP2 molecules not covalently attached to DNA are extracted by heparin to a soluble fraction, while TOP2ccs are resistant to heparin and can be analyzed by immunoblotting of the insoluble pellet fraction. As shown in Figure 2—figure supplement 1, the point mutations present in clones CXR #A6 and CXR #A1 decrease the amount of TOP2Acc following ETP treatment.
In parallel, using the same genetic approach, we isolated F14R HAP1 clones resistant to a lethal concentration of 30 nM F14512, a potent and selective TOP2A poison (Bombarde et al., 2017). Targeted sequencing of TOP2A cDNA in five F14R clones revealed five different TOP2A mutations, confirming that TOP2A is the main mediator of F14512 cytotoxicity (Bombarde et al., 2017). One mutation was found in the transducer domain (W414L), while the four others lay in the DNA-binding domain (Y590H, P593S, G777D, and P890L; Figure 1C). All of these mutations were different from the ones found in the CXR clones, which suggests that, despite both acting through TOP2A, CX-54161 and F14512 affect TOP2A activity differently Strikingly, cell survival assays clearly demonstrated that all of the F14R cells were cross-resistant to CX-5461 (Figure 2D), while showing no resistance to nocodazole (Supplementary file 1), with IC50 values for CX-5461 3.6- to 16.8-fold higher than the IC50 value of CX-5461 in WT cells (Figure 2D and Supplementary file 1), supporting the concept that the TOP2A mutations found in F14R clones also conferred resistance to the G4 ligand CX-5461.
Importantly, population doubling time determination of five clones carrying different TOP2A mutations and flow cytometry analyses of 5-ethynyl-2′-deoxyuridine (EdU) incorporation in CXR #A6 and F14R #C1 cells revealed no significant differences in the proliferation rate of TOP2A mutated clones as compared to WT cells (Figure 2—figure supplement 2A, B).
In addition, the survival of WT and TOP2A mutated clones to camptothecin, a DNA TOP1 poison, or calicheamicin, a radio-mimetic compound, was not significantly different. Since camptothecin induces toxic replication-associated DSBs, which are mainly repaired by HR, while calicheamicin induces DSB in all cell cycle phases, which are repaired by both NHEJ and HR (Elmroth et al., 2003; Furuta et al., 2003), these results indicate that DSB repair mechanisms are proficient in TOP2A mutated cells (Figure 2—figure supplement 2C). Supporting this conclusion, the sensitization to calicheamicin induced by an inhibitor of the NHEJ factor DNA-PK was similar in WT and TOP2A mutated cells (Figure 2—figure supplement 2D) and the kinetics of DSB repair after X-ray irradiation measured by flow cytometry analysis of the γH2AX signal were identical between WT and mutant cells (Figure 2—figure supplement 2E). Finally, western blot analysis showed that TOP1 levels were unaffected in TOP2A mutated cells (Figure 2—figure supplement 2F). Altogether, these results indicate that the resistance to CX-5461 observed in CXR and F14R cells is not associated with major changes in proliferation or DNA repair mechanisms.
To investigate whether CX-5461 resistance induced by TOP2A mutations extends to other G4 stabilizers, we tested CXR mutants for their cross-resistance to pyridostatin (PDS), one of the best characterized G4 ligands so far (Rodriguez et al., 2008). In cells, PDS treatment, similar to CX-5461, induces a rapid accumulation of DSBs, but in contrast to CX-5461, PDS does not affect RNA-Pol I activity (Rodriguez et al., 2012). Cell survival assays showed that all CXR clones were cross-resistant to PDS (IC50 values for PDS 2.7- to 4.4-fold higher than the IC50 value of PDS on WT cells) but with lower resistance indexes than those observed for CX-5461 and F14512 (Figure 2E and Supplementary file 1). More remarkably, cross-resistance studies established that all F14R clones were also resistant to PDS (Figure 2F and Supplementary file 1). The role of TOP2A in the cytotoxicity of CX-5461 and PDS was further supported by the marked resistance to both CX-5461, PDS, and ETP (as a control) conferred by small-interfering RNA-mediated depletion of TOP2A in HAP1 cells, with IC50 values similar to those observed for mutant cells (Figure 2G).
In cells, stabilization of G4 by ligands has been shown to provoke cell growth modifications relying on DNA and RNA transaction alterations (Varshney et al., 2020). The acridine derivative RHPS4, a potent G4 ligand, affects both telomere maintenance and mitochondrial DNA replication (Falabella et al., 2019; Salvati et al., 2007), while two selective bisquinolinium G4 ligands, 360A and PhenDC3, impact on telomeres, gene expression, and RNA metabolism (Dumas et al., 2021; Halder et al., 2012; Pennarun et al., 2005). Through performing cell viability with RPHPS4 and PhenDC3, we revealed that the TOP2A mutant cells do not have significantly different modifications in response to these two chemically unrelated ligands as compared to WT cells (Figure 2—figure supplement 2G).
Altogether, these results indicate that TOP2A is specifically involved in the cytotoxic effect induced by the two G4 stabilizers CX-5461 and PDS.
Rapid accumulation of DNA DSBs upon G4 ligand treatment depends on TOP2 activity
In human cells, short treatments with PDS or CX-5461 induce rapid production of the DSB markers γH2AX and 53BP1 foci (Rodriguez et al., 2012; Xu et al., 2017; Figure 3A, Figure 3—figure supplement 1A, and Figure 3—figure supplement 2A). G4-dependent γH2AX foci production is increased in cells incubated with the DNA-PK inhibitor NU7441 (DNA-PKi; Figure 3—figure supplement 2B), indicating that a substantial number of G4-induced DNA breaks are repaired through the DNA-PK-dependent NHEJ pathway, the major DSB repair pathway in human cells (Pannunzio et al., 2018). Considering the main contribution of DSBs to the cytotoxic effect of several anticancer agents and having shown that TOP2A activity determines the cytotoxic effect of G4 ligands (Figures 1 and 2), we evaluated the role of TOP2A in DSB production upon PDS and CX-5461 treatments. First, in CXR cells carrying different mutations in the TOP2A protein, clones CXR #A1 (F85I), #A2 (S654N), and #A6 (intron retention), γH2AX production was significantly reduced as compared to control cells (Figure 3A). We assessed the role of TOP2 activities on DSBs production induced by PDS and CX-5461 treatments in HeLa cells by studying the impact of the TOP2 catalytic inhibitor BNS-22 (Kawatani et al., 2011; Figure 3B, Figure 3—figure supplement 1B, and Figure 3—figure supplement 2C). Pre-incubation with BNS-22 significantly decreased the number of γH2AX foci in cells treated with both G4 ligands, thereby demonstrating the major role of TOP2 catalytic activity in the production of DNA breaks following G4 ligand treatments. To evaluate the respective contribution of TOP2A and TOPB2B proteins in the formation of G4-dependent DSBs, we analyzed γH2AX production in HeLa cells transfected by siRNA against each TOP2 as compared to cells transfected with control siRNA. Immunofluorescence studies showed that in TOP2A-depleted cells γH2AX production was abolished upon PDS treatment and strongly reduced upon CX-5461 treatment (Figure 3C, Figure 3—figure supplement 1C, and Figure 3—figure supplement 2D). In contrast, TOP2B depletion did not impact γH2AX production upon PDS treatment, while it reduced γH2AX production upon CX-5461 treatment to a similar extent to that induced by TOP2A depletion. Moreover, simultaneous siRNA-mediated knock-down of TOP2A and TOP2B further decreased γH2AX production by CX-5461 when compared to the effect of separate siRNAs (Figure 3C and Figure 3—figure supplement 1C). Altogether, these results indicate that TOP2A activity is the main effector of DSBs by clastogenic G4 ligands in human cells, with TOP2B variably contributing depending on the ligand class.
TOP2-dependent DSBs induced by CX-5461 and PDS are associated with G4 structures in cells
To investigate whether TOP2-dependent DSBs induced by CX-5461 and PDS are associated with G4 structures, we adapted the heparin pre-extraction protocol to immunofluorescence in order to visualize the formation of covalent TOP2A-DNA complex (TOP2Acc) in cells. To validate this novel assay, TOP2A and γH2AX signals were quantified upon ETP treatment. As shown in Figure 4A, B, ETP provoked a strong increase in TOP2A signal that correlated with a similar increase in γH2AX foci (Figure 4A, B and Figure 4—figure supplement 1A, B). Using the same approach, we observed a significant increase of TOP2Acc and γH2AX signals in cells treated with CX-5461 and PDS (Figure 4A, B and Figure 4—figure supplement 1A, B), supporting the idea that these G4 ligands mediate their cytotoxic effect through stabilizing TOP2Acc.
The bisquinolinium compound 360A is a potent and selective G4 binder that exhibits a lack of clastogenic activity in most human cells as compared to CX-5461 and PDS (Figure 4A and Figure 4—figure supplement 1A). Assuming that the observed formation of TOP2Acc upon CX-5461 and PDS treatment results from the stabilization of G4 structures, we hypothesized that the binding of 360A to the same structures would impair that of CX-5461 and PDS and thereby inhibit the formation of TOP2Acc in cells treated with the latter clastogenic ligands. As shown in Figure 4B (and Figure 4—figure supplement 1B), 360A pre-treatment strongly reduced the induction of TOP2Acc by both clastogenic G4 ligands, without impacting on the TOP2Acc signal induced by ETP. As expected, inhibition by 360A pre-treatment of TOP2Acc formation resulted in a significant reduction in the γHA2X signal induced by clastogenic G4 ligands, without impacting on the formation of DSBs in ETP-treated cells (Figure 4A and Figure 4—figure supplement 1A).
To correlate TOP2Acc with stabilization of G4 structures in cells treated by CX-5461 and PDS, we used high-resolution imaging (3D-SIM) to perform colocalization studies of TOP2Acc with G4 as visualized with the selective anti-G4 antibody, BG4 (Biffi et al., 2013). As shown in Figure 4C (and Figure 4—figure supplement 1C ), CX-5461 and PDS treatments induced a significant increase in the percentage of BG4 signals colocalizing (Figure 4D) with TOP2Acc and induced a significant increase in the percentage of TOP2Acc signals colocalizing with G4 structures. Furthermore, analysis of BG4 signals showed a strong increase of G4 structures in all the treated conditions (Figure 4—figure supplement 2). It is noteworthy that the increase in TOP2Acc-G4 colocalizations was not observed in ETP-treated cells, showing that the percentage of colocalization does not result from a global increase of TOP2Acc or BG4 signals induced by the different treatments (Figure 4—figure supplement 2). The increased BG4 signal observed in ETP-treated cells probably results from a global increase of DNA supercoiling by ETP, which would promote G4 formation as shown in several other studies (Sekibo and Fox, 2017; Sun, 2010; Sun and Hurley, 2009).
Altogether, our data support the conclusion that TOP2-dependent DSBs induced by CX-5461 and PDS are associated with G4 structures in cells.
TOP2-dependent DSBs induced by G4 ligands are transcription-dependent
We observed that PDS-induced DSB markers appeared in all cell cycle phases, indicating that PDS-induced DNA damage is not strictly S-phase dependent (Figure 5—figure supplement 1A; Rodriguez et al., 2012). Using the DNA-base analog EdU to visualize DNA synthesis, we confirmed that DNA damage induction by PDS and CX-5461 also occurs outside of S phase (Figure 5A and Figure 5—figure supplement 2A; Rodriguez et al., 2012). Furthermore, comparison of the number of G4 ligand-induced γH2AX foci in EdU-negative cells and in the total cell population indicated that DNA replication-independent processes were very efficient for the production of DSBs by PDS (Figure 5A and Figure 5—figure supplement 2A). Accordingly, DSB production by PDS in HeLa cells was strongly reduced upon inhibition of RNA-Pol II-dependent transcription by 5,6-dichloro-1-b-D-ribofuranosylbenzimidazole (DRB), an inhibitor of critical phosphorylations in the RNA-Pol-II C-terminal domain (Yankulov et al., 1995; Figure 5B and Figure 5—figure supplement 2B). Importantly, immunofluorescence analysis with the G4-specific antibody BG4 showed that inhibition of RNA-Pol II-dependent transcription did not impede the accumulation of G4 structures caused by PDS treatment (Figure 5—figure supplement 1B).
DRB also reduced DSB production upon CX-5461 treatment, albeit to a lesser extent than with PDS (Figure 5B and Figure 5—figure supplement 2B). Thus, we next investigated the contribution of RNA-Pol I-dependent transcription to γH2AX production upon CX-5461 treatment. BHM21, a DNA intercalator showing a preferential binding to GC-rich sequences (Peltonen et al., 2014) but not to G4 structures (Xu et al., 2017), is a potent RNA-Pol I inhibitor that causes the proteasome-dependent degradation of RPA194, the large catalytic subunit protein of RNA-Pol I holoenzyme (Peltonen et al., 2014). In HeLa cells, pre-treatment with BMH21 significantly reduced the formation of γH2AX signals induced by CX-5461 treatment (Figure 5C and Figure 5—figure supplement 2C). Remarkably, the concomitant inhibition of RNA-Pol II and RNA-Pol I-dependent transcription strongly reduced DSB production by this compound (Figure 5D and Figure 5—figure supplement 2D).
Altogether, these data support a key role of transcription in the production of DSBs following the stabilization of G4 structures by PDS and CX-5461.
In human cells, CX-5461 and PDS have been shown to provoke a rapid accumulation of R-loops (De Magis et al., 2019; Sanij et al., 2020; Shen et al., 2017; Duquette et al., 2004), which have been associated with the production of transcription-dependent DSBs and genomic instability (De Magis et al., 2019). Thus, in order to investigate the impact of R-loops in the DNA damage production by PDS and CX-5461, we used a U2OS cell line expressing the Escherichia coli RNaseHI, which cleaves and removes R-loops, under the control of a doxycycline-inducible promoter (Britton et al., 2014). Upon RNaseHI expression, we observed a significant decrease of DNA damage signals induced by PDS and CX-5461 (Figure 5E and Figure 5—figure supplement 2E), indicating that the production of R-loop structures contributes to the formation of DNA breaks following G4 stabilization. Altogether, these results indicate that transcription elongation plays a key role in the production of DNA damage by both G4 ligands.
TOP2-dependent DSBs induced by G4 stabilizers are countered by TOP1 and factors promoting transcription elongation
Since the formation of G4 structures is linked to negative supercoiling, RNA-Pol II pausing, and R-loop formation, three transcription-dependent processes involving TOP1 activity (Kim and Jinks-Robertson, 2017), we evaluated the role of TOP1 enzyme in the formation of DSBs following treatment with G4 ligands. Strikingly, shRNA-mediated depletion of TOP1 protein in human cells caused a significant increase in DNA damage induced by CX-5461 and PDS (Figure 6A and Figure 6—figure supplement 1A). Furthermore, immunofluorescence studies showed that the increase in DNA damage signals in TOP1-depleted cells was dependent on TOP2 activity since BNS22 pre-treatment in these cells blocked the induction of γH2AX signals by CX-5461 and PDS (Figure 6—figure supplement 2A). Confirming the role of TOP2, DNA break formation upon PDS treatment in TOP1 knock-down cells was also impaired by TOP2A depletion (Figure 6—figure supplement 2B).
EdU staining revealed that enhanced DNA damage production in TOP1 knock-down cells is not restricted to S phase of the cell cycle (Figure 6—figure supplement 3B). Moreover, DRB pre-treatment abrogated PDS-induced DNA damage in TOP1-depleted cells, indicating that they were fully dependent on transcription (Figure 6—figure supplement 3A). Finally, the depletion of TOP1 in HeLa cells caused a significant increase in the cytotoxic effect of both PDS and CX-5461 that was reverted by pre-treatment with DRB or DRB plus BMH21 for PDS and CX-5161, respectively (Figure 6B and Figure 6—figure supplement 1B). Altogether, these results support the hypothesis that a key factor in the mechanism of TOP2-mediated DNA damage produced by clastogenic G4 ligands is the accumulation of topological stresses provoked by transcription that are countered by TOP1. Consistent with these findings, immunofluorescence analysis showed that TOP1 knock-down provoked a significant increase of the BG4 signal in human cells (Figure 6C, Figure 6—figure supplement 1C, and Figure 6—figure supplement 3C) with no cumulative effect of PDS and TOP1 depletion on BG4 signal. We concluded that the accumulation of transcription-dependent negative supercoiling resulting from TOP1 depletion has a major impact on G4 formation. In line with this result, kinetic studies of γH2AX production in PDS-treated cells showed that TOP1 depletion significantly accelerates the formation of PDS-induced DNA damage signals compared to control cells (Figure 6D and Figure 6—figure supplement 1D), indicating that TOP1 depletion facilitates DSB production by G4 ligands.
During transcription, TOP1 activity is enhanced by the BRD4-dependent phosphorylation of RNA-Pol II CTD (Baranello et al., 2016). BRD4 also associates with the positive transcription elongation factor (PTEF) to promote transcription and release RNA-Pol II from paused sites (Jang et al., 2005; Yang et al., 2005). In cells, BRD4 depletion provokes R-loop accumulation that drives transcription-dependent DNA damage (Kim et al., 2019; Lam et al., 2020). Interestingly, we showed through immunofluorescence analysis that the inhibition of BRD4 activity by the potent inhibitor JQ1 provoked a significant increase of DNA breaks induced by PDS (Figure 6E and Figure 6—figure supplement 1E), indicating that inhibition of transcription elongation facilitates the formation of DSBs upon G4 ligand treatment. Altogether, our results strongly suggest that TOP1 antagonizes G4 ligand action through a transcription-dependent mechanism that is related to RNA polymerase stalling.
Discussion
In this study, using an unbiased genetic approach, we identified the TOP2A protein as the main effector of the cytotoxicity of two clastogenic G4 ligands: CX-5461 and PDS. Our study highlights the strength of the genetic approach we applied here, relying on chemical mutagenesis in a haploid background (Forment et al., 2017). Indeed, despite TOP2A being essential in proliferating cells (Akimitsu et al., 2003; Carpenter and Porter, 2004), we were able to readily isolate CX-5461-resistant clones carrying TOP2A mutations allowing direct identification of its crucial role in DSB induction upon G4 stabilization. This would not have been possible through loss of function screens, for example, using CRISPR/Cas9 or insertional mutagenesis. In addition, we produced novel mutations of TOP2A conferring resistance to both CX-5461 and F14512 that could be useful to obtain insights into TOP2A biology. It is noteworthy that these nine-point mutations, separating TOP2A essential function from its role in G4-induced DSB production, are broadly distributed throughout the TOP2A coding sequence and functional domains. This indicates that, as observed for TOP2 poisons in a yeast complementation system (Blower et al., 2019), resistance to CX-5461 can be obtained through different protein modifications and is not strictly dependent on the catalytic function. This is especially true for the mutation in the CXR #A6 clone, which resulted in the expression of a TOP2A lacking its terminal nuclear localization sequence and is therefore sequestered in the cytoplasm in interphase cells (Figure 1—figure supplement 1C), although it can access the DNA during mitosis. It is noteworthy that such a mutant would have been difficult to devise and express at proper levels in complementation experiments since TOP2A overexpression is toxic (McPherson and Goldenberg, 1998). Four of the eight amino acid changes in TOP2A identified in this work (P593S, S654N, L703I, and P890L) have been previously described and confer Vosaroxin resistance to complemented Saccharomyces cerevisiae (Blower et al., 2019). Vosaroxin is a quinolone derivative that acts as a DNA intercalator and a TOP2 inhibitor, further supporting the hypothesis that some of the TOP2A point mutations identified here alter TOP2A activity.
Differential contribution of TOP2 proteins to cytotoxicity induced by PDS and CX-5461
A major finding of our work is the differential contribution of TOP2 isoforms to cytotoxicity induced by both clastogenic G4 ligands. Indeed, while Bruno et al., 2020 have previously shown a major role for the TOP2A protein in the induction of DSBs upon CX-5461 treatment, the impact of the TOP2A activity on PDS-induced DSBs has not been clearly reported. Here, we show that TOP2A single-point mutations found in both CXR and F14R cells are sufficient to confer resistance to both CX-5461 and PDS. In addition, the major role of TOP2A in resistance to both G4 ligands was confirmed through an RNA silencing approach in HAP1 and HeLa cells. Although some slight differences in the response to G4 ligands were observed between HAP1 and HeLa cells, detailed analysis in HeLa cells of the relative contribution of the two TOP2 isoforms to DNA break production following G4 ligand treatments also supports a major role for TOP2A in the formation of DSBs by CX-5461 and PDS. Indeed, while TOP2A depletion strongly impacted on DSB production in response to PDS and CX-5461, TOP2B protein was only partially involved in DNA break production in response to CX-5461.
The main role of TOP2A in producing transcription-associated DNA breaks in response to PDS is at first intriguing since TOP2B is considered as the principal TOP2 responsible for the resolution of topological stress associated with transcription (Pommier et al., 2016; Nitiss, 2009a; Madabhushi, 2018; Austin et al., 2018). In contrast, TOP2A is believed to resolve mainly topological constraints associated with replication and chromosome segregation, in line with its increased expression during S and G2 phases (Pommier et al., 2016; Nitiss, 2009a; Akimitsu et al., 2003; Carpenter and Porter, 2004). However, some studies also implicate TOP2A in transcription (Mondal and Parvin, 2001; Ray et al., 2013; Yu et al., 2017) and its activity is required for maximal production of transcription-dependent DNA breaks induced by ETP, a potent poison of TOP2 (Tammaro et al., 2013). Moreover, genome-wide analysis of TOP2A cleavage sites shows a significant enrichment of TOP2A on highly transcribed loci (Zhang et al., 2006). Interestingly, elevated transcription levels have been shown to promote G4 formation that are favored by negative superhelicity caused by the progression of RNA-Pol complexes through DNA (Zheng et al., 2017).
Transcription drives TOP2-dependent G4-ligand-induced DSBs
Another important finding of our work is the major impact of transcription on the formation of TOP2-mediated DNA breaks following G4 ligand treatment. In this study, we show that DRB treatment, which specifically blocks RNA-Pol II transcriptional elongation, completely abrogates DSB formation induced by PDS and significantly reduces DNA break production induced by CX-5461. Interestingly, we found that DNA breaks induced by CX-5461 are strongly reduced when both RNA-Pol II and RNA-Pol I activities are inhibited, indicating an important contribution of rDNA transcription in the cellular response to this ligand. Consistent with our results, active transcription of rDNA repeats has been found to enhance the sensitivity of cells to CX-5461 and DNA damage production (Son et al., 2020). In cells, topological stresses induced by transcription are mainly resolved by the TOP1 protein (Kim and Jinks-Robertson, 2017). Here, we show evidence for a major role of the TOP1 protein in countering DNA break formation by CX-5461 and PDS. Consistent with our data, TOP1 depletion in yeast drives genomic instability at highly transcribed G4-forming sequences (Yadav et al., 2014; Yadav et al., 2016).
We also demonstrate that impairing RNA-Pol II progression by inhibiting the transcriptional elongation-promoting BRD4 protein increases the number of DSBs induced by PDS. Interestingly, the inhibition of BRD4 and TOP1 activities as well as the reduction of RNA-Pol elongation has been shown to promote R-loops (Shivji et al., 2018; Tuduri et al., 2009; Zhang et al., 2017; Kim et al., 2019; El Hage et al., 2010; Li et al., 2015), a DNA-RNA secondary structure associated with the formation of transcription-dependent DNA breaks (Aguilera, 2002; Gaillard et al., 2013; Gómez-González and Aguilera, 2019; Aguilera and García-Muse, 2012). In this study, we show that the expression of the RNaseH1 protein, which resolves R-loop structures, provokes a significant decrease of DSBs induced upon G4-ligand treatments. In human cells, G4-forming sequences are highly correlated with R-loop-forming regions (Puget et al., 2019; Zheng et al., 2017; Zhang et al., 2020); therefore, our results suggest that some of the DSBs induced upon G4 stabilization by G4 ligands are associated with R-loop formation. Consistent with our data, De Margis et al. have recently reported that DNA damage induced by G4 ligands in human cells is mediated through R-loop formation (De Magis et al., 2019).
Why are G4 ligands toxic?
In cells, G4 ligands have been shown to affect cell growth through different mechanisms, altering telomere stability, replication, transcription, RNA metabolism, and mitochondrial maintenance (Varshney et al., 2020). Although cell alterations induced by these ligands have been linked to stabilization of G4 structures, their differential impact on different targets is probably related to their chemical structure, affecting their molecular activity and/or localization. For instance, the main impact of RHPS4 on mitochondrial DNA is related to its ability to translocate to the mitochondria, which is most probably due to its positive charge (Falabella et al., 2019). Our study, together with very recent reports (Bruno et al., 2020; Olivieri et al., 2020), indicates that cytotoxic effects induced by CX-5461 and PDS mainly rely on the formation of DNA breaks through a TOP2-dependent mechanism, arguing for the presence in these molecules of structural determinants that may block TOP2 activity during its active cycle. Supporting this assumption, competition experiments with 360A showed that stabilization of G4 structures is not sufficient to provoke TOP2-dependent DNA breaks. Of note, CX-5461, which was initially identified as a potent and selective RNA-Pol I inhibitor (Haddach et al., 2012), is a quarfloxin derivative, a small compound originally derived from the fluoroquinolone family that shows dual TOP2 inhibitor and G4-binding activities (Brooks and Hurley, 2010). Although the global TOP2 poisoning activity of quarfloxin derivatives has been shown to be attenuated by increased selectivity for G4 structures (Brooks and Hurley, 2010), we assume in our model that the TOP2 poisoning activity of CX-5461 is restricted to G4-forming regions.
In cells, inhibition of DNA-PKcs activity dramatically increases the number of DNA break signals in PDS and CX-5461-treated cells, demonstrating that an important number of G4 ligand-induced DNA breaks are repaired through the NHEJ pathway (Xu et al., 2017; McLuckie et al., 2013). This result suggests that most G4 ligand-induced breaks are efficiently repaired by NHEJ, while some TOP2-dependent DNA breaks refractory to repair probably cause the toxicity of these molecules.
On the basis of the results obtained in this study, and very recent reports from Bruno et al., 2020; Olivieri et al., 2020, we propose a model in which G4 ligands CX-5461 and PDS act as ‘G4-dependent TOP2 poisons’ (Figure 7). In this model, the interaction of both compounds with DNA is facilitated by DNA topological stress provoked by transcription. G4 stabilization by G4 ligands in transcriptionally active loci would provoke sustained RNA-Pol arrest, mobilizing topoisomerase enzymes to resolve topological stresses that at some loci may be poisoned in the vicinity of G4. Our model unifies the topoisomerase poisoning and G4-binding properties of these molecules in the new concept of DNA structure-driven topoisomerase poisoning at G-rich transcribed sequences.
Materials and methods
Cell lines
Request a detailed protocolAll the cell lines used in this study were obtained from ATCC that certified their identity. DAPI staining analysis was used to confirm the absence of Mycoplasma contamination during the course of experiments.
Cell culture conditions and treatments
Request a detailed protocolAll culture media were provided by Gibco and were supplemented with 10% fetal bovine serum (Eurobio), 100 U/mL penicillin (Gibco), and 100 µg/mL streptomycin (Gibco). Cells were grown in a humidified atmosphere with 5% CO2 at 37°C. HeLa and U2OS cells were grown in Dulbecco’s Modified Eagle Medium. HAP1 cells were cultured with Iscove's Modified Dulbecco's Medium. RPE1-hTERT cells were cultured with RPMI Media 1640 buffered with 0.3% Na(CO3)2. Expression of shTOPI in HeLa cells was induced with 5 µg/mL doxycycline for 5 days before treatments and for 4 days before siTOP2A transfections (described below) and maintained during siRNA-mediated protein depletion. RNaseHI-mCherry expression in U2OS cells was induced with 2.5 µg/mL doxycycline for 14 hr.
For selection of HAP1-resistant clones, haploid HAP1 were isolated using cell sorting and 100.106 haploid HAP1 were mutagenized by treatment with 300 µg/mL EMS (Sigma-Aldrich) for 3 days. After a 1-week recovery, 0.5–1.106 cells were seeded in 140 mm dishes. Plates were treated twice at a 1-week interval with 0.3 µM CX-5461 or 30 nM F14512 for 4 days. Around 10 days after the second treatment, individual clones (CXR and F14R) were isolated and used for further studies.
For immunofluorescence studies in U2OS, PDS was used at 20 µM for 8 hr. For immunofluorescence studies in HeLa and HAP1 cells, PDS and CX-5461 were used respectively at 20 μM and 0.2 µM for 4 hr. The G4 ligand 360A (20 µM) was added to cells 3 hr prior to PDS, CX-5461, and ETP (3.75 µM) treatments and renewed for the duration of the treatment. EdU treatments were performed at 10 μM at the same time as PDS and CX-5461 treatments and were maintained for the duration of the experiments. For inhibitors, DNA-PKi (2 μM, NU7441) and JQ1 (2 µM, (+/-)-JQ1) were added to cells 1 hr prior to PDS and or CX-5461 treatments. BNS22 (5 µM) was added to cells 30 min prior to PDS and or CX-5461 treatments. The transcription inhibitors DRB and BMH21 were used at 100 μM and 0.2 µM, respectively, and added to cells 1 hr prior to treatments. All inhibitors remained on cells for the duration of the experiment. For flow cytometry studies, HAP1 cells were incubated 1 hr with 10 µM EdU or exposed to 10 Gy X-ray irradiation (exposure time of 11 min 53 s) using a calibrated irradiation system (RX-650; Faxitron) fitted with a 0.5 mm aluminum filter for soft X-rays.
Chemical compound | Furnisher | CAS number |
---|---|---|
BMH21, N-[2-(Dimethylamino)ethyl]−12-oxo-2H-benzo[g]pyrido[2,1-b]quinazoline-4-carboxamide | Sigma-Aldrich | 896705-16-1 |
BNS-22 | Sigma-Aldrich | 1151668-24-4 |
Calicheamicin γ1 was generously gifted by P.R. Hamann (Wyeth Research, Pearl River, NY, USA) | ||
Camptothecin | Sigma-Aldrich | 7689-03-4 |
CX-5461 | Selleckchem | 1138549-36-6 |
DNA-PKi – NU7441 | Tocris | 503468-95-9 |
Doxycycline hydrochloride | Sigma-Aldrich | 10592-13-9 |
DRB, 5,6-dichloro-1-b-D-ribofuranosylbenzimidazole | Sigma-Aldrich | 54-85-0 |
EdU, 5-Ethynyl-2'-deoxyuridine | Thermo Fisher Scientific | A10044 |
Etoposide (ETP) | Sigma-Aldrich | 33419-42-0 |
F14512 | Pierre Fabre Laboratories | 866874-63-7 |
JQ1, (+/-)-JQ1 | Sigma-Aldrich | 1268524-69-1 |
Nocodazole | Sigma-Aldrich | 31430-18-9 |
PhenDC3 | Sigma-Aldrich | 929895-45-4 |
Pyridostatin | Sigma-Aldrich | 1085412-37-8 |
RHPS4 | Selleckchem | 390362-78-4 |
Plasmid and cell constructions
Request a detailed protocolU2OS cells conditionally expressing the RNaseH1-mCherry fusion protein were previously described in Britton et al., 2014. To construct HeLa cells conditionally expressing shRNAs against TOP1 mRNA, HeLa were infected with pLV-tTR-KRAB-Red and pLVTHM shTOP1 lentiviral particles. Individual clones from the transduced cell population were then isolated and selected for their capacity to downregulate TOP1 expression under treatment by the tetracycline analog doxycycline. pLV-tTR-KRAB-Red is a lentiviral vector encoding the transcriptional repressor tTR-KRAB fused to the DsRed fluorescent protein. pLVTHM is a lentiviral vector allowing conditional expression of an shRNA of interest under the control of the H1 promoter and the tetracycline operator/repressor system (TetO/TetR). pLVTHM vector allowing conditional expression of an shRNA against TOP1 was obtained by inserting duplex oligonucleotides (5′-CGCGTCCCCGGACTCCATCAGATACTATTTCAAGAGAATAGTATCTGATGGAGTCCTTTTTGGAAAT-3′ and 5′- CGATTTCCAAAAAGGACTCCATCAGATACTATTCTCTTGAAATAGTATCTGATGGAGTCCGGGGA-3′) between MluI and ClaI restriction sites in the pLVTHM plasmid. Transfection of HEK-293T cells (kindly provided by Genethon, Evry, France) with pLV-tTR-KRAB-Red or pLVTHM, and preparation of high titer lentiviruses pseudotyped with VSV-G protein have been performed as previously described (Delenda, 2004). pLV-tTRKRAB-red and pLVTHM were a gift from Didier Trono (Addgene plasmid # 12250; http://n2t.net/addgene:12250; RRID:Addgene_12250 and Addgene plasmid # 12247; http://n2t.net/addgene:12247; RRID:Addgene_12247) (Wiznerowicz and Trono, 2003).
RNA interference
Request a detailed protocolHeLa cells and HAP1 cells were respectively seeded at 250,000 cells and 400,000 cells per well in a 6-well plate. siRNAs oligonucleotides (see table below) were transfected twice (24 and 48 hr after seeding) at 50 nM final concentration per well with Lipofectamine RNAiMAX Reagent (Thermo Fisher Scientific) according to the manufacturer’s recommendations. For TOP2A and TOP2B co-depletion, each siRNA was used at the final concentration of 25 nM. Cells were split 24 hr after the second-round transfection for immunodetection, immunoblotting, and viability assays, and were treated 24 hr after being seeded.
Target | Name | Sequence or reference manufacturer | Manufacturer |
---|---|---|---|
Luciferase | siCtrl | 5′-CUUACGCUGAGUACUUCGATT-3′ | Eurofins Genomics |
TOP1 | siTOP1.2 | 5′-GGACUCCAUCAGAUACUAUTT-3′ | Eurofins Genomics |
TOP2A | siTOP2A | 5′-CGUAGGCUGUUUAAAGAAATT-3′ | Eurofins Genomics |
TOP2A | siTOP2A5 | 5′-GAAUAACCAUAGAAAUGAATT-3′ | Qiagen |
TOP2A | siTOP2A7 | 5′-GCGUGGUCAAAGAGUCAUUTT-3′ | Qiagen |
TOP2B | siTOP2B | 5′-GGGCUAGGAAAGAAGUAAATT-3′ | Qiagen |
RNA-seq
Request a detailed protocolRNA-seq was performed at the GeT-PlaGe core facility, INRA Toulouse, from total RNA prepared with the RNeasy Plus Mini Kit (Qiagen) according to the manufacturer’s instructions. RNA-seq libraries were prepared according to Illumina’s protocols using the Illumina TruSeq Stranded mRNA sample prep kit. Briefly, mRNAs were selected using poly-dT beads. Then, RNAs were fragmented and adaptors ligated. Eleven cycles of PCR were applied for library amplification. Library quality was assessed using a Fragment Analyser System (Agilent), and libraries were quantified by Q-PCR using the Kapa Library Quantification Kit (Roche). RNA-seq experiments were performed on an Illumina HiSeq3000 using a paired-end read length of 2 × 150 bp.
RNA-seq alignment and SNP prediction and filtering
Request a detailed protocolRead quality was checked within the ng6 environment (Mariette et al., 2012) using fastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and Burrows-Wheeler Aligner (BWA) (Li and Durbin, 2009) to search for contamination. The reads were cleaned with cutadapt v1.8.3 and aligned against hg38 reference human genome with STAR v2.5.2b (Dobin et al., 2013). Expression levels were computed with featureCount (Liao et al., 2014) using Ensembl annotation. Alignments were deduplicated with samtools rmdup and reads not uniquely mapped removed. Then GATK v3.5 base quality score recalibration was applied (McKenna et al., 2010). Indel realignment, SNP, and INDEL discovery were performed with HaplotypeCaller using standard hard filtering parameters according to GATK Best Practices recommendations for RNA-seq. Finally, variants were annotated using snpEff v4.3T (Cingolani et al., 2012). A Python script was used to select protein coding variants specific to CXR clones as compared to wild-type HAP1, with a minimal allele frequency of 0.9 and a depth greater than 10 reads. Among these variants, we selected variants resulting in frameshifts, mis- and non-sense mutations as compared to the reference human genome hg38. Cytoscape v3.2.0 (Shannon et al., 2003) was used to identify genes found mutated in several CXR clones. Upon TOP2A identification as a common gene mutated in five CXR clones, IGV v2.4.15 was used to scrutinize alignment data and revealed two TOP2A mutations missed by the analysis: for CXR #A2, the point mutation S654N; and, for CXR #A6, a mutation of the first nucleotide of the last intron leading to intron retention. Clone clustering under Cytoscape, based on shared mutated genes, suggested a common origin for clones CXR #A1, #A3, #A5, and #B4 (multiple common mutations). RNA-seq data from wild-type HAP1 and CXR clones have been deposited on SRA with the project ID PRJNA637883.
Targeted sequencing of TOP2A cDNA from HAP1 clones
Request a detailed protocolTotal RNAs were extracted from wild-type or F14R HAP1 with the RNeasy Plus Mini Kit (Qiagen) according to the manufacturer’s instructions. TOP2A cDNA was produced from these RNAs with the Superscript III First-Strand kit (Thermo Fisher Scientific) according to the manufacturer’s instructions and using the TOP2A-Rv primer. The resulting TOP2A cDNAs was amplified in four overlapping fragments using the primer pairs [TOP2A-F1, TOP2A-R1], [TOP2A-F2, TOP2A-R2], [TOP2A-F3, TOP2A-R3], and [TOP2A-F4, TOP2A-Rv] and sequenced using the same primers except for the last fragment for which the TOP2A-R4 sequencing primer was also used.
Name | Oligonucleotide sequence (5′ to 3′) |
---|---|
TOP2A-F1 | GTCGCTTTCAGGGTTCTTGAGCC |
TOP2A-R1 | TGGCATGTTGATCCAAAGCTCTTGG |
TOP2A-F2 | TGGTGTTGCAGTAAAAGCACATCAGG |
TOP2A-R2 | GCAACCTTTACTTCTCGCTTGTCATTCC |
TOP2A-F3 | TCCTGAGGATTACTTGTATGGACAAACTACC |
TOP2A-R3 | GCCTTCACAGGATCCGAATCATATCCC |
TOP2A-F4 | GGCTCCTAGGAATGCTTGGTGC |
TOP2A-R4 | TCATCTGGGAAATGTGTAGCAGGAGG |
TOP2A-Rv | GCTTCAGGTAACTTTAAAACCAGTCTTGG |
Heparin-based extraction for TOP2cc immunodetection
Request a detailed protocolThis method was adapted from de Campos-Nebel et al., 2016. For immunofluorescence, HeLa cells were seeded as described below. Cells were washed with ice-cold PBS and incubated twice for 5 min on ice with CSK buffer (10 mM PIPES pH 7.0, 100 mM NaCl, 300 mM sucrose, 3 mM MgCl2) containing 0.7% Triton X-100 and 20 U/mL Heparin (Sigma-Aldrich). Cells were washed with ice-cold PBS and fixed on ice 15 min with paraformaldehyde 2% in PBS. Then, immunofluorescence was performed as described below without the permeabilization step. For immunoblotting, HAP1 cells were seeded at 1.5 106 cells in 6 cm dishes 24 hr prior ETP treatment (200 µM, 1 hr). After treatment, cells were harvested with trypsin and washed with cold PBS. After gentle centrifugation, cells were resuspended in lysis buffer (150 mM NaCl, 1 mM EDTA, 0.5% IGEPAL CA-630, 2X HALT Protease and Phosphatase Inhibitor Cocktail [Thermo Fisher Scientific] 20 mM Tris-HCl, pH 8.0) complemented with 100 U/mL Heparin (H3393, Sigma-Aldrich) and incubated on ice for 15 min. Then lysates were centrifugated at 15,000 rpm at 4°C for 5 min and pellets were resuspended with lysis buffer. In order to facilitate migration on polyacrylamide gel, a sonication was performed to degrade DNA present within the extracts. Protein concentrations were determined by measuring absorbance at 280 nm (Nanodrop) and heparin-based extracts were diluted with denaturing lysis buffer (120 mM Tris-HCl pH 6.8, 4% SDS, 20% glycerol). Western blot was performed as described below.
Cell lysis and western blotting
Request a detailed protocolWhole-cell extracts were prepared from PBS-washed pellet lysed with denaturing lysis buffer (120 mM Tris-HCl pH 6.8, 4% SDS, 20% glycerol) and 10 strokes through a 24G needle. Protein concentrations were determined by measuring absorbance at 280 nm (Nanodrop). For loading, an equal volume of a solution of 0.01% bromophenol blue and 200 mM dithiothreitol was added to the extracts, then boiled at 95°C for 5 min. About 80 μg of denatured proteins were loaded for each condition and separated on standard or Stain-Free gradient 4–12% polyacrylamide TGX pre-cast gels (Bio-Rad) and transferred onto nitrocellulose membranes (0.45 µm pore size, Bio-Rad or Protran, GE Healthcare). Before blocking (incubation 1 hr at room temperature in PBS 0.1% Tween-20, 5% non-fat dry milk), Ponceau S staining or UV exposition of membrane (for Stain-Free gels) was used to confirm homogeneous loading. The membrane was successively probed with primary antibody and appropriate goat secondary antibodies coupled to horseradish peroxidase (described in table below). A ChemiDoc CCD imager (Bio-Rad) was used to acquire pictures of the stain-free total protein staining and the chemiluminescence signal after membrane incubation with adequate peroxidase substrate (Clarity ECL, Bio-Rad). Digital data were processed and quantified using ImageJ software.
Target | Dilution | Species | Class | Reference | Manufacturer |
---|---|---|---|---|---|
αTubulin | 1:25,000 | Mouse | Monoclonal | T5168 | Sigma-Aldrich |
KU80 (Dobin et al., 2013) | 0.2 µg/mL | Mouse | Monoclonal | MA5-12933 | Thermo Fisher Scientific |
KU70 (N3H10) | 0.2 µg/mL | Mouse | Monoclonal | MA5-13110 | Thermo Fisher Scientific |
TOP1 | 1:1000 | Rabbit | Monoclonal | EPR5375 | Abcam |
TOP2A (aa1352-1493) | 1 µg/mL | Mouse | Monoclonal | GTX35137 | Genetex |
TOP2A (aa1500-C-term) | 1 µg/mL | Rabbit | Polyclonal | A300-054B | Bethyl Laboratories |
TOP2B | 0.2 µg/mL | Rabbit | Polyclonal | A300-950A | Bethyl Laboratories |
Anti-rabbit HRP-coupled | 1:10,000 | Goat | Polyclonal | 111-035-003 | Jackson Immunoresearch |
Anti-mouse HRP-coupled | 1:10,000 | Goat | Polyclonal | 115-035-003 | Jackson Immunoresearch |
Immunofluorescence
Request a detailed protocolRPE1-hTERT cells were seeded, treated, and stained in the same conditions than HeLa cells. HeLa cells and HAP1 cells were seeded in 24-wells plate at respectively 100,000 cells/well and 25,000 cells/well on #1.5 glass coverslips (VWR, #631-0150). HeLa cells and HAP1 cells were respectively treated 24 hr and 48 hr later, and then fixed with paraformaldehyde 2% in PBS at room temperature (10 min for HeLa cells, 15 min for HAP1 cells), washed with PBS, and permeabilized for 15 min at room temperature with 10 mM Tris-HCl pH 7.5, 120 mM KCl, 20 mM NaCl, 0.1% Triton-X 100. In EdU- treated cells, cells were washed with PBS and EdU detection reaction was performed with Click-iT RNA Alexa Fluor Imaging Kit according to the manufacturer’s recommendation, with 2 µM Alexa Fluor 594 or 648 azide for 30 min at room temperature. Then, cells were washed with PBS and incubated for about 1 hr at room temperature in blocking buffer (20 mM Tris-HCl pH 7.5, 150 mM NaCl, 2% BSA, 0.2% fish gelatin, 0.1% Triton-X 100) prior to incubation overnight at 4°C with primary antibody diluted in blocking buffer (dilutions shown in table below). For BG4 immunodetection, blocking buffer were complemented 0.3 μg/μL of RNAse A (David et al., 2019). Cells were then washed with PBS 0.1% Tween-20 and incubated with appropriate secondary goat antibody coupled to Alexa Fluor 488 or 594 diluted in blocking buffer (dilutions shown in table below) for 1 hr at room temperature. For TOP2Acc and BG4 co-immunodetection, both primary anti-BG4 antibody and appropriate secondary antibody were successively added to cells prior to incubations to primary anti-TOP2A antibody and its appropriate secondary antibody. At last, cells were washed with PBS 0.1% Tween-20 and stained with 0.1 μg/mL DAPI for 20 min at room temperature, and coverslips were mounted with Vectashield mounting medium (Vector Laboratories). Nuclear γH2AX foci, 53BP1 foci, BG4 foci, TOP2Acc foci, and EdU-integrated density staining overlapping with DAPI staining were quantified with ImageJ software. Nuclear DAPI-integrated density staining was quantified with ImageJ software and correlated to nuclear EdU-integrated density to determine cell cycle phase for each cell as described in Roukos et al., 2015. Quantifications of nuclear γH2AX, 53BP1, and TOP2Acc foci induced by G4 ligands or ETP are represented normalized to non-treated (NT) conditions. For TOP2Acc imaging, images were acquired with a Zeiss Elyra 7 3D Lattice SIM super-resolution microscope fitted with a 63× objective (PLANAPO NA 1.4, Zeiss) and dual sCMOS cameras (pco.edge). 3D-SIM reconstructions were performed with Zen Black 2.3 (Zeiss). For TOP2Acc and BG4 co-immunodetection, images were obtained by performing a maximum intensity projection of 20 3D-SIM Z-stacks (interval 0.091 µM) with Zen Blue 3.3 (Zeiss). Quantification of foci and colocalization events was done manually. Zen Blue 3.3 was used to adjust brightness and contrast of corresponding micrographs as well as image cropping.
Target | Dilution (µg/mL) | Species | Class | Reference | Manufacturer |
---|---|---|---|---|---|
53BP1 | 1.3 | Mouse | Monoclonal | MAB-3803 | Millipore |
γH2AX (Phospho S139) | 0.7 | Rabbit | Monoclonal | 81299 | Abcam |
BG4 | 0.25 | Mouse | Monoclonal | Ab00174-1.1 | Absolute antibody |
TOP2A | 0.5 | Rabbit | Polyclonal | A300-054B | Bethyl Laboratories |
TOP2A | 0.2 | Mouse | Monoclonal | GTX35137 | Genetex |
Anti-rabbit 488 | 2 | Goat | Polyclonal | A11008 | Thermo Fisher Scientific |
Anti-mouse 594 | 2 | Goat | Polyclonal | A11005 | Thermo Fisher Scientific |
Anti-mouse 488 | 2 | Goat | Polyclonal | A11001 | Thermo Fisher Scientific |
Anti-rabbit 594 | 2 | Goat | Polyclonal | A11012 | Thermo Fisher Scientific |
Viability assay (SRB)
Request a detailed protocolHAP1 cells were seeded in 96-flat-wells plate at 3500 cells per well. Serial dilutions of various compounds were realized allowing same solvent concentration for each condition, and cells were treated 24 hr after seeding. After 3 days, HAP1 cells were fixed for 1 hr at 4°C by addition of 10% trichloroacetic acid to a 3.33% final concentration, before being washed with tap water and dried overnight. Cells were stained by incubation 30 min at room temperature in a 1% acetic acid solution containing 0.057% sulforhodamin B, then cells were washed with 1% acetic acid and dried overnight. Finally, 200 μL of a 10 mM Tris-base solution was added, plates were agitated for 1 hr at room temperature, and SRB levels were measured by absorbance at 490 nm using μQuant microplate spectrophotometer (Bio-Tek Instruments). Percentages of cell viability are expressed after normalization relative to NT controls. For characterization of wild-type and drug-resistant HAP1, the IC50 (50% inhibitory concentration) was computed for each drug and cell lines with the GraphPad Prism v8 software using a nonlinear regression to a four-parameter logistic curve.
Cell proliferation assay
Request a detailed protocolHAP1 cells were seeded in 6-well plates at 40,000 cells/well. The occupied area (% confluency) was monitored every hour for 74 hr using a live-imaging Incucyte Zoom system (Sartorius). Analysis of the 7–72 hr exponential part of the occupied area vs. time exponential curve was used to compute the doubling time of wild-type and resistant HAP1 cells using GraphPad Prism v8. The analysis was repeated three times, and the scatter dot-blot shows the three values, the mean, and the SD from the mean. A one-way ANOVA test was used and revealed no significant difference between the doubling time of resistant clones and the WT HAP1.
Clonogenic assay
Request a detailed protocolClonogenic assay was performed as described by Bombarde et al., 2017. Briefly, after transfection with siRNA, HeLa cells was seeded at low density (250 cells/well) the day before treatment, pre-incubated with 100 µM DRB, 0.2 µM BMH21, or dimethylsulfoxide (DMSO) for 1 hr and treated for 4 hr with PDS or CX-5461 in the presence of transcription inhibitor or DMSO before being replaced by fresh medium. After 10–15 days, cells were stained with crystal violet and the colonies were counted (at least 50 colonies were counted for each condition per experiment). Data were normalized to the NT conditions. The IC50 was computed in each condition with the GraphPad Prism v8 software using a nonlinear regression to a four-parameter logistic curve.
Flow cytometry
Request a detailed protocolBriefly, HAP1 cells were collected by trypsination at the end of the treatment, washed with PBS 1% BSA, and fixed for 15 min with paraformaldehyde 2% in PBS at room temperature. Cells were washed with PBS 1% BSA, incubated for 30 min with PBS 0.2% Triton X-100, and washed again with PBS 1% BSA. If necessary, EdU detection reaction was performed as described above. If not, cells were incubated for 1 hr with primary anti-γH2AX antibody diluted in PBS 0.1% Tween-20 5% BSA (1:1000, ref ZMS05636, Sigma-Aldrich), washed with PBS 1% BSA, and incubated 30 min with PBS 0.1% Tween-20 5% BSA containing appropriate secondary antibodies coupled to Alexa Fluor 488 (1:200). Lastly, cells were washed with PBS 1% BSA and incubated in PBS containing 250 µg/mL RNAse A and 2 µg/mL DAPI. A BD LSR II flow cytometer (Becton Dickinson) was used to analyze a minimum of 30,000 cells. Data were analyzed and formatted using FlowJo v8.8.7.
Statistical analyses
Request a detailed protocolAll results represent at least three independent experiments. Statistical analyses were performed with the GraphPad Prism software (version 8). For γH2AX and 53BP1 quantifications analyses, multiple unpaired t tests (without corrections for multiple comparisons) were performed between pairs of conditions. For BG4 quantifications analyses, results of at least three independent experiments were pulled together and unpaired Welch’s t tests were performed between pairs of conditions. On all figures, significant differences between specified pairs of conditions are shown by asterisks (*p-value<0.05; **p-value<0.01; ***p-value<0.001; ****p-value<0.0001). NS is for nonsignificant difference.
Data availability
RNA-seq data from wild-type HAP1 and CXR clones have been deposited on SRA with the project ID PRJNA637883.
-
NCBI BioProjectID PRJNA637883. RNA-seq sequencing data from individual CX-5461-resistant HAP-1 clones.
References
-
The connection between transcription and genomic instabilityThe EMBO Journal 21:195–201.https://doi.org/10.1093/emboj/21.3.195
-
MRN, CtIP, and BRCA1 mediate repair of topoisomerase II-DNA adductsJournal of Cell Biology 212:399–408.https://doi.org/10.1083/jcb.201504005
-
Topoisomerase-mediated chromosomal break repair: an emerging player in many gamesNature Reviews Cancer 15:137–151.https://doi.org/10.1038/nrc3892
-
TOP2B: the first thirty yearsInternational Journal of Molecular Sciences 19:2765.https://doi.org/10.3390/ijms19092765
-
DNA damage triggers SAF-A and RNA biogenesis factors exclusion from chromatin coupled to R-loops removalNucleic Acids Research 42:9047–9062.https://doi.org/10.1093/nar/gku601
-
Targeting MYC Expression through G-QuadruplexesGenes & Cancer 1:641–649.https://doi.org/10.1177/1947601910377493
-
Quadruplex DNA: sequence, topology and structureNucleic Acids Research 34:5402–5415.https://doi.org/10.1093/nar/gkl655
-
Construction, characterization, and complementation of a conditional-lethal DNA topoisomerase IIalpha mutant human cell lineMolecular Biology of the Cell 15:5700–5711.https://doi.org/10.1091/mbc.e04-08-0732
-
High-throughput sequencing of DNA G-quadruplex structures in the human genomeNature Biotechnology 33:877–881.https://doi.org/10.1038/nbt.3295
-
Emerging roles for R-loop structures in the management of topological stressJournal of Biological Chemistry 295:4684–4695.https://doi.org/10.1074/jbc.REV119.006364
-
R-Loops as cellular regulators and genomic threatsMolecular Cell 73:398–411.https://doi.org/10.1016/j.molcel.2019.01.024
-
CNBP controls transcription by unfolding DNA G-quadruplex structuresNucleic Acids Research 47:7901–7913.https://doi.org/10.1093/nar/gkz527
-
Lentiviral vectors: optimization of packaging, transduction and gene expressionThe Journal of Gene Medicine 6:S125–S138.https://doi.org/10.1002/jgm.501
-
The DNA cleavage reaction of topoisomerase II: wolf in sheep's clothingNucleic Acids Research 37:738–748.https://doi.org/10.1093/nar/gkn937
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
G-Quadruplexes in RNA biology: recent advances and future directionsTrends in Biochemical Sciences 46:270–283.https://doi.org/10.1016/j.tibs.2020.11.001
-
G4 motifs correlate with promoter-proximal transcriptional pausing in human genesNucleic Acids Research 39:4975–4983.https://doi.org/10.1093/nar/gkr079
-
Genome-wide genetic screening with chemically mutagenized haploid embryonic stem cellsNature Chemical Biology 13:12–14.https://doi.org/10.1038/nchembio.2226
-
Site-specific cleavage of a DNA hairpin by topoisomerase II. DNA secondary structure as a determinant of enzyme recognition/cleavageJournal of Biological Chemistry 269:7719–7725.https://doi.org/10.1016/S0021-9258(17)37346-5
-
Transcription-associated genome instabilityChemical Reviews 113:8638–8661.https://doi.org/10.1021/cr400017y
-
Transcription-mediated replication hindrance: a major driver of genome instabilityGenes & Development 33:1008–1026.https://doi.org/10.1101/gad.324517.119
-
Discovery of CX-5461, the first direct and selective inhibitor of RNA polymerase I, for Cancer therapeuticsACS Medicinal Chemistry Letters 3:602–606.https://doi.org/10.1021/ml300110s
-
G-quadruplex structures mark human regulatory chromatinNature Genetics 48:1267–1272.https://doi.org/10.1038/ng.3662
-
DNA G-quadruplexes in the human genome: detection, functions and therapeutic potentialNature Reviews Molecular Cell Biology 18:279–284.https://doi.org/10.1038/nrm.2017.3
-
DrugTargetSeqR: a genomics- and CRISPR-Cas9-based method to analyze drug targetsNature Chemical Biology 10:626–628.https://doi.org/10.1038/nchembio.1551
-
The interplay between G-quadruplex and transcriptionCurrent Medicinal Chemistry 26:2898–2917.https://doi.org/10.2174/0929867325666171229132619
-
The roles of DNA topoisomerase iiβ in transcriptionInternational Journal of Molecular Sciences 19:1917.https://doi.org/10.3390/ijms19071917
-
26 S proteasome-mediated degradation of topoisomerase II cleavable complexesJournal of Biological Chemistry 276:40652–40658.https://doi.org/10.1074/jbc.M104009200
-
G-quadruplex DNA as a molecular target for induced synthetic lethality in cancer cellsJournal of the American Chemical Society 135:9640–9643.https://doi.org/10.1021/ja404868t
-
Induction of apoptosis by deregulated expression of DNA topoisomerase IIalphaCancer Research 58:4519–4524.
-
G-quadruplex-R-loop interactions and the mechanism of anticancer G-quadruplex bindersNucleic Acids Research 48:11942–11957.https://doi.org/10.1093/nar/gkaa944
-
Site-Specific cleavage by topoisomerase 2: a mark of the core centromereInternational Journal of Molecular Sciences 19:534.https://doi.org/10.3390/ijms19020534
-
Existence and consequences of G-quadruplex structures in DNACurrent Opinion in Genetics & Development 25:22–29.https://doi.org/10.1016/j.gde.2013.10.012
-
DNA topoisomerase II and its growing repertoire of biological functionsNature Reviews Cancer 9:327–337.https://doi.org/10.1038/nrc2608
-
Targeting DNA topoisomerase II in Cancer chemotherapyNature Reviews Cancer 9:338–350.https://doi.org/10.1038/nrc2607
-
Nonhomologous DNA end-joining for repair of DNA double-strand breaksJournal of Biological Chemistry 293:10512–10523.https://doi.org/10.1074/jbc.TM117.000374
-
Small molecule BMH-compounds that inhibit RNA polymerase I and cause nucleolar stressMolecular Cancer Therapeutics 13:2537–2546.https://doi.org/10.1158/1535-7163.MCT-14-0256
-
Roles of eukaryotic topoisomerases in transcription, replication and genomic stabilityNature Reviews Molecular Cell Biology 17:703–721.https://doi.org/10.1038/nrm.2016.111
-
A novel small molecule that alters shelterin integrity and triggers a DNA-damage response at telomeresJournal of the American Chemical Society 130:15758–15759.https://doi.org/10.1021/ja805615w
-
Small-molecule-induced DNA damage identifies alternative DNA structures in human genesNature Chemical Biology 8:301–310.https://doi.org/10.1038/nchembio.780
-
Cell cycle staging of individual cells by fluorescence microscopyNature Protocols 10:334–348.https://doi.org/10.1038/nprot.2015.016
-
Telomere damage induced by the G-quadruplex ligand RHPS4 has an antitumor effectJournal of Clinical Investigation 117:3236–3247.https://doi.org/10.1172/JCI32461
-
The effects of DNA supercoiling on G-quadruplex formationNucleic Acids Research 45:12069–12079.https://doi.org/10.1093/nar/gkx856
-
rDNA chromatin activity status as a biomarker of sensitivity to the RNA polymerase I transcription inhibitor CX-5461Frontiers in Cell and Developmental Biology 8:568.https://doi.org/10.3389/fcell.2020.00568
-
In vitro footprinting of promoter regions within supercoiled plasmid DNAMethods in Molecular Biology 613:223–233.https://doi.org/10.1007/978-1-60327-418-0_14
-
Topoisomerase II contributes to DNA secondary structure-mediated double-stranded breaksNucleic Acids Research 48:6654–6671.https://doi.org/10.1093/nar/gkaa483
-
The regulation and functions of DNA and RNA G-quadruplexesNature Reviews Molecular Cell Biology 21:459–474.https://doi.org/10.1038/s41580-020-0236-x
-
Using transcriptome sequencing to identify mechanisms of drug action and resistanceNature Chemical Biology 8:235–237.https://doi.org/10.1038/nchembio.779
-
DNA topoisomerase II as a target for Cancer chemotherapyCancer Investigation 20:570–589.https://doi.org/10.1081/CNV-120002156
-
Human DNA topoisomerase IIbeta binds and Cleaves four-way junction DNA in vitroNucleic Acids Research 27:984–992.https://doi.org/10.1093/nar/27.4.984
-
A protease pathway for the repair of topoisomerase II-DNA covalent complexesJournal of Biological Chemistry 281:35997–36003.https://doi.org/10.1074/jbc.M604149200
Article and author information
Author details
Funding
Agence Nationale de la Recherche (ANR-17-CE18-0002-01)
- Madeleine Bossaert
- Linh-Trang Nguyên
- Patrick Calsou
- Sébastien Britton
- Dennis Gomez
Agence Nationale de la Recherche (ANR-16-CE11-0006-01)
- Angélique Pipier
- Jean-Francois Riou
- Eric Defrancq
- Patrick Calsou
- Sébastien Britton
- Dennis Gomez
Cancéropôle Grand Ouest (Emergence funding "CX-Break")
- Madeleine Bossaert
- Angélique Pipier
- Linh-Trang Nguyên
- Patrick Calsou
- Sébastien Britton
- Dennis Gomez
Ligue Contre le Cancer (Equipe Labellisée 2018)
- Madeleine Bossaert
- Angélique Pipier
- Linh-Trang Nguyên
- Patrick Calsou
- Sébastien Britton
- Dennis Gomez
Agence Nationale de la Recherche (ANR-10-INBS-09)
- Céline Noirot
- Remy-Felix Serre
- Olivier Bouchez
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was funded by grants from ANR (ANR-17-CE18-0002-01 and ANR-16-CE11-0006-01), Cancéropôle GSO Emergence funding ‘CX-Break,’ and La Ligue Nationale Contre le Cancer (Equipe Labellisée 2018). We thank Julia Coates and Steve Coates for proofreading the manuscript. We are grateful to Emmanuelle Näser, Antonio Peixoto, and Elodie Vega for technical assistance and the Imaging Core Facility TRI-IPBS, supported by ITMO Cancer (Alliance Nationale pour les Sciences de la Vie et de la Santé, National Alliance for Life Sciences and Health) within the Framework of the Cancer Plan. We thank P Mailliet for providing 360A. We also thank J Coates and S Coates for the careful proofreading and language editing of the manuscript. We thank ‘Région Midi‐Pyrénées’ for supporting the ‘Toulouse Réseaux Imagerie’ platform. Patrick Calsou is a researcher from INSERM. This work was performed in collaboration with the GeT core facility, Toulouse, France (http://get.genotoul.fr), and was supported by France Génomique National infrastructure, funded as part of ‘Investissement d’avenir’ program managed by Agence Nationale pour la Recherche (contract ANR-10-INBS-09).
Copyright
© 2021, Bossaert et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,333
- views
-
- 602
- downloads
-
- 52
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Cancer Biology
Many cells in high glucose repress mitochondrial respiration, as observed in the Crabtree and Warburg effects. Our understanding of biochemical constraints for mitochondrial activation is limited. Using a Saccharomyces cerevisiae screen, we identified the conserved deubiquitinase Ubp3 (Usp10), as necessary for mitochondrial repression. Ubp3 mutants have increased mitochondrial activity despite abundant glucose, along with decreased glycolytic enzymes, and a rewired glucose metabolic network with increased trehalose production. Utilizing ∆ubp3 cells, along with orthogonal approaches, we establish that the high glycolytic flux in glucose continuously consumes free Pi. This restricts mitochondrial access to inorganic phosphate (Pi), and prevents mitochondrial activation. Contrastingly, rewired glucose metabolism with enhanced trehalose production and reduced GAPDH (as in ∆ubp3 cells) restores Pi. This collectively results in increased mitochondrial Pi and derepression, while restricting mitochondrial Pi transport prevents activation. We therefore suggest that glycolytic flux-dependent intracellular Pi budgeting is a key constraint for mitochondrial repression.
-
- Cancer Biology
- Genetics and Genomics
Uveal melanoma (UM) is a rare melanoma originating in the eye’s uvea, with 50% of patients experiencing metastasis predominantly in the liver. In contrast to cutaneous melanoma, there is only a limited effectiveness of combined immune checkpoint therapies, and half of patients with uveal melanoma metastases succumb to disease within 2 years. This study aimed to provide a path toward enhancing immunotherapy efficacy by identifying and functionally validating tumor-reactive T cells in liver metastases of patients with UM. We employed single-cell RNA-seq of biopsies and tumor-infiltrating lymphocytes (TILs) to identify potential tumor-reactive T cells. Patient-derived xenograft (PDX) models of UM metastases were created from patients, and tumor sphere cultures were generated from these models for co-culture with autologous or MART1-specific HLA-matched allogenic TILs. Activated T cells were subjected to TCR-seq, and the TCRs were matched to those found in single-cell sequencing data from biopsies, expanded TILs, and in livers or spleens of PDX models injected with TILs. Our findings revealed that tumor-reactive T cells resided not only among activated and exhausted subsets of T cells, but also in a subset of cytotoxic effector cells. In conclusion, combining single-cell sequencing and functional analysis provides valuable insights into which T cells in UM may be useful for cell therapy amplification and marker selection.