1. Introduction
Radiation therapy allows cancer cells elimination through the formation of DNA damage. DNA double-stranded breaks (DSBs) can be repaired by homologous recombination (HR) [
1] or non-homologous end joining (NHEJ) [
2]. Single-stranded DNA (ssDNA) lesions are subject to repair by nucleotide excision repair [
3] and base excision repair [
4], while DNA mismatch repair functions to correct replication errors in newly synthesized DNA [
5]. During the S phase of the cell cycle unrepaired SSBs stalls the progression of DNA replication and causes fork reversal [
6]. Translesion DNA synthesis (TLS) emerged as a mechanism of DNA damage tolerance (DDT) that enables the completion of DNA replication leaving the unrepaired damage for future elimination [
7].
During TLS, replicative DNA polymerase is temporarily replaced by a specialized TLS polymerase that has the ability to replicate across DNA lesions. TLS polymerases can use the damaged template for synthesis and bypass DNA lesions, may incorporate the non-complementary nucleotides thereby facilitating replication fork progression [
8]. These polymerases include B-family Pol ζ and Y-family Pol η, Pol κ, Pol ι and Rev1.
Pol ι efficiently bypasses a variety of DNA lesions including lesions induced by oxidative stress but possesses very low accuracy of nucleotide incorporation [
9,
10]. Pol ι efficiently incorporates deoxynucleotides opposite the 3' T of [6-4] PPs, abasic sites [
11], and opposite
N2-adducted guanines [
11], products of cytosine oxidation [
12] and 8-oxoguanine (8-oxoG) [
13], the most abundant and mutagenic oxidative lesion in DNA. Pol ι also plays TLS-independent functions in DDT.
Several SNPs in
POLI (encoding Polι) have been linked to the development of melanoma, prostate cancer, lung adenocarcinoma and squamous cell carcinoma [
14,
15,
16].
POLI expression was elevated in breast cancer [
13], bladder cancer [
17] and oesophageal squamous cell carcinoma (ESCC) [
18]. At the same time, Pol ι inhibited tumorigenesis and mutagenesis in some cancers, despite its extreme replication infidelity. Polι deficiency in mice leads to increased susceptibility to UV-induced skin tumors [
19,
20]. Pol ι also protects from chemically induced lung cancer. Polι was found to be the
Par2 (pulmonary adenoma resistance 2) gene responsible for adenoma and adenocarcinoma susceptibility in mice [
21,
22,
23]. In addition, the activation of Polι triggers the Erk and JNK signaling pathways, which play a pivotal role in promoting invasion and metastasis via EGFR-ERK-mediated epithelial to mesenchymal transition [
24,
25]. Ultimately, this leads to the increase in proliferation and a grim prognosis for individuals with ESCC [
18,
25,
26]. However, whether Polι plays a role in proliferation and ionizing radiation-induced DNA damage response (DDR) of NSCLC cells remains unclear.
Recently, a novel primase-polymerase called PrimPol has been identified in eukaryotic cells [
27]. PrimPol promotes DNA repriming downstream of a lesion or replication block to facilitate fork restart following ultraviolet (UV) damage and is important for inter-strand cross-link traverse [
28]. PrimPol synthesizes DNA primers
de novo, restarting replication in the nucleus and mitochondria after damaged DNA regions or secondary structures—G quadruplexes and R loops [
28,
29,
30,
31,
32,
33,
34,
35]. PrimPol also exhibits properties of a translesion DNA polymerase, effectively incorporating nucleotides opposite a number of small DNA lesions, including 8-oxoguanine, O
6-methylguanine, 5-formyluracil, and apurine-apyrimidine sites
in vitro [
36,
37,
38]. PrimPol also quite effectively incorporates complementary dCMPs opposite the intra-strand cisplatin crosslink
in vitro [
39]. Deletion of the
Polι gene in a mouse embryonic fibroblast cell line slows down the replication of nuclear and mitochondrial DNA under physiological conditions and causes sensitivity to a number of genotoxic agents: UV radiation, methyl methanesulfonate, hydroxyurea, benz[a]pyrene and cisplatin [
29,
32,
34,
40,
41,
42,
43].
PRIMPOL −/− cells were more sensitive to UV-C irradiation in colony survival assays than wild-type and Pol η-deficient cells [
40]. Mutations in
PRIMPOL have also been found in human cancers [
44,
45,
46,
47,
48,
49] and compiled in Genomic Data Commons Data Portal, GDC [
50] and Catalogue of Somatic Mutations in Cancer, COSMIC databases [
51]. The essential role of
PRIMPOL in ensuring the smooth progression of DNA replication forks highlights the tremendous potential of inhibiting this non-essential mammalian enzyme. This compelling strategy offers a promising avenue for effectively treating various conditions by disrupting DNA synthesis.
The lack of a particular type of DNA polymerase is anticipated to have a profound effect on the DDR. This is expected to increase the replication stress caused by DNA damage [
52,
53], leading to targeting differentiated and cancer stem cells/stem cell-like cells (CSCs) [
54] by potential activation of apoptosis. However, the impact of PRIMPOL and POLI elimination in NSCLC cells exposed to ionizing radiation has not been elucidated.
In this study, we aimed to evaluate the impact of PRIMPOL and POLI elimination on human A549 NSCLC cell line exposed to X-ray irradiation stress. We observed significant differences in cell survival, apoptosis, DNA repair and the expression of cancer stem cells (CSCs) markers in PRIMPOL and POLI knockout cells compared to wild-type cells.
3. Discussion
PrimPol plays a crucial role in facilitating the resumption of DNA replication after encountering DNA lesions, such as DNA interstrand cross-links [
28]. This remarkable ability is demonstrated through its involvement in DNA replication-restart pathways [
41,
63,
64,
65]. Pol iota (Polι) seems to be primarily responsible for regulating a DDT pathway choice that determines how DNA damage is tolerated [
66]. This pathway choice involves a competitive interaction between PrimPol-mediated repriming and fork reversal. It should be noted that this function is distinct from Pol iota's role in DNA replication and DNA damage repair [
9]. Instead, it serves as an important protective enzyme, even though it does play a significant role in TLS. Together, these data suggest that the significance of specialized DNA polymerases in the DDR extends beyond their participation in TLS. TLS is generally regarded as mechanistically simple and straightforward, but error-prone DDT pathway as it has high potential for mutagenesis [
67]. These include dormant origin firing, homologous recombination (HR) and specialized DNA polymerases involved in TLS bypass of replication-stalling lesions [
68,
69]. Heterogeneous fidelity of TLS polymerases together with different combinations of polymerases/lesions/sequence contexts can result in distinct sets of mutated genes and thus the development of specific cancers [
70]. Recent studies indicated that TLS primase/polymerase PrimPol is possibly associated with the development and chemosensitivity of many cancers [
32,
71]. In present study, we aimed to evaluate the impact of PRIMPOL and POLI genes on metastatic potential and radiosensitivity of human NSCLC cells.
As was previously shown by others, the ability to migrate, structure and mechanics of cancer cells are linked directly to their metastatic potential [
72,
73,
74,
75,
76]. Recently, machine learning algorithm-based single-cell analysis revealed possible participation of
PRIMPOL in various crucial biological functions such as cell invasion, cell cycle regulation, apoptosis, and hypoxia response [
56].
Recent evidence has shown a compelling connection between the ability of cancer cells to migrate through 8 µm pores and their capacity to indent and display metastatic behavior [
77]. Hence, we employed this method to assess the motility and invasiveness of both wild type and polymerase knock-out cells in response to a serum gradient while confined within a porous matrix of Boyden chambers. Our data has uncovered a significant decline in a vital biophysical (mechanobiological) characteristic of cancer cells - their ability to withstand migration while being confined to a specific area for three consecutive days. This decline was observed in both knockout (KO) cell lines, as shown in
Figure 1, when compared to their syngeneic wild type counterparts. It is true that both KO cell lines had a decrease in S-phase cells after 72 hours of cultivation (Fig. 2b). However, it is highly improbable that cell cycle disruptions are the cause of the observed reduced migration.
Figure 3 (a-c) illustrates the remarkable increase in baseline apoptosis in KO cells. The levels reach a peak nearly twice as high as those in their parental counterparts. On the seventh day of cell cultivation, the frequency of basal apoptosis increases significantly in both knockout cell lines, as shown in
Figure 3d. This finding can be pinpointed as the possible cause behind the observed migratory decline effect. Thus, our current study offers experimental evidence affirming the AI-based algorithm's prediction of a positive correlation between
PRIMPOL expression and migration/invasiveness, as well as cell cycle progression [
56]. Additionally, we have demonstrated an inverse relationship between
PRIMPOL expression and the level of apoptosis in NSCLC cells.
Furthermore, we have effectively presented similar discoveries and connections regarding the role of
POLI deletions in these essential biological functions of the same NSCLC cells. It was previously demonstrated that Polι promotes cell migration and invasion of ESCC and breast cancer cells [
18,
24,
25,
78]. Our data strongly supports the effectiveness of small interfering RNA (siRNA) targeting the
POLI gene in inhibiting invasion of the A549 lung cancer cell line [
79]. It is important to mention that this treatment has had contrasting effects on apoptosis as compared to the results we obtained from knocking out
POLI gene (Fig. 3d) of the same cell line. The use of siRNA to partially silence this gene had little impact on the cellular MTT-detected activity, which was mistakenly attributed to apoptosis by these authors. Our groundbreaking study reveals that completely knocking out
POLI gene significantly enhances the occurrence of true apoptosis, establishing a clear inverse correlation with the activity of the
POLI gene.
Although
PRIMPOL and
POLI knockouts augment basal apoptosis, their effect on reproductive cell survival following exposure to ionizing radiation (IR) was increased (Fig. 4). Notably, this increase reached statistical significance only in
PRIMPOL-/- cells after irradiation with 6 Gy, indicating their increased DDT. While our current study did not prioritize it, the data we collected do not rule out the possibility of anastasis. Anastasis is a remarkable cellular recovery phenomenon that saves dying cells from the edge of death [
80]. The unique response to IR stress, characterized by activated metastasis, epithelial mesenchymal transition (EMT), and DNA damage repair mechanisms promoting anastasis [
81], may still occur differentially in
PRIMPOL and
POLI knockouts. Thus, our data is in line with the clinical observation that the expression of
PRIMPOL is decreased in lung adenocarcinoma when compared to normal tissues. This decrease was significantly linked to patients' prognosis [
56].
The absence of specialized DNA polymerases can have a significant impact on various aspects of DDR. It is highly likely that without these enzymes, there would be an increase in replication stress caused by DNA damage [
52,
53]. To evaluate whether the efficiency of DNA DSBs repair involved into increased DDT of
PRIMPOL -/- and
POLI -/- cell clones, we evaluated γH2AX and pATM foci kinetics.
ATM, ATR, and DNA-PK are three phosphoinositide 3-kinase (PI3K)-related kinases that phosphorylate H2AX in response to DNA damage [
82]. While ATM is considered as a major physiological mediator of H2AX phosphorylation in response to DSB formation [
83], ATR phosphorylates H2AX in response to single-stranded DNA breaks and during replication stress, such as replication fork arrest [
84]. We observed no discrepancies in the number of pATM foci between the cell lines being compared. This compelling finding suggests that the absence of
PRIMPOL and
POLI does not have any impact on the cellular capacity to repair 'true' DNA DSBs.
The levels of γH2AX are a crucial indicator of DNA damage throughout the entire organism. Examining these levels can provide invaluable insights into the degree of replication stress in the absence of specialized DNA polymerases. Therefore, we examined the accumulation of γH2AX in cancer cells lacking either
PRIMPOL or
POLI and exposed to ionizing radiation (IR), a commonly used stimulus that induces replication stress. We observed a significant disparity in the number of γH2AX foci between the compared cells (Fig.5a). We observed a decrease in γH2AX levels in cells that lacked
POLI, whereas the absence of
PRIMPOL led to a significant increase in both the baseline and residual γH2AX foci numbers induced by ionizing radiation (Fig. 5a). Of note,
PRRIMPOL knock-out also augments the Rad51 foci, albeit during the first four hours after irradiation (Fig. 5c). Upon irradiation PrimPol is able to bind to replication protein A (RPA) [
85]. RPA is the single-stranded DNA (ssDNA)-binding protein and one of the proteins required for the maintenance of genetic information [
86]. Since the pATM foci number remained unchanged among populations (Fig. 5b), the remarkable increase in the number of γH2AX foci in
PRIMPOL-/- cells could be linked to the known deficiencies in fork restart after replication arrest [
85].
POLI expression in Esophageal Squamous Cell Carcinoma (ESCC) was negatively correlated with overall survival of ESCC patients with postoperative adjuvant radiotherapy, suggesting its role in ESCC radioresistance [
87]. Recently, the requirement of
POLI for optimal γH2AX and 53BP1 accumulation has been shown [
66]. The study conducted by Sabrina F. Mansilla et al. yielded similar results, demonstrating a reduction in γH2AX and 53BP1 levels following treatment of Polι-depleted cells with cisplatin [
66]. Of note,
POLI knock-out also mitigated the accumulation of pan-nuclear Rad51 foci during the first 4 and 24 h after irradiation (Fig. 5c). Our data indicates that
POLI deletion might be directly associated with either translational or post-translational down-regulation of RAD51 as was shown previously [
66,
87]. Important, RAD51 is one of three well-known proteins (behind BRCA1, BRCA2) that negatively regulate the
PRIMPOL expression. The role of Rad51 in homologous recombination is well-established. It is worth mentioning that this protein is also linked to fork reversal (FR) [
88] and is a favored binding partner of POLι in non-differentiated cancer stem-like cells [
54]. Consequently, the concurrent reduction of Rad51 and elimination of POLι are highly likely to stimulate a greater reliance on PrimPol-dependent repriming [
89]. This process might significantly impair the p53-POLι-mediated FR-DDT pathway in cell populations that resemble CSCs. If that's true, then it supports the idea that POLI has a critical influence on the selection of the differential DDT pathway [
54]. Such a change in the DDT pathway could potentially explain why
POLI -/- cells have a lower clonogenic survival rate than
PRIMPOL -/- cells following exposure to IR stress in our study.
In the absence of
POLI,
PRIMPOL to take over and resulted in enhanced repriming, DNA replication acceleration, and bypassing of the checkpoint activation during the S phase. Consequently, this led to a surge in chromosome instability during the M phase [
66]. Our findings (Fig. 2c,d) support the observation of G2/M arrest instead of S phase arrest after exposure to IR. According to our data,
POLI knocked-out cells likely have a remarkable ability to amplify replication stress after being exposed to IR. This finding challenges the conventional understanding of other specialized DNA polymerases, which typically counter this stress instead.
Cancer stem-like cells (CSCs) are a subpopulation of tumor cells that can drive tumor initiation and recurrence [
90].
PRIMPOL mutations have been discovered in human cancers, but what's even more fascinating is its over-expression in CD 133-positive CSC-like glioblastoma subtypes [
91]. These points to potential functions in combating replication stress in diseased tissues. Since PrimPol seems to have crucial functions in sustaining RF progression, despite not being a necessary enzyme in mammals, the development of PrimPol inhibitors to interrupt DNA synthesis could potentially offer a viable approach for treating various conditions. Here, we present evidence that a dosage of 6 Gy IR could trigger a remarkable surge in the proportion of CD133+ and CD133+/CD44+ populations in
PRIMPOL-/- cells. The data indicate that their greater clonogenic survival (Fig.4) could be linked to the existence of these CSC-like cells. Additional research is necessary to understand better the mechanisms that determine the choice of the DDT pathway and the subsequent induction of CSC-like cells by
PRIMPOL-/- cells in response to IR stress. It is critical to exercise caution when developing inhibitors for specialized DNA polymerases in cancer treatment, as our data shows a strong indication for potential complications.
4. Materials and Methods
4.1. Cell Cultures and Cultivation conditions
The PRIMPOL and POLI knocked out A549 cells were kindly provided by Dr. A. Makarova (Institute of Gene Biology of Russian Academy of Sciences, Moscow, Russia) and Prof. D.O. Zharkov (ICBFM SB RAS, Novosibirsk, Russia). Cells were obtained using the CRISPRrispr/Cas9 system and the pSpCas9(BB)-2A-GFP vector as previously described [
55]. The absence of PrimPol and Polι in cells was verified by immunoblotting (
Figure S1) Cells were cultured in DMEM with high glucose content GlutaMAX™ Supplement (Gibco, USA), 10% fetal bovine serum, 100 U/mL penicillin, and 100 U/mL streptomycin at 37°C and 5% CO
2.
4.2. Trans-Well Migration Assay Using Boyden Chambers
Briefly, cells were serum starved overnight and 104 cells were seeded in 0,1 ml of serum-free cell-type specific media in the upper compartment of transwell inserts (EDM Milli-pore, Billerica, MA) with 8 µm pores. Serum-supplemented growth media was placed in the lower chamber, serving as a chemoattractant. Cells were incubated in the wells for 72 h at 37 °C and 5% CO2. The cells from the upper and lower compartments were collected, by incubation with 0,2 and 0,3 ml of trypsin-EDTA solution (0.25%; Biological Industries, Cromwell, CT), respectively, for 5 min at 37 °C and under 5% CO2. After trypsin neutralization with serum-containing growth media the cells collected from upper and lower chamber were counted using hemocytometer. For each cell type 6–9 independent experiments were performed, and the percentage of cancer cells from all cell lines that cross through 8 µm Boyden chamber membranes was determined.
4.3. Apoptosis Analysis
To quantify the proportion of early-stage apoptotic cells the commercial kit “Vybrant Apoptosis Assay Kit #4” with YO-PRO-1 and PI for Flow Cytometry (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA; catalog number: V13243) was used. The cells were stained according to the supplemented manufacturer protocol. Immediately (0 hours) or 48 hours post-radiation, cells were collected and washed in cold phosphate-buffered saline (PBS). A total of 1 μL of YO-PRO-1 stock solution and 1 μL of PI stock solution were added to each 1 mL containing 1 × 106 cells. Cells were incubated on ice for 20–30 min and analyzed by flow cytometry (BD FACSCalibur, Becton Dickinson, San Jose, CA, USA) using 488 nm excitation with green fluorescence emission for YO-PRO-1 (i.e., 530/30 bandpass) and red fluorescence emission for PI (i.e., 610/20 bandpass). A total of 50,000 events were acquired for each sample and analyzed with FlowJo™ Software (Becton Dickinson, San Jose, CA, USA).
4.4. Colony Formation and Soft Agar Assay
Cells were exposed to 0 Gy, 2 Gy, 4 Gy, and 6 Gy of X-ray irradiation, plated on 60 mm Petri dishes in 150, 500, 1000, and 2000 cells/well, respectively and incubated at 37°C 5% CO2. After two weeks, cells were fixed with 100% methanol for 15 min at room temperature, followed by Giemsa staining for 15 min. Only colonies containing more than 50 cells were counted. Plating efficiency (PE) and survival fractions (SF) were calculated using the following equations:
4.5. γH2AX, pATM and Rad51 Foci Analysis
Cells were seeded in 96-well plates, fixed 1–24 h after irradiation in 4% paraformaldehyde for 15 min at room temperature and permeabilized in 0.3% Triton-X 100 (in PBS, pH 7.4) supplemented with 2% bovine serum albumin (BSA). Cells were incubated for 1 h at room temperature with primary mouse monoclonal antibody against γH2AX (dilution 1:200, clone JBW 301, Cat. # 05-636, Merck-Millipore, Burlington, VT, USA), primary mouse monoclonal antibody against pATM (dilution 1:200, clone 10H11.E12, Merck Millipore, Burlington, VT, USA) and with primary rabbit polyclonal antibody against Rad51 (dilution 1:200, Merck Millipore, Burlington, VT, USA) diluted in PBS with 1% BSA and 0.3% Triton-X 100. After several rinses with PBS, cells were incubated for 1 h with secondary antibodies IgG (H + L) goat anti-mouse (Alexa Fluor 555 conjugated, dilution 1:800; Cat. # A-21424, Merck-Millipore, Burlington, VT, USA), goat anti-rabbit (Alexa Fluor 488 conjugated, dilution 1:500; Cat. # A-11008, Merck Millipore, Burlington, VT, USA) diluted in PBS (pH 7.4) with 1% BSA. Cells were imaged using Nikon Eclipse Ni-U microscope (Nikon, Tokyo, Japan) equipped with a high-definition camera ProgResMFcool (Jenoptik AG, Jena, Germany). Filter sets used were UV-2E/C (340–380 nm excitation and 435–485 nm emission), B-2E/C (465–495 nm excitation and 515–555 nm emission), and Y-2E/C (540–580 nm excitation and 600–660 nm emission). A total of 300–400 cells were imaged for each data point. Foci were counted by manual scoring.
4.6. Analysis of the Proportion of Cancer Stem Cells by Flow Cytometry
Cells were collected by trypsinization, washed in ice-cold PBS (pH = 7.4) and 1 × 106 cells per sample were incubated with Anti-CD133 Antibody, Alexa Fluor® 488 conjugated (MAB4310X, Sigma-Aldrich, Darmstadt, Germany) and with monoclonal Anti-CD44−PE antibody (SAB4700187, Sigma-Aldrich, Darmstadt, Germany) for 30 min at 4 °C. Cells were analyzed by flow cytometry (BD FACSCalibur, Becton Dickinson, San Jose, CA, USA). A total of 50,000 events were acquired for each sample and the proportion of positive cells was analyzed with BD CellQuest Pro 5.1 software (Becton Dickinson, San Jose, CA, USA).
4.7. Cell Cycle Analysis by Flow Cytometry
Exponentially growing cells were irradiated (6 Gy) and incubated for 48 and 72 hours post-radiation. Cells were then harvested and 1 × 106 cells per sample were resuspended in ice-cold PBS. Cells were fixed in ice-cold 70% ethanol for 30 min at 4 °C and then kept at −20 °C until analyzed. Before analysis, specimens were washed twice with PBS, and resuspended in 0.5 mg/ml PI containing 50 μl of 100 μg/ml RNase for 30 min in the dark. Cells were analyzed by flow cytometry (BD FACSCalibur, Becton Dickinson, San Jose, CA, USA). A total of 50,000 events were acquired for each sample and the percentage of cells in the different phases of the cell cycle was analyzed with FlowJo™ Software (Becton Dickinson, San Jose, CA, USA).
4.8. Statistics
Statistics were performed using GraphPad Prism 9.0.2.161 (GraphPad Software, San Diego, CA, USA) software. Statistical significance was tested using the Student t-test. The results are represented as means ± SEM of more than three independent experiments. Significance levels were denoted by asterisks: * p < 0.05, ** p < 0.01, *** p < 0.001, ****p < 0.0001.