Beta human papillomavirus 8E6 promotes alternative end joining
Abstract
Double strand breaks (DSBs) are one of the most lethal DNA lesions in cells. The E6 protein of beta-human papillomavirus (HPV8 E6) impairs two critical DSB repair pathways: homologous recombination (HR) and non-homologous end joining (NHEJ). However, HPV8 E6 only delays DSB repair. How DSBs are repaired in cells with HPV8 E6 remains to be studied. We hypothesize that HPV8 E6 promotes a less commonly used DSB repair pathway, alternative end joining (Alt-EJ). Using CAS9-based Alt-EJ reporters, we show that HPV8 E6 promotes Alt-EJ. Further, using small molecule inhibitors, CRISPR/CAS9 gene knockout, and HPV8 E6 mutant, we find that HPV8 E6 promotes Alt-EJ by binding p300, an acetyltransferase that facilitates DSB repair by HR and NHEJ. At least some of this repair occurs through a subset of Alt-EJ known as polymerase theta dependent end joining. Finally, whole genome sequencing analysis showed HPV8 E6 caused an increased frequency of deletions bearing the microhomology signatures of Alt-EJ. This study fills the knowledge gap of how DSB is repaired in cells with HPV8 E6 and the mutagenic consequences of HPV8 E6 mediated p300 destabilization. Broadly, this study supports the hypothesis that beta-HPV promotes cancer formation by increasing genomic instability.
Editor's evaluation
This article reports useful data on how human papillomavirus 8E6 protein regulates DSB repair pathways in human cells. The data support the claim that 8E6 promotes alternative end-joining through binding and destabilizing the p300 acetyltransferase, showing the involvement of PARP-1-dependent alternative end-joining and to a lesser degree DNA Polymerase theta-dependent alternative end-joining.
https://doi.org/10.7554/eLife.81923.sa0Introduction
Beta genus human papillomaviruses (beta-HPVs) are ubiquitous and transiently infect cutaneous epithelia in the general population (Bouwes Bavinck et al., 2018; Gheit, 2019; Neale et al., 2013). Beta-HPVs, including type 8 (HPV8), are associated with nonmelanoma skin cancer (NMSC) in immunocompromised individuals including people with a rare genetic disorder epidermodysplasia verruciformis and organ transplant recipients (Bouwes Bavinck et al., 2018; Bouwes Bavinck et al., 2001; Dell’Oste et al., 2009; Sichero et al., 2019). However, the contribution of beta-HPV infections to NMSC in the general population is unclear. Because beta-HPV infections are not persistent in immunocompetent people, they are hypothesized to promote cancer formation by making UV-induced DNA damage more mutagenic (Hufbauer et al., 2015; Spriggs and Laimins, 2017; Wendel and Wallace, 2017). In support of this hypothesis, our group and others have used in vitro and in vivo systems to demonstrate that the E6 protein from HPV8 (8E6) impairs DNA repair (Hufbauer et al., 2015; Rollison et al., 2019; Wallace et al., 2012).
The interaction with and destabilization of p300 is a key mechanism by which 8E6 hinders DNA repair (Dacus and Wallace, 2021; Howie et al., 2011). P300 is an acetyltransferase that regulates transcription by chromatin remodeling (Ait-Si-Ali et al., 2000; Dutto et al., 2018; Goodman and Smolik, 2000). By binding p300, 8E6 decreases the abundance of multiple DNA repair proteins including ATR, ATM, BRCA1, and BRCA2 (Wallace et al., 2015; Wallace et al., 2013; Wallace et al., 2012). This lowers the activation of ATM and ATR signaling, decreasing the cellular response to UV-damaged DNA (Snow et al., 2019; Wallace et al., 2012). The limited ability to repair UV damage increases the frequency with which UV causes replication forks to collapse into double strand breaks (DSBs) in DNA (Pathania et al., 2011; Rapp and Greulich, 2004). Cells have multiple DSB repair mechanisms. Homologous recombination (HR) is minimally mutagenic, but restricted to in S/G2 phase when the sister chromatids can serve as homologous templates (Orthwein et al., 2015; Shibata et al., 2011). Whenever possible cells use HR to fix DSBs as it allows them to avoid mutations. When HR is inhibited (by cell cycle position, mutation to repair factors, or artificially), non-homologous end joining (NHEJ) is used to repair DSBs (Wang et al., 2018). NHEJ can occur throughout the cell cycle, as it does not require a homologous template. However, because NHEJ generates and ligates blunt ends to fix a DSB, it is more mutagenic than HR (Mao et al., 2008). We have shown that 8E6 attenuates both HR and NHEJ by degrading p300 (Hu et al., 2020; Hu and Wallace, 2022c; Wallace et al., 2015).
Notably, 8E6 does not block the initiation of NHEJ or HR. NHEJ initiation occurs when DNA-dependent protein kinase catalytic subunit (DNA-PKcs) complexes at DSBs and activates itself by autophosphorylation (pDNA-PKcs) (Davis et al., 2014; Jiang et al., 2015). This occurs readily in the presence of 8E6. However, 8E6 prevents the resolution of pDNA-PKcs repair complexes and attenuates other downstream steps in NHEJ (Hu et al., 2020). Similarly, 8E6 allows the HR pathway to initiate, before hindering the resolution of RAD51 repair complexes (Wallace et al., 2015). We recently demonstrated that cells respond to 8E6-associated inhibition of NHEJ by trying to complete HR during G1 (Hu et al., 2022a; Hu et al., 2022b). This ultimately leads to persistent unresolved RAD51 repair complexes formed during G1.
Thus, currently there is a detailed understanding of how 8E6 causes DSB repair to fail, but less is known about how DSB repair occurs in cells that express 8E6. When NHEJ and HR fail, another mutagenic repair pathway known as alternative end joining (Alt-EJ) is tasked with completing DSB repair (Iliakis et al., 2015). Here, we use reporter constructs and small molecule inhibitors of DNA repair factors to demonstrate that 8E6 promotes DSB repair by Alt-EJ and to show that the use of Alt-EJ is the indirect result of initiating NHEJ when the pathway cannot be completed. We also employ whole genome sequence analysis of cells expressing 8E6 and passage matched empty vector control expressing cells to determine the frequency with which 8E6 promotes mutations with the characteristics of repair by Alt-EJ. These observations address a key knowledge gap in the field. By promoting DSB repair by Alt-EJ, 8E6 increases the risk of mutations associated with DSBs while allowing cells to avoid the apoptosis that would be associated with an unrepaired DSB. This is consistent with the proposed mechanism by which beta-HPV infections are hypothesized to promote NMSC.
Results
8E6 promotes DSB repair by Alt-EJ
8E6 delays rather than abrogates DSB repair, but it is unclear how these lesions are repaired as 8E6 hinders the completion of HR and NHEJ. Because Alt-EJ serves as another DSB repair mechanism should HR and NHEJ fail, we hypothesized that DSBs were instead repaired by Alt-EJ. To test this, we examined previously described telomerase (N/TERT) immortalized human foreskin keratinocyte (HFK) expressing vector control (HFK LXSN) and 8E6 (HFK 8E6) (Bedard et al., 2008). An established reporter cassette where a 46 nt insertion disrupts a GFP open reading frame was used to measure Alt-EJ (Bhargava et al., 2018; Tsai et al., 2020). Figure 1A describes Alt-EJ that requires resection (imbedded) and Alt-EJ that occurs independent of resection (terminal). The end joining is mediated by a 4 nt microhomology (ACGG). Transient transfection of this reporter into HFK LXSN and HFK 8E6 demonstrated that 8E6 increased terminal and imbedded Alt-EJ (Figure 1B–E and Figure 1—figure supplement 1). Transfection efficiency varies between independent experiments, but no significant difference was observed among groups within each experiment (Figure 1—figure supplement 1).
To further confirm that 8E6 promoted Alt-EJ, we examined DSB repair in cells where Alt-EJ was blocked by a small molecule inhibitor against PARP-1, an established component of the Alt-EJ pathway (1 µM of olaparib) (Kötter et al., 2014). We confirmed that PARP-1 inhibition blocked Alt-EJ using the reporter system described above (Figure 2A–D and Figure 2—figure supplement 1). DSBs were induced by growth in media containing zeocin (10 µg/mL, 10 min), a radiation mimetic reagent (Tsukuda and Miyazaki, 2013). H2AX phosphorylated at Serine 139 or pH2AX is used as a standard DSB marker (Rogakou et al., 1998). PARP-1 inhibition did not significantly alter DSB repair in HFK LXSN cells (Figure 2E–F). This is consistent with the established view that most DSB repair occurs via either HR or NHEJ, with Alt-EJ serving as a fail-safe machinery should these pathways fail (Iliakis et al., 2015; Patterson-Fortin and D’Andrea, 2020; Simsek and Jasin, 2010). In contrast, PARP-1 inhibition significantly delayed DSB repair in HFK 8E6 cells (Figure 2E–F).
Alt-EJ is an umbrella term that includes any end joining event that does not require DNA-PKcs activity. While most Alt-EJ relies on PARP-1, a subset of Alt-EJ can occur when PARP-1 activity is impaired. Specifically, loss of PARP-1 activity only reduces polymerase theta mediated end joining (TMEJ) by two- to fourfold (Luedeman et al., 2022). To begin determining the extent that 8E6 promoted TMEJ, we used immunoblots to determine if 8E6 increased the abundance of polymerase theta (POLθ). This analysis did not find any notable changes in POLθ abundance (Figure 3A–B). We next used a small molecule inhibitor of POLθ (ART558) to determine if HFK 8E6 cells were reliant on POLθ activity, using an MTT assay to detect viability. There were no significant differences in viability between HFK LXSN and HFK 8E6 cells when grown in media containing a gradient (0–50 μM) of ART558 concentrations (Figure 3C) However, 8E6 impairs apoptosis by degrading BAK, which complicates the interpretation of these data. To more directly define the extent that 8E6 promotes DSB repair by TMEJ, we determined the extent that POLθ inhibition impaired DSB repair. For this analysis, we induced DSBs by growth media containing zeocin (10 µg/mL, 10 min). Cells were then switched to media containing DMSO (solvent control), 1 μM 558, or 5 μM ART 558. DSBs were detected by immunofluorescence microscopy of pH2AX24 hr later. These concentrations of ART558 were analyzed because they caused the largest differences in viability between HFK 8E6 and HFK LXSN cells. ART558 delayed DSB repair in HFK 8E6 cells but not in HFK LXSN cells, indicating that 8E6 promoted repair by TMEJ (Figure 3D–E). However, although statistically significant, the increased persistence of pH2AX foci was a much lower magnitude than what PARP1 inhibition caused. These experiments were also repeated with higher concentrations of ART558, but no differences in the persistence of pH2AX was detected (not shown). We interpret these data as evidence that 8E6 inhibits multiple types of DNA-PKcs-independent end joining, including but not limited to TMEJ. For simplicity, we will only use the term Alt-EJ moving forward.
DNA-PK inhibition forces 8E6 expressing cells to use Alt-EJ more frequently
NHEJ initiation blocks DSB repair by other pathways, including Alt-EJ (Patterson-Fortin and D’Andrea, 2020; Simsek and Jasin, 2010; Wang et al., 2006). 8E6 does not prevent the initiation of NHEJ as autophosphorylated DNA-PKcs foci form readily in HFK 8E6 cells (Hu et al., 2020). Instead, 8E6 blocks the completion of NHEJ. Thus, we hypothesized that preventing HFK 8E6 cells from initiating NHEJ would force them to repair DSBs via Alt-EJ. To test this, we determined the frequency of Alt-EJ in HFK LXSN and HFK 8E6 cells in the presence of a small molecule inhibitor of DNA-PKcs (1 µM of NU7441) to block NHEJ initiation. As expected, DNA-PKcs inhibition increased both imbedded and terminal Alt-EJ in HFK LXSN cells (Figure 4 and Figure 4—figure supplement 1). However, the increase in HFK LXSN cells did not reach the levels of Alt-EJ in mock-treated HFK 8E6 cells. DNA-PKcs inhibition resulted in a further significant increased use of Alt-EJ in HFK 8E6 cells. Together these data support the conclusion that 8E6 promotes the use of Alt-EJ and that the use of the pathway is further enhanced if cells are not allowed to initiate NHEJ.
When HFK 8E6 cells are allowed to initiate NHEJ, many of the pDNA-PKcs repair complexes that form persist for over 24 hr (Hu et al., 2020). HR factors (e.g., RAD51) are then recruited to these unresolved repair foci (Hu et al., 2022a). However, HR is ultimately unable to repair the lesions. We hypothesized that forcing HFK 8E6 cells to use Alt-EJ by preventing initiation of NHEJ would lead to more efficient DSB repair by allowing them to avoid these abortive attempts at DSB repair. NU7441 (1 mM) was used to block NHEJ initiation and DSBs were detected using pH2AX. As we have shown before, DSBs were more persistent in mock-treated 8E6 cells than LXSN cells. Supporting our hypothesis, DNA-PKcs inhibition made DSB repair more efficient (less pH2AX) in HFK 8E6 cells (Figure 5A–B).
We next determined the extent that DNA-PKcs inhibition in HFK 8E6 cells resulted in increased resolution of RAD51 foci. As we have previously reported (Hu et al., 2022a), DNA-PKcs inhibition increased the persistence of RAD51 in HFK LXSN (Figure 5C–D). A similar increase in RAD51 persistence was also seen in mock-treated HFK 8E6 cells. However, DNA-PKcs inhibition increased the rate of RAD51 resolution in HFK 8E6 (Figure 5C–D). To determine if inhibition of a later step in the NHEJ pathway also increased the efficiency of DSB repair in HFK 8E6 cells, we used a small molecule inhibitor of ligase IV (1 µM of SCR7) to block a near terminal step in NHEJ. Ligase IV inhibition delayed DSB repair in both HFK LXSN and HFK 8E6 cells (Figure 5E–F). Thus, restoration of DSB repair in HFK 8E6 specifically requires inhibition of an early NHEJ step, rather than inhibition of a later step in the pathway.
DNA-PKcs inhibition prevents 8E6 from causing RAD51 foci to form during G1
The persistent RAD51 repair complexes that form in HFK 8E6 cells occur during G1 (Hu et al., 2022a). These observations and the ones described in Figure 4 led us to hypothesize that blocking NHEJ initiation via DNA-PKcs inhibition would prevent 8E6 from causing RAD51 foci to form in G1. To test this, we detected RAD51 foci and cyclin E (G1 marker) after the induction of DSB by zeocin. DNA-PKcs inhibition increased the frequency of HFK LXSN in G1 that contained RAD51 foci (Figure 6A–B). Without DNA-PKcs inhibition, 8E6 increased the frequency of cells in G1 that had RAD51 foci. However, DNA-PKcs inhibition prevented 8E6 from promoting the formation of RAD51 repair complexes during G1. As cyclin E is also expressed during early S phase, this examination was repeated using cyclin A as a marker of cells in S/G2 (Figure 6—figure supplement 1). These experiments confirmed our observations using cyclin E to determine cell cycle position, providing further evidence that DNA-PKcs inhibition prevented 8E6 from promoting the formation of RAD51 foci during G1 (Figure 6—figure supplement 1). We used flow cytometry as a final determinant of cell cycle position, using NUCLEAR-ID Red DNA staining to select cells in G1 based on DNA content and then determined the frequency with which these cells stained for RAD51 (Figure 6—figure supplement 2). Consistent with our hypothesis, DNA-PKcs inhibition increased the frequency of RAD51 in G1 in HFK LXSN and prevented 8E6 from promoting RAD51 in G1 (Figure 6C–D).
DNA-PKcs inhibition does not promote HR in cells with 8E6
The data above demonstrate that the attenuation of DSB repair by 8E6 can be overcome by inhibiting DNA-PKcs. DNA-PKcs inhibition also increases the resolution of RAD51, suggesting that DNA-PKcs inhibition may prevent 8E6 from attenuating HR. To test this, we measured HR efficiency using an established HR reporter, described in Figure 7A (Pierce et al., 1999). In vector control U2OS cells (U2OS LXSN), NU7441 increases HR efficiency (Figure 7B–C). This is consistent with the established idea that NHEJ and HR compete for access to DSBs (Bunting et al., 2010; Chapman et al., 2012; Orthwein et al., 2015). As previously reported, 8E6 in U2OS cells (U2OS 8E6) decreased HR efficiency (Wallace et al., 2015). However, DNA-PKcs inhibition did not prevent 8E6 from hindering HR (Figure 7B–C). Consistent with a p300-dependent mechanism (see next section for more details), U2OS cells expressing 8E6 with the residues responsible for binding p300 deleted (U2OS 8E6Δ132–136) were similar to vector control (Figure 7B–C). These cell lines have been previously described (Hu et al., 2020; Wallace et al., 2015). U2OS cells are routinely used to probe 8E6 biology because 8E6 retains its ability to alter DNA repair in these cells (Hu et al., 2020; Wallace et al., 2015).
8E6 promotes Alt-EJ by destabilizing p300
8E6 delays DSB repair by binding/destabilizing p300, leading us to hypothesize that the residues of 8E6 that facilitate the interaction with p300 were important for the phenotypes described thus far (Hu et al., 2022a, Hu et al., 2020; Wallace et al., 2015). To test this, we examined U2OS LXSN, U2OS 8E6, and U2OS 8E6Δ132–136 cells. Consistent with the p300-dependent mechanism, 8E6Δ132–136 shows similar Alt-EJ frequency with U2OS LXSN (Figure 7A–B). The deletion of these residues from 8E6 has been shown to prevent some but not all aspects of 8E6 biology. As a result, we examined p300 knockout N/TERT immortalized HFK cells. P300 knockout led to increases in both terminal and imbedded Alt-EJ (Figure 8C–D). To further confirm p300 dependence, we used a small molecule inhibitor of p300 (1 µM of CCS1477) to block p300 activity. Consistently, CCS1477 increased both terminal and imbedded Alt-EJ (Figure 8E–F).
DNA-PKcs inhibition prevents the formation of RAD51 foci in G1 that is caused by loss of p300
We next confirmed that the ability of 8E6 to allow RAD51 foci to form in G1 was blocked by DNA-PKcs inhibition in U2OS cells (Figure 9A–B). Consistent with a p300-dependent mechanism, RAD51 foci were more likely to occur in G1 when DNA-PKcs was inhibited in U2OS LXSN and U2OS 8E6Δ132–136 cells. Further confirming a p300-dependent mechanism of action, HFK without p300 displayed an increased frequency of RAD51 staining in G1 that could be overcome by DNA-PKcs inhibition (Figure 9C–D). As a final confirmation of the p300 dependence of this phenotype, we treated N/TERT HFKs with 1 µM of CCS1477 to block p300 activity. Inhibition of p300 alone increased the frequency of RAD51 staining in G1 as did DNA-PKcs inhibition alone (Figure 9E–F). However, RAD51 staining in G1 was not increased by their dual application.
Alt-EJ deletions occur at a higher frequency in 8E6 expressing cells
HFK cells were transfected with the 8E6 expression vector or empty vector and equally passaged under selective pressure to establish HFK 8E6 versus control HFK LXSN cell lines. These were subjected to whole genome sequencing to determine the type and frequency of mutations accumulated within each cell type. This yielded 933,163,972 high-quality paired-end sequence reads for LXSN and 624,472,714 reads for 8E6 (Table 1). Normalizing LXSN reads to the same number as 8E6 and alignment to the human genome yielded 92.4% genome coverage for both and slightly more read depth for LXSN (mean 45.4 for LXSN compared to 44.6 for 8E6) (Table 2). Variant calling including single-nucleotide polymorphisms (SNPs) and insertions or deletions (indels) were also more for LXSN cells than 8E6 expressing cells with one variation every 3 kb in LXSN compared to one every 3333 bases in 8E6 (Table 3).
Repair by Alt-EJ can frequently result in indels with characteristic sequence signatures of flanking microhomology. Using previously established algorithms to identify such characteristic indels (McIntosh et al., 2020), we quantified the accumulation of Alt-EJ DNA scars in 8E6 cells compared to LXSN control cells. Counting mutations unique to both LXSN and 8E6 cells, this quantified the frequency of short and long deletions with increasing stretches of microhomology up to 20 bp (Tables 4–6, Tables 7 and 8); small insertions (<5 bp) resulting from templated synthesis in trans previously reported in Drosophila (Tables 4 and 7, Table 9); or large insertions (≥18 bp) (Tables 4 and 10) resulting from templated synthesis in cis (snapback synthesis) described in mouse embryonic stem (ES) cells. Deletions represented the most abundant Alt-EJ mutations accumulating in both cell lines. Despite a larger number of annotated indels unique to LXSN (168,005 for LXSN compared to 104,322 for 8E6), this revealed a significantly higher frequency of accumulated short deletions (2–29 bp) in the 8E6 cell line (86.75% in 8E6 compared to 82.33% in LXSN, p-value = 2.2e-16). The frequency of short Alt-EJ type deletions was most significant for deletions with short (≤10 bp) stretches of flanking microhomology (Table 4). This trend was similar for long deletions (30–500 bp), albeit not statistically significant (Table 5, Table 6). Accumulating insertions both small and large bearing Alt-EJ signatures were far less in frequency in both 8E6 and LXSN cell lines though significantly more in LXSN cells (Tables 7 and 8) consistent with more non-Alt-EJ mutations annotated in LXSN.
Discussion
We have previously shown that 8E6 attenuates the two most prominent DSB pathways (HR and NHEJ) (Hu et al., 2020; Wallace et al., 2015). However, 8E6 delays rather than abrogates DSB repair, leaving the question of how DSBs are repaired in cells expressing 8E6. Here, we show that 8E6 promotes DSB repair via Alt-EJ (Figure 10). Because Alt-EJ can be accomplished using different sets of repair factors, we used small molecule inhibitors to determine which Alt-EJ components were most essential for DSB repair in 8E6 expressing cells. This analysis demonstrated that in HFK 8E6 cells, DSBs become three- to fourfold more persistent when PARP1 is inhibited but only ~1.5-fold more persistent when POLθ is inhibited. Thus, 8E6 promotes DSB repair that is dependent on PARP1 and to a lesser extent on POLθ.
Moreover, we show that HFK 8E6 cells can be further induced to repair DSBs via Alt-EJ by the inhibition of an early step during NHEJ (DNA-PKcs) but not a later step in the pathway (Ligase IV). The increased use of Alt-EJ induced by DNA-PKcs inhibition prevented 8E6 from generating previously described DSB repair defects including the formation of RAD51 foci in G1 and delayed DSB repair. We also show that DNA-PKcs inhibition does not cause an increase in HR. Further, we provide mechanistic insight into these phenomena by showing that 8E6 promotes Alt-EJ via p300 degradation and that DNA-PKcs inhibition can prevent RAD51 foci from forming in G1 because of p300 loss. Finally, we provide whole genome sequence analysis that demonstrates 8E6 expression results in a significantly higher frequency of short deletions bearing microhomology signatures of Alt-EJ. These Alt-EJ deletions appear more frequently as short deletions (2–29 bp) bearing short stretches of microhomology (2–10 bp) (Table 4).
There are some limitations to this study. This includes potential off target activity of inhibitors. Further complicating the analysis, many DNA repair factors (e.g., PARP1) are involved in more than one repair pathways. We also acknowledge that U2OS may contain mutations (such as an ATR mutation) that change the mechanisms of DSB repair pathway choice (Elbakry et al., 2021; Juhász et al., 2018). However, the known ways that 8E6 alters DSB repair are consistent in U2OS and primary (HFK) cells, suggesting that results obtained in U2OS cells can be extrapolated to other cell types (Hu et al., 2020; Wallace et al., 2015). 8E6 was also expressed with the other HPV8 early genes as would be the case during a natural infection. We are currently working to determine the extent that these other early genes augment or suppress 8E6-induced changes in DSB repair.
Many of the findings in this study are consistent with the existing hypotheses about DSB repair and the potential for beta-HPV infections to promote NMSCs (Gheit, 2019; Kremsdorf et al., 1982; Pfister, 2003). For example, our data in vector control (LXSN) cells show that Alt-EJ is increased when NHEJ is inhibited and that DSB repair is not significantly delayed by inhibition of Alt-EJ. This is compatible with the idea that Alt-EJ is primarily used when there are defects in either HR or NHEJ (Deriano and Roth, 2013; Simsek and Jasin, 2010; Truong et al., 2013). Similarly, by showing that 8E6 promotes Alt-EJ, a mutagenic DSB repair pathway, we provide evidence in support of the idea that transient β-HPV infections may promote tumorigenesis by causing mutations (Tommasino, 2019; Viarisio et al., 2018; Wendel and Wallace, 2017). We also show that p300 restricts DSB repair by Alt-EJ and that at least some of the DSB repair defects caused by p300 loss or 8E6 expression can be overcome by inhibition of DNA-PKcs. Does DNA-PKcs inhibition represent a feasible approach to block the increased mutagenesis associated with 8E6 expression? Or on the contrary, does DNA-PKcs inhibition promote more mutations generated by Alt-EJ?
Induction of Alt-EJ is not limited to HPV8 as a recent report demonstrated that a different HPV protein from a different genus of HPV (HPV16 E7) also increases the use of Alt-EJ (Leeman et al., 2019). Further, HPV positive head and neck squamous cell carcinomas and other cancers with downregulated TGF-b signaling have an elevated frequency of mutations with signatures of repair by Alt-EJ (Liu et al., 2018; Liu et al., 2021). We show that 8E6 increases the frequency of short deletions with very short stretches of flanking microhomology. Thus, the ability to promote Alt-EJ seems to have evolved in two separate genes in the HPV family (once in HPV16 E7 and once in HPV8 E6) (Leeman et al., 2019). Why would these viruses both evolve ways to promote Alt-EJ? Perhaps, this can be linked to their ability to impair HR and/or NHEJ (Hu et al., 2020; Wallace et al., 2015). Unrepaired DSBs are highly lethal, thus we infer a strong selective pressure for HPV (a non-lytic virus) to find alternative way(s) to repair DSBs.
The data presented here invoke other interesting thoughts. Do HPV8 infections leave Alt-EJ signatures during natural infections? If so, can Alt-EJ signatures be used to provide evidence that naturally occurring beta-HPVs cause mutations? The ability to identify mutations caused by past transient beta-HPV infections would provide the long-sought after evidence that these infections permanently harm host cells.
Materials and methods
Cell culture and reagents
Request a detailed protocolPrimary HFKs were derived from neonatal human foreskins. Immortalized human foreskin keratinocytes (N/TERT HFK) provided by Michael Underbrink (University of Texas Medical Branch). Both HFKs were grown in EpiLife medium (MEPICF500, Gibco), supplemented with 60 µM calcium chloride (MEPICF500, Gibco), human keratinocyte growth supplement (MEPICF500, Gibco), and 1% penicillin-streptomycin (PSL02-6X100ML, Caisson). Both keratinocyte cell lines were derived from different donors. They were maintained at low passage number (<20 passages) for all experiments described in this manuscript. Their identity was confirmed by morphology and growth in selective media. Only keratinocytes are viable in EpiLife media.
U2OS were maintained in DMEM supplemented with 10% FBS and 1% penicillin-streptomycin. Zeocin (J67140-XF, Alfa Aesar) was used to induce DSBs (10 µg/mL, 10 min). Their identity was confirmed by the presence of the DR-GFP construct used to measure HR. This reporter cassette only exists in these cells within our lab. 8E6 expression (or lack thereof) was confirmed by rtPCR. All cell lines in this study underwent regular mycoplasma testing. NU7441 (S2638, Selleckchem) was used to inhibit DNA-PKcs phosphorylation (1 µM) and verify the pDNA-PKcs antibody. siRNA DNA-PKcs was used to further validate pDNA-PKcs antibody. KU55933 (Sigma-Aldrich, SML1109) was used to validate RAD51 antibody as previously described (Bakr et al., 2015). CCS1477 (CT-CCS1477, Chemietek) was used to inhibit p300 activity (1 µM). ART558 (HY-141520, MedChem Express) was used to inhibit Pol Theta activity. Alt-EJ plasmids (#113619, #113620, #113625, #113626, Addgene) were used to measure Alt-EJ efficiency.
Immunofluorescence microscopy
Request a detailed protocolCells were seeded onto either 96-well glass-bottom plates and grown overnight. Cells treated with zeocin (10 μg/mL, 10 min) were fixed with 4% paraformaldehyde. Then, 0.1% Triton-X was used to permeabilize the cells, followed by blocking with 3% bovine serum albumin. Cells were then incubated with the following antibodies: RAD51 (ab1837, Abcam, 1:200), cyclin E (4132S, Cell Signaling), pH2AX S139 (9718S, Cell Signaling), and cyclin A (ab39, Abcam). The cells were washed and stained with the appropriate secondary antibodies: Alexa Fluor 594 (red) goat anti-rabbit (A11012, Thermo Fisher Scientific), Alexa Fluor 488 (green) goat anti-mouse (A11001, Thermo Fisher Scientific). After washing, the cells were stained with 10 µM DAPI in PBS and visualized with the Zeiss LSM 770 microscope. Images were analyzed using the ImageJ techniques previously described (Murthy et al., 2018). Cyclin E intensity was measured for each cell. Average cyclin E intensity of cells grown in media without growth factor for 4 hr was used to define the threshold of cyclin E positive.
Flow cytometry
Request a detailed protocolCells were collected from 6 cm plates, at about 80–90% confluence, by using trypsinization. Cells were washed with cold PBS and fixed with 95% cold ethanol for 10 min at –20°C. Cells were stained with anti-RAD51 antibody (ab1837, Abcam, 1:100) and Alexa Fluor 488 goat anti-mouse (A11001, Thermo Fisher Scientific). After washing, cells were resuspended in 200 µL PBS and NUCLEAR-ID Red DNA stain (ENZ-52406, Enzo Life Science), and incubated in the dark at room temperature for 30 min. Samples were analyzed by a BD Accuri C6 Plus Flow Cytometer.
Immunoblotting
Request a detailed protocolAfter being washed with ice-cold PBS, cells were lysed with RIPA Lysis Buffer (VWRVN653-100ML, VWR Life Science), supplemented with Phosphatase Inhibitor Cocktail 2 (P5726-1ML, Sigma) and Protease Inhibitor Cocktail (B14001, Bimake). The Pierce BCA Protein Assay Kit (89167-794, Thermo Fisher Scientific) was used to determine protein concentration. Equal protein lysates were run on Novex 3–8% Tris-acetate 15 Well Mini Gels (EA03785BOX, Invitrogen) and transferred to Immobilon-P membranes (IPVH00010, Thermo Fisher Scientific). Membranes were then probed with the following primary antibodies: GAPDH (sc-47724, Santa Cruz Biotechnologies, 1:1000) P300 (sc-48343, Santa Cruz Biotechnologies). CAS9 (65832S, Cell Signaling). Pol Theta (PA5-69577, Thermo Fisher Scientific). After exposure to the matching HRP-conjugated secondary antibody, cells were visualized using SuperSignal West Femto Maximum Sensitivity Substrate (34095, Thermo Fisher Scientific).
Transfection and Alt-EJ assay
Request a detailed protocolHFK cells were plated in 3 mL of complete growth medium in a 6 cm plate. Cells were used at 80% confluency. Two µg of plasmids were diluted in 200 µL Xfect transfection reagent (631317, Takara). The mixture was incubated at room temperature for 15 min. The transfection mixture was added to each well drop-wise and incubated for 48 hr at 37°C. Cells were harvested for flow cytometry analysis.
Whole genome sequencing analysis
Request a detailed protocolHFK cells were transfected with the 8E6 expression vector or empty vector and equally passaged under selective pressure to establish HFK 8E6 versus control HFK LXSN cells. DNA was extracted using Trizol (Invitrogen), libraries were prepared using the Illumina DNA PCR-Free prep kit and sequenced with paired 300-bp v1.5 reads on a Novaseq 6000. All protocols were performed using the standard procedures provided by the manufacturers. Sequences have been deposited in the NCBI SRA database with accession number (PRJNA 856469).
Variant calling and Alt-EJ mutational analysis
Request a detailed protocolWhole genome sequences for passage matched LXSN and 8E6 cell lines were analyzed for indel detection using the best practices for variant discovery analysis outlined by the Broad Institute (DePristo et al., 2011). Briefly, reads from the sequenced samples were preprocessed for quality control by removing adapter sequences and low-quality reads (<Q20). Quality paired-end reads from LXSN were normalized to the same number of reads from 8E6 (624,472,714) and both sets were mapped to the human reference genome (GRCh37/hg19) using Burrows-Wheeler Aligner (BWA v0.7.17) (Li and Durbin, 2010). Aligned reads were sorted and duplicate reads were removed using Picard-tools (v2.25.5). Raw genomic variants, including SNPs and DNA insertions and deletions (indels), were identified using the GATK HaplotypeCaller algorithm in GATK (v4.1.9.0) (‘Genomics in the Cloud [Book]’, n.d.). The identified indels were extracted and filtered for quality with GATK VariantFiltration (QUAL score normalized by allele depth: QD <2.0, Phred-scaled probability: FS >200, and Symmetric odds ratio test: SOR >10), followed by base quality score recalibration using GATK BaseRecalibrator and applied to the aligned bam files. A second round of variant calling (GATK Haplotype Caller) was performed using recalibrated (analysis-ready) SNPs and indels further filtered for quality. SnpEff 5.0 was used to annotate the genetic variants.
Alt-EJ mutational signatures were counted as previously described (McIntosh et al., 2020). Briefly, indels previously identified by the 1000 Genomes Project or shared between LXSN and 8E6 cell lines were removed. Alt-EJ mutational signatures bearing flanking microhomology were then determined in R for small deletions (2–29 bp) and long deletions (30–500 bp), small insertions (<5 bp), and large snap-back insertions (≥18 bp). For analysis of deletions with microhomology, the sequence adjacent to each deletion was extracted based on the human genome (R package BSgenome.Hsapiens. UCSC.hg19) and annotated according to the number of contiguous matching nucleotides between the deletion and the adjacent sequence, counted from the beginning of each sequence (microhomology). Different microhomology length thresholds (ranging from 2 to 20 bp) were then used to filter the results into Alt-EJ matching deletions unique to each LXSN and E6 cell line.
To detect synthesis-dependent small insertions, we extracted ±15 nt surrounding the insertion (from the hg19 human genome) and identified whether there was a repeated sequence either 5’ or 3’ to the insertion that matched the full insertion and also at least 2 nt on either side. Matches also needed to have a gap of at least 1 nt between the two sequences.
For synthesis-dependent snapback (large) insertions, only the insertion sequence was considered. Each sequence was compared to its own reverse-complement sequence and any repeats of ≥7 nt also separated by at least 4 nt were labeled as snapback-like. To avoid spurious false negatives, only insertions of at least (2×7)+4 = 18 nt were considered.
Cell viability assay
Request a detailed protocolTen thousand cells/well were seeded on a 96-well plate and grown for 24 hr. Zeocin treatment (10 μg/mL, 10 min) was applied 24 hr after seeding. Then ART558 dose series was added and incubated for 48 hr at 37°C. Forty-eight hr after treatment, 10 µL/well of MTT solution (10 mg/mL) was added for 24 hr. Subsequently, wells were incubated with 100 µL solubilization solution for 24 hr and the optical density measured at 640 nm. Doses for ARTT558 : 0, 1, 5, 10, 15, 20, 25, 50.
Statistical analysis
Request a detailed protocolAll values are represented as mean ± standard error (SE). Statistical differences between groups were measured by using Student’s t-test. p-Values in all experiments were considered significant at less than 0.05. Statistical significance of Alt-EJ mutational frequencies between LSXN and 8E6 in each instance (small or large deletions or insertions) were determined by Fisher’s exact t-test using prop.test in R.
Data availability
Sequences have been deposited in the NCBI SRA database with accession number (PRJNA 856469).
-
NCBI BioProjectID PRJNA856469. Beta human papillomavirus 8E6 promotes alternative end-joining.
References
-
Involvement of ATM in homologous recombination after end resection and Rad51 nucleofilament formationNucleic Acids Research 43:3154–3166.https://doi.org/10.1093/nar/gkv160
-
ConferenceHuman papillomavirus infection and skin cancer risk in organ transplant recipientsThe Journal of Investigative Dermatology. Symposium Proceedings. pp. 207–211.https://doi.org/10.1046/j.0022-202x.2001.00048.x
-
Human papillomavirus and posttransplantation cutaneous squamous cell carcinoma: a multicenter, prospective cohort studyAmerican Journal of Transplantation 18:1220–1230.https://doi.org/10.1111/ajt.14537
-
Brca1 modulates the autophosphorylation status of DNA-PKcs in S phase of the cell cycleNucleic Acids Research 42:11487–11501.https://doi.org/10.1093/nar/gku824
-
High beta-HPV DNA loads and strong seroreactivity are present in epidermodysplasia verruciformisThe Journal of Investigative Dermatology 129:1026–1034.https://doi.org/10.1038/jid.2008.317
-
CREBBP and p300 lysine acetyl transferases in the DNA damage responseCellular and Molecular Life Sciences 75:1325–1338.https://doi.org/10.1007/s00018-017-2717-4
-
Cbp/P300 in cell growth, transformation, and developmentGenes & Development 14:1553–1577.https://doi.org/10.1101/gad.14.13.1553
-
Alternative end-joining repair pathways are the ultimate backup for abrogated classical non-homologous end-joining and homologous recombination repair: implications for the formation of chromosome translocationsMutation Research. Genetic Toxicology and Environmental Mutagenesis 793:166–175.https://doi.org/10.1016/j.mrgentox.2015.07.001
-
Characterizing DNA repair processes at transient and long-lasting double-strand DNA breaks by immunofluorescence microscopyJournal of Visualized Experiments 1:57653.https://doi.org/10.3791/57653
-
Human papillomavirus load in eyebrow hair follicles and risk of cutaneous squamous cell carcinomaCancer Epidemiology, Biomarkers & Prevention 22:719–727.https://doi.org/10.1158/1055-9965.EPI-12-0917-T
-
Xrcc3 promotes homology-directed repair of DNA damage in mammalian cellsGenes & Development 13:2633–2638.https://doi.org/10.1101/gad.13.20.2633
-
Dna double-stranded breaks induce histone H2AX phosphorylation on serine 139The Journal of Biological Chemistry 273:5858–5868.https://doi.org/10.1074/jbc.273.10.5858
-
Factors determining DNA double-strand break repair pathway choice in G2 phaseThe EMBO Journal 30:1079–1092.https://doi.org/10.1038/emboj.2011.27
-
Beta human papillomavirus and associated diseasesActa Cytologica 63:100–108.https://doi.org/10.1159/000492659
-
Alternative end-joining is suppressed by the canonical NHEJ component XRCC4-ligase IV during chromosomal translocation formationNature Structural & Molecular Biology 17:410–416.https://doi.org/10.1038/nsmb.1773
-
Hpv and skin carcinogenesisPapillomavirus Research 7:129–131.https://doi.org/10.1016/j.pvr.2019.04.003
-
Rnf8 has both Ku-dependent and independent roles in chromosomal break repairNucleic Acids Research 48:6032–6052.https://doi.org/10.1093/nar/gkaa380
-
Dna fragmentation caused by an overdose of ZeocinJournal of Bioscience and Bioengineering 116:644–646.https://doi.org/10.1016/j.jbiosc.2013.05.004
-
Parp-1 and Ku compete for repair of DNA double strand breaks by distinct NHEJ pathwaysNucleic Acids Research 34:6170–6182.https://doi.org/10.1093/nar/gkl840
-
Dissection of DNA double-strand-break repair using novel single-molecule forcepsNature Structural & Molecular Biology 25:482–487.https://doi.org/10.1038/s41594-018-0065-1
-
Loss of genome fidelity: beta hpvs and the DNA damage responseFrontiers in Microbiology 8:2250.https://doi.org/10.3389/fmicb.2017.02250
Article and author information
Author details
Funding
National Institute of General Medical Sciences (P20GM130448)
- Nicholas Wallace
National Institutes of Health (NCI R15 CA242057 01A1)
- Nicholas Wallace
U.S. Department of Defense (CMDRP PRCRP CA160224 (NW))
- Nicholas Wallace
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We appreciate KSU-CVM Confocal Core for our immunofluorescence microscopy, Michael Underbrink for providing the TERT-immortalized HFKs, and Jeremy M Stark for providing plasmids of the alt-EJ reporter assay. This research was supported by the National Institute of General Medical Sciences (NIGMS) of the National Institutes of Health (NIH, P20GM130448); NIH Research Enhancement Award (NCI R15 CA242057 01A1); US Department of Defense (CMDRP PRCRP CA160224 [NW]); and Johnson Cancer Research Center of Kansas State University. We acknowledge the University of Kansas Medical Center Genomics Core for your core support services in publications containing Genomics Core generated data. We would also like to note the following NIH supporting grants – Kansas Intellectual and Developmental Disabilities Research Center (NIH U54 HD 090216), the Molecular Regulation of Cell Development and Differentiation – COBRE (P30 GM122731-03) – the NIH S10 High-End Instrumentation Grant (NIH S10OD021743) and the Frontiers CTSA grant (UL1TR002366) at the University of Kansas Medical Center, Kansas City, KS 66160. We would also like to acknowledge Dr Massimo Tommasino. His kindness and influence will be missed.
Copyright
© 2023, Hu 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
-
- 1,110
- views
-
- 137
- downloads
-
- 5
- 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
-
- Cancer Biology
- Cell Biology
PML, a multifunctional protein, is crucial for forming PML-nuclear bodies involved in stress responses. Under specific conditions, PML associates with nucleolar caps formed after RNA polymerase I (RNAPI) inhibition, leading to PML-nucleolar associations (PNAs). This study investigates PNAs-inducing stimuli by exposing cells to various genotoxic stresses. We found that the most potent inducers of PNAs introduced topological stress and inhibited RNAPI. Doxorubicin, the most effective compound, induced double-strand breaks (DSBs) in the rDNA locus. PNAs co-localized with damaged rDNA, segregating it from active nucleoli. Cleaving the rDNA locus with I-PpoI confirmed rDNA damage as a genuine stimulus for PNAs. Inhibition of ATM, ATR kinases, and RAD51 reduced I-PpoI-induced PNAs, highlighting the importance of ATM/ATR-dependent nucleolar cap formation and homologous recombination (HR) in their triggering. I-PpoI-induced PNAs co-localized with rDNA DSBs positive for RPA32-pS33 but deficient in RAD51, indicating resected DNA unable to complete HR repair. Our findings suggest that PNAs form in response to persistent rDNA damage within the nucleolar cap, highlighting the interplay between PML/PNAs and rDNA alterations due to topological stress, RNAPI inhibition, and rDNA DSBs destined for HR. Cells with persistent PNAs undergo senescence, suggesting PNAs help avoid rDNA instability, with implications for tumorigenesis and aging.
-
- Cancer Biology
- Computational and Systems Biology
Assay for Transposase-Accessible Chromatin sequencing (ATAC-Seq) is a widely used technique to explore gene regulatory mechanisms. For most ATAC-Seq data from healthy and diseased tissues such as tumors, chromatin accessibility measurement represents a mixed signal from multiple cell types. In this work, we derive reliable chromatin accessibility marker peaks and reference profiles for most non-malignant cell types frequently observed in the microenvironment of human tumors. We then integrate these data into the EPIC deconvolution framework (Racle et al., 2017) to quantify cell-type heterogeneity in bulk ATAC-Seq data. Our EPIC-ATAC tool accurately predicts non-malignant and malignant cell fractions in tumor samples. When applied to a human breast cancer cohort, EPIC-ATAC accurately infers the immune contexture of the main breast cancer subtypes.