DNA double-strand breaks in human induced pluripotent stem cell reprogramming and long-term in vitro culturing
© The Author(s). 2017
Received: 21 September 2016
Accepted: 24 February 2017
Published: 21 March 2017
Human induced pluripotent stem cells (hiPSCs) play roles in both disease modelling and regenerative medicine. It is critical that the genomic integrity of the cells remains intact and that the DNA repair systems are fully functional. In this article, we focused on the detection of DNA double-strand breaks (DSBs) by phosphorylated histone H2AX (known as γH2AX) and p53-binding protein 1 (53BP1) in three distinct lines of hiPSCs, their source cells, and one line of human embryonic stem cells (hESCs).
We measured spontaneously occurring DSBs throughout the process of fibroblast reprogramming and during long-term in vitro culturing. To assess the variations in the functionality of the DNA repair system among the samples, the number of DSBs induced by γ-irradiation and the decrease over time was analysed. The foci number was detected by fluorescence microscopy separately for the G1 and S/G2 cell cycle phases.
We demonstrated that fibroblasts contained a low number of non-replication-related DSBs, while this number increased after reprogramming into hiPSCs and then decreased again after long-term in vitro passaging. The artificial induction of DSBs revealed that the repair mechanisms function well in the source cells and hiPSCs at low passages, but fail to recognize a substantial proportion of DSBs at high passages.
Our observations suggest that cellular reprogramming increases the DSB number but that the repair mechanism functions well. However, after prolonged in vitro culturing of hiPSCs, the repair capacity decreases.
KeywordsHuman induced pluripotent stem cells DNA double-strand breaks γH2AX 53BP1 Long-term in vitro culture DNA repair
Human induced pluripotent stem cells (hiPSCs) hold great promise for clinical applications because of their potential to differentiate into all three embryonic germ layers [1–3]. To use hiPSCs in cell therapy or disease modelling , it is fundamental that they possess an intact genome. Much research has been performed in the field of genome maintenance in mouse and human embryonic stem cells (hESCs). However, less is known about the causes of genomic aberrations and the functionality of repair mechanisms in hiPSCs . In general, the genomic instabilities in hiPSCs may be introduced: 1) by pre-existing mutations in source cells; 2) during reprogramming; and 3) during in vitro expansion of the hiPSCs. It has been reported that at least 50% of the single-nucleotide variations in hiPSCs pre-existed in the source cells . The process of reprogramming itself represents a serious risk of mutation acquisition. Primarily, deletions of tumour suppressor genes were observed during reprogramming . Using episomal vectors may lower the risk of reprogramming-associated genome changes . The culturing of pluripotent stem cells (PSCs) in vitro is probably the main cause of the accumulation of genomic instabilities, resulting from the adaptation to culture conditions and clonal selection during passaging. The data reported by Taapken et al.  support this idea, indicating that the types and frequency of karyotypic abnormalities are similar between hESCs and hiPSCs. In contrast, the results of Laurent et al.  revealed slight differences in the distribution of subchromosomal variations between hESCs and hiPSCs. Interestingly, in their study , prolonged in vitro culturing was associated with oncogene duplication.
One of the key techniques for monitoring DNA integrity is the detection of DNA double-strand breaks (DSBs). DSBs are a severe type of DNA damage that may cause irreversible changes in the genomic content of the cell. They are induced by internal factors such as the by-products of cell metabolism or replication stress, or by external factors such as exposure to irradiation or chemical agents . A damaged cell may arrest the cell cycle until the lesions are repaired. If the DNA damage is not successfully repaired, apoptosis is commonly induced to prevent the propagation of chromosomal aberrations. The repair of DSBs is executed either by fast non-homologous end-joining (NHEJ) or more precise homologous recombination (HR). Both mechanisms contribute to DSB repair in a cell cycle-specific manner. NHEJ occurs at all phases of the cell cycle but is primarily responsible for DSB repair in the G1 stage. HR occurs predominantly in the late S and G2 phases . Published data suggest that pluripotent cells exert stronger genomic protection and can repair DNA lesions more efficiently than differentiated somatic cells [5, 12–14]. However, a strong DNA protective mechanism may cause the pluripotent cells to be more prone to apoptosis.
Various DNA damage-response proteins have been used as markers of DSBs, including phosphorylated histone H2AX (known as γH2AX) and p53-binding protein 1 (53BP1) . The generation of γH2AX foci at the site of DNA lesions precedes the formation of 53BP1 foci [16–18]. Several studies have confirmed that 53BP1 functions exclusively in NHEJ and that it inhibits the 5′ end resection needed for HR [19–21]. In contrast, γH2AX influences both NHEJ and HR . The foci formation of γH2AX is dependent on the cell cycle phase [22–24]. S/G2 phase cells exhibit more γH2AX foci than do cells in G1 phase because of replication-related DSBs. Cell-cycle dependency has not been observed for 53BP1 .
In the present study, we compared the genomic integrity of fibroblasts and pluripotent stem cells. We used fluorescence microscopy to visualize the DSBs recognized by γH2AX and 53BP1 in three hiPSC lines and one hESC line at low or high passage numbers and in one line of source cell fibroblasts. Each hiPSC line is unique and represents a different reprogramming approach, as described in the Methods. We also aimed to detect differences in the ability to recognize DSBs artificially induced by γ-irradiation and their decrease over time. The measurements were conducted with respect to cell cycle stage, and the data were analysed separately for the G1 and G2/S phases. Thus, we aimed to elucidate genomic stability during hiPSC generation and in vitro culturing.
hiPSC generation and cell culture
Human dermal fibroblasts (hDFs; kindly provided by the National Tissue Centre Inc., Brno, Czech Republic) and CD34+ haematopoietic progenitors (blood sample kindly provided by the Department of Internal Medicine, Haematology and Oncology, Masaryk University, and University Hospital Brno, Czech Republic) were used as source cells for the generation of hiPSCs as described in Šimara et al., 2014 [26, 27]. For this study, we used the hiPSC line CBIA-3–CD34+ haematopoietic progenitors reprogrammed by the Sendai virus (CytoTuneTM-iPS Reprogramming Kit; Thermo Fisher Scientific, Waltham, MA, USA), hiPSC line CBIA-5–fibroblasts reprogrammed by the Sendai virus, and hiPSC line CBIA-7–fibroblasts reprogrammed by episomal vectors (Epi5™ Episomal iPSC Reprogramming Kit; Thermo Fisher Scientific). The CCTL-14 hESC line  was a kind gift from the Department of Histology and Embryology (Faculty of Medicine, Masaryk University, Brno, Czech Republic). All three hiPSC lines and the hESC line were maintained in the form of colonies on irradiated mouse embryonic fibroblast feeder cells (MEFs; 2.5 × 105 cells per 3.5-cm dish) in DMEM/F12 (1:1) supplemented with 20% knock-out serum replacement, 2 mM l-glutamine, 100 μM non-essential amino acids, 1% penicillin/streptomycin, 0.1 mM 2-mercaptoethanol, and 10 ng/ml basic fibroblast growth factor (bFGF) (all from Thermo Fisher Scientific). The medium was changed daily. Markers of pluripotency (Oct-3/4, Sox2, Nanog, and SSEA4) were detected as described previously , and all three hiPSC lines were positive for all of these markers (Additional file 1: Figure S1). A teratoma formation assay confirmed that the hiPSCs could differentiate into all three germ layers (Additional file 2: Figure S2).
The following cell passage numbers (p) were used (low and high): CBIA-3 at p27 and p76, CBIA-5 at p19 and p65, CBIA-7 at p25 and p67, and CCTL-14 at p30 and p302. hDFs were used at p6. No CD34+ blood progenitors used to generate CBIA-3 were available for DSB analysis.
Prior to irradiation, the hiPSCs and hESCs were feeder depleted by culturing on a Geltrex® matrix for 3 days. Essential 8™ medium (Thermo Fisher Scientific) was changed daily. The cells were then irradiated by ionizing radiation (IR; 0.5 Gy/min; 137Cs; 1 and 4 Gy) to induce DSBs and fixed in 4% paraformaldehyde at 0.5, 2, and 6 h after irradiation.
The dose of 1 Gy was selected for the experiments based on published results [12, 29] and our DSB count measurement after 1 Gy or 4 Gy irradiation (data not shown). The peak value of the foci number was recorded 0.5 h after irradiation; therefore, this time point was selected for the study of the functionality of the repair system.
Immunocytochemical staining was used to visualize the DSBs and distinguish between the cell cycle stages G1 and S/G2. Four hours before fixation, a nucleoside analogue of thymidine, 5-ethynyl-2′-deoxyuridine (EdU; Thermo Fisher Scientific), was added at a final concentration of 10 μM. The cells were fixed in 4% paraformaldehyde and permeabilized in 0.2% Triton-X (both from Sigma-Aldrich, St. Louis, MO, USA). Overnight incubation with primary antibodies against γH2AX (Biolegend, San Diego, CA, USA) and 53BP1 (Santa Cruz Biotechnology, Dallas, TX, USA) was followed by 1 h of incubation with a secondary antibody conjugated with Alexa 555 (Cell Signaling Technology, Danvers, MA, USA). The samples were stained with the Click-iT® EdU Alexa Fluor® 488 Imaging Kit (Thermo Fisher Scientific) to visualize EdU according to the manufacturer’s instructions. Finally, the nuclei were stained with Hoechst dye (BisBenzemide H33342; 1 μg/ml; Sigma-Aldrich).
Fluorescence microscopy and image analysis
Fluorescent signals were detected using the Zeiss Axiovert 200 M system (Carl Zeiss, Oberkochen, Germany). The images were captured using a CoolSNAP HQ2 CCD camera in the wide-field mode (Photometrics, Tucson, AZ) at –30 °C. Thirty 3-μm slices were acquired in each field at a resolution of 1392 × 1040 pixels. The pixel size of the images was 124 × 124 nm. Between 500 and 1000 cells were analysed from each microscopic slide. Two slides, γH2AX and 53BP1, were prepared from each sample and each time point.
Acquiarium software, developed by our group, was used to acquire and analyse the images . Acquiarium is open source software available for download at our group’s official website (http://cbia.fi.muni.cz/projects/acquiarium.html). During the analysis, individual cells in the field of view were first cropped manually. Next, the nucleus of each cell, stained with Hoechst dye, was recognized automatically using the “Find objects (hysteresis thresholding)” plugin. We used the Gaussian filter in the preprocessing phase (with sigma = 1), the threshold was calculated using the two-level Otsu method, and we defined the minimum size of an object to exclude the parts of adjacent cells. This plugin defined the area in which we counted γH2AX or 53BP1 foci. For this purpose, we employed the eMax algorithm described in  using the parameters sigma = 1, a spot height threshold of 80, and a maximum spot size of 800, which we set empirically. The EdU signal was quantified based on the total intensity calculated in the nucleus. The threshold for the separation of EdU-negative (G1) and EdU-positive (S/G2) cells was computed in MATLAB (Mathworks) using the Otsu method.
Flow cytometry analysis
To assess early apoptosis, cells were stained with Annexin-V fluorescein isothiocyanate (FITC) and 7-amino-actinomycin D (7-AAD; BD Via-Probe) in Annexin-V binding buffer (Miltenyi Biotec, Bergisch Gladbach, Germany). From each sample, approximately 1 × 105 cells were processed. All of the samples were measured using a BD FACS Canto II flow cytometer (Becton-Dickinson). BD FACSDiva (Becton-Dickinson) software was used for the data analysis.
For each time point, approximately 1 × 106 cells were lysed in RIPA buffer. The total protein concentration was assessed using a Pierce™ BCA Protein Assay Kit (Thermo Fisher Scientific). Laemmli buffer was added, and the samples were separated by SDS-PAGE. The proteins were transferred onto polyvinylidene fluoride membranes, and the membranes were blocked with 5% milk in TBS-Tween for 1 h. The membranes were then incubated with a 1:1000 dilution of PARP and GAPDH primary antibodies (both from Cell Signaling Technology) in TBS-Tween with 5% milk at 4 °C overnight. Subsequently, the membranes were incubated with the secondary antibody (1:5000 anti-rabbit HRP; Cell Signaling Technology) for 1 h at room temperature, and the blots were developed using the Clarity™ Western ECL Substrate (Bio-Rad Laboratories, Hercules, CA, USA) according to the manufacturer’s instructions.
Comparison of two data sets was performed using Student’s t test. Multi-group assays were analysed by a one-way analysis of variance (ANOVA) in conjunction with Tukey’s test. A level of P ˂ 0.05 was considered to be statistically significant.
Discrimination between the cell cycle phases using EdU increases the accuracy of analysing DNA lesions
Taken together, this series of experiments demonstrates the robustness of the approach that we have developed to visually discriminate between G1 and S/G2 cells in situ. Our data show that, using this technique, we can identify changes in cell cycle progression. In the context of cell cycle-associated differences in numbers of γH2AX and 53BP1 foci, this approach is extremely useful and was employed for all the following analyses in this study. The Acquiarium software also represents an extremely valuable tool for complex and automated microscope image analysis.
Reprogramming is accompanied by increased numbers of γH2AX and 53BP1 foci, but this trend is reversed with prolonged in vitro culturing
Number of foci per cell
Non-irradiated control cells
Irradiated cells (1Gy)
G1 (EdU-negative cells)
S/G2 (EdU-positive cells)
G1 (EdU-negative cells)
S/G2 (EdU-positive cells)
No. of foci ± SEM
No. of cells
No. of foci ± SEM
No. of cells
No. of foci ± SEM
No. of cells
No. of foci ± SEM
No. of cells
1.13 ± 0.11
29.9 ± 1.32
19.89 ± 0.36
45.49 ± 0.88
1.49 ± 0.07
19.64 ± 1.23
15.55 ± 0.22
19.20 ± 0.91
5.56 ± 0.22
22.98 ± 0.38
21.87 ± 0.34
28.71 ± 0.34
2.13 ± 0.12
8.76 ± 0.35
19.01 ± 0.31
23.1 ± 0.42
4.42 ± 0.20
17.91 ± 0.41
8.01 ± 0.27
7.97 ± 0.35
1.46 ± 0.07
4.73 ± 0.18
14.45 ± 0.29
13.66 ± 0.31
5.69 ± 0.23
25.1 ± 0.47
28.6 ± 0.52
41.00 ± 0.55
4.03 ± 0.18
15.3 ± 0.49
17.84 ± 0.26
23.32 ± 0.30
2.98 ± 0.13
10.98 ± 0.26
11.57 ± 0.27
24.45 ± 0.39
1.56 ± 0.07
4.71 ± 0.21
8.04 ± 0.17
10.84 ± 0.20
5.89 ± 0.25
23.77 ± 0.45
26.48 ± 0.44
36.71 ± 0.52
2.73 ± 0.13
10.42 ± 0.31
20.05 ± 0.35
21.94 ± 0.38
2.56 ± 0.13
11.78 ± 0.25
12.41 ± 0.26
12.73 ± 0.30
1.51 ± 0.07
6.5 ± 0.27
16.9 ± 0.32
17.62 ± 0.34
4.45 ± 0.20
12.83 ± 0.31
20.99 ± 0.34
30.37 ± 0.57
2.73 ± 0.14
10.28 ± 0.36
16.77 ± 0.27
19.13 ± 0.39
6.66 ± 0.29
19.30 ± 0.49
16.26 ± 0.28
19.38 ± 0.34
4.75 ± 0.21
14.34 ± 0.52
21.49 ± 0.30
24.20 ± 0.42
hiPSCs lose their DNA repair capacity after prolonged maintenance in vitro
The above experiments demonstrated that, under normal culture conditions, hDFs, hiPSCs, and hESCs all have characteristic numbers of γH2AX and 53BP1 foci. However, based on these measurements we cannot resolve whether this is due to differences in the level of “spontaneous” DNA damage, DNA repair capability (recognition of DNA lesions), or both. It is understood that the amount of DSBs in cultured cells caused by γ-irradiation is about the same for the same dose of irradiation, regardless of the type of cell. With this holding true, the numbers of γH2AX and 53BP1 foci detected in cells irradiated by the same dose of γ-rays should then reflect the capability of the DNA repair machinery to recognize DSBs rather than the level of DNA damage. In the following series of experiments, we built on this presumption to study the DNA repair efficiency of human pluripotent stem cells. We irradiated the respective cells (hDFs, hiPSCs, and hESCs) with the same dose of γ-rays (1 Gy) and then determined the number of γH2AX and 53BP1 foci at three different time points after irradiation (0.5, 2, and 6 h). It has previously been shown that the levels of γH2AX and 53BP1 loaded onto chromatin usually reach a maximum at approximately 15–30 min after ionizing irradiation [32–34]. Based on this data, we used 30 min as the starting point. Two additional time points (2 and 6 h) then provided information on how DNA repair is sustained.
Figure 3c and d show the numbers of γH2AX and 53BP1 foci 30 min after γ-irradiation in EdU-negative cells. As expected for normal cells, hDFs exhibited a dramatic increase to 19.9 and 15.6, respectively. This represents an 18-fold (for γH2AX) and 10-fold (for 53BP1) increase over their numbers in non-irradiated controls, which indeed mirrors a massive initiation of DNA repair pathways (Fig. 4b). Surprisingly, although the numbers of both types of foci were higher in non-irradiated hiPSCs (irrespective of their passage number) than in hDFs (see the previous section), this was not the case for irradiated hiPSCs. Specifically, at 30 min after irradiation, low-passage hiPSCs produced 21.9 to 28.6 γH2AX foci and 17.8 to 20.1 53BP1 foci, thus always exceeding the corresponding numbers observed in hDFs. In contrast, high-passage hiPSCs produced only 8.0 to 12.4 γH2AX foci and 8.0 to 16.9 53BP1 foci. In other words, in hiPSCs, their prolonged passaging dramatically diminished the numbers γH2AX and 53BP1 foci induced by γ-rays to levels below or similar to those observed in hDFs.
As described in the previous section, the numbers of “spontaneous” γH2AX and 53BP1 foci were, for all cell types and categories studied here, always higher in S/G2 (EdU-positive) then in G1 (EdU-negative) cells. In hiPSCs, the fold-increase ranged from three-times in high-passage CBIA-5 cells (53BP1) to 4.6-times in high-passage CBIA-7 cells (γH2AX), and the changes were consistently statistically significant (Fig. 4a). This overall trend was also retained in γ-irradiated cells (at 30 min after irradiation); however, the actual fold-increase (S/G2 versus G1) was much lower, in four cases showing either no changes or statistically insignificant changes (for both γH2AX and 53BP1 foci in high-passage CBIA-3 and CBIA-7 hiPSCs) (Fig. 4a). Specifically, for hiPSCs, the fold-increase ranged from none to 2.1-times (24.6/11.6) for γH2AX foci in high-passage CBIA-5 cells. The percent increase of foci (both γH2AX and 53BP1) after treatment with 1 Gy was higher in the cells in G1 phase than in those in S/G2 phase (Fig. 4b). Taken together, this set of data reveals that a high level of spontaneous DNA damage (replicative stress occurring in S/G2 phase) dramatically distorts the outcome of γ-irradiation as measured by the numbers of γH2AX and 53BP1 foci.
We also analysed hESCs in parallel to hiPSCs to determine whether the studied phenomena are associated with de-differentiation to pluripotency or with the pluripotency per se. Overall, the differences between irradiated low- and high-passage hESCs were much less pronounced than those in hiPSCs. This conclusion is substantiated by the numbers of: 1) γH2AX foci in low- and high-passage EdU-negative hESCs (21.0 versus 16.3); 2) 53BP1 foci in low- and high-passage EdU-negative hESCs (16.8 versus 21.5) (Fig. 3c and d); 3) γH2AX foci in low- and high-passage EdU-positive hESCs (30.4 versus 19.4); and 4) 53BP1 foci in low- and high-passage EdU-positive hESCs (19.1 versus 24.2) (Fig. 4a). It is of note that in high-passage hESCs (both EdU-negative and EdU-positive) the numbers of 53BP1 foci (but not of γH2AX foci) even increased compared to those typical for low-passage hESCs. Finally, the steepness of the time-course curves indicated that the decrease was more similar to hDFs than to hiPSCs (Fig. 5c and d).
A DNA molecule is unstable and subject to internal and external harmful factors. Correct functioning of the DNA repair mechanisms is, therefore, essential for the maintenance of genomic integrity. In the field of hiPSC research, only cells with an intact genome can be used for clinical application. Unfortunately, the generation and expansion of hiPSCs in vitro causes genomic instability. In our research, we monitored the amount of DSBs, either spontaneous or irradiation-induced, in three lines of hiPSCs (CBIA-3, CBIA-5, and CBIA-7) at low or high passage numbers, as well as in original source cell fibroblasts (hDFs). One hESC line (CCTL-14) was also examined. Our goal was to shed light on the reaction of the cells to reprogramming and on the prolonged in vitro culturing of pluripotent cells. Here, we focused specifically on the kinetics of DSB generation and repair, cell cycle speed changes, triggering of apoptosis, and cell viability. Special attention was paid to the cell cycle phase of individual cells.
We selected two markers of DSBs—the phosphorylated histone variant γH2AX and its binding partner, the DNA repair mediator 53BP1. Fluorescence microscopy was chosen to detect these proteins because it offers two main advantages over other methods such as Western blotting. First, the expression of 53BP1 does not change; only its localization at DNA damage sites is affected. Second, analysis at the single-cell level assures a higher sensitivity and allows for the discrimination between cells at various cell cycle phases. We employed EdU staining, which discriminates between the G1 and S/G2 phases of the cell cycle. By incorporating a nucleoside analogue of thymidine (EdU) into the DNA during replication, only cells in the S or G2 stage are labelled positive . The images were analysed using Acquiarium software. This software allows us to reliably determine the foci number together with the intensity of the EdU signal for each individual cell and to analyse the data from hundreds of cells per sample on a large scale. Our method for separating EdU-negative (G1 phase) cells from EdU-positive (S/G2 phase) cells is based on plotting the EdU intensity levels in histograms and using the Otsu method to find the threshold. Using this method, we revealed a longer cell cycle in somatic cells compared to pluripotent cells, which is in accordance with previously published data [36–39] and justifies the use of this method for cell cycle discrimination on the single-cell level. This approach also assures consistency among samples.
While counting the numbers of γH2AX and 53BP1 foci, it is of utmost importance to know exactly which phase of the cell cycle each individual cell is in at the moment. Our data show that cells in S/G2 phase contain more γH2AX and 53BP1 foci than cells in G1 phase and that this difference is more pronounced in non-irradiated controls. These foci emerge due to replication stress during S phase [12, 22, 23, 40]. The replication-related foci play a critical role in the comparison of DSB numbers, especially between different cell types. As long as the cells have a similar cell-cycle length (e.g. pluripotent cells at a similar passage number), the number of DSBs could be compared relatively well without using cell cycle discrimination. However, the following factors influence the cell cycle speed and should be considered: 1) pluripotent cells have been reported to have a shorter cell cycle than differentiated somatic cells [36–39]; 2) pluripotent cells at high passages may have an increased rate of proliferation [39, 41, 42]; and 3) irradiation induces cell cycle arrest through checkpoints [43–45]. We analysed the foci number separately for the EdU-negative and EdU-positive groups to eliminate the replication stress bias. Our data indicate a higher percent increase of foci upon γ-irradiation of cells in G1 phase, which are not burdened by replication-related foci. The cell-cycle dependency was confirmed for both γH2AX and 53BP1 markers. In general, fewer foci were detected for 53BP1 than for γH2AX, suggesting that 53BP1 is a less sensitive DSB marker with a lower capacity to recognize DSBs than γH2AX. It is known that the HR pathway plays a pivotal role during hiPSC generation , and 53BP1 promotes the NHEJ repair pathway while inhibiting the HR pathway [19–21]. In contrast, γH2AX influences both the NHEJ and HR pathways, and 53BP1 does not bind to all of the γH2AX foci [10, 11, 18].
Similar research was performed by Suchánková et al. , who measured the formation of γH2AX- and 53BP1-positive nuclear bodies in relation to cell cycle stages. They used genetically modified HeLa-Fucci cells, which are able to express RFP-Cdt1 in the G1 phase and GFP-geminin in the S/G2/M phases to discriminate among the cell cycle phases. They observed a higher number of γH2AX-positive repair foci in the G2 phase than in the G1 phase for both non-irradiated and γ-irradiated (5 Gy) HeLa cells. In contrast to our work, they did not observe such a difference for 53BP1. It is of note that different cell types as well as a different radiation dose (1 Gy) were used in our study compared to Suchankova et al., and it has been previously published that foci formation upon ionizing radiation may vary between cell types and is radiation-dose dependent [32, 34, 47, 48].
In our study, we worked with three unique hiPSC lines that were derived from two independent cell types (dermal fibroblasts and blood cells) and reprogrammed either by the Sendai virus or episomal vectors. This selection of samples enables us to generalize our conclusions for hiPSCs to a certain extent. To avoid bias caused by replication-related foci, we further analysed γH2AX and 53BP1 foci numbers only in the G1 (EdU-negative) subgroup. Our results indicate that spontaneously occurring DSBs are best recognized by both markers in hiPSCs at low passage, while fewer foci were observed in hiPSCs at high passage and in source fibroblasts. A low foci number in fibroblasts, therefore, is increased significantly after reprogramming into hiPSCs (either by Sendai virus or episomal vectors) and then decreases again after long-term in vitro passaging. Our results are consistent with recently reported data showing that H2AX plays a critical role in iPSC generation. Gonzáles et al. reported an increase in γH2AX during the cellular reprogramming of mouse embryonic fibroblasts independent of viral integration . The HR pathway was confirmed to be essential for the error-free repair of DSBs in both genome-integrating and non-integrating reprogramming. The importance of H2AX at the early stage of reprogramming was also suggested by Wu et al. . Our observations markedly resemble the results of copy number variation (CNV) measurements by Hussein et al. . They concluded that most de novo-formed CNVs are present in early-passage hiPSCs, while fewer CNVs are found in late-passage hiPSCs and fibroblasts. There is a strong connection between CNVs and DSBs because deletions in subtelomeric regions have been shown to be highly sensitive to DSBs and are the major cause of chromosomal instability [51, 52]. Similar results were published by Laurent et al.  who reported a higher frequency of CNVs in pluripotent samples than in non-pluripotent samples and noticed that some of the deletions receded from the population over long-term passaging. Taken together, their data suggest that genomic instability is highest in low-passage hiPSCs, and CNVs vanish during multiple clonal-based passages because most of the mutations do not provide any advantage. However, certain growth-advantageous mutations—for example, defects in genes controlling the cell cycle—may be fixed in the population .
The abovementioned findings imply that more DSBs at low passages are detected as a consequence of reprogramming stress and disappear as the hiPSCs are adapted to in vitro conditions and clonally selected. However, the irradiation experiments revealed that the high-passage hiPSCs cannot recognize DSBs as effectively as hDF source cells, particularly by γH2AX. The lack of ability to recognize the irradiation-induced DSBs was also obvious in all three high-passage hiPSCs lines in the time-course study. These data suggest that hiPSCs lose their repair capacity over multiple passages in vitro. Similar results were published by Zhang et al. on one mouse iPSC line . They confirmed the compromised DNA damage repair capacity of iPSCs compared with the respective source cells after γ-irradiation treatment but did not focus on the length of the in vitro culturing of iPSCs. For potential clinical applications, the length of in vitro culturing time should be reduced to as short as possible. However, a certain amount of time in vitro is unavoidable because of the reprogramming process itself, cell expansion, and clearance of the remaining reprogramming factors (viral particles or vectors).
Of note, low- and high-passage hiPSCs displayed similar apoptotic responses upon γ-irradiation. PARP cleavage peaked 2 h after irradiation, which led to an increase in early apoptosis after an additional 4 h in most of the hiPSC lines. These data suggest that, despite differences in DSB recognition, both low- and high-passage hiPSCs exert DNA protection mechanisms that trigger apoptosis in reaction to γ-irradiation. Increased apoptosis was not observed in somatic hDFs or in the hESC line CCTL-14, suggesting their lower sensitivity to DNA damage.
This study addresses the question of hiPSC capability to repair their DNA using three independent lines of hiPSCs. It shows for the first time that: 1) reprogramming to pluripotency increases the number of DNA double-strand breaks (DSBs) as recognized by the γH2AX and 53BP1 proteins; 2) these DSBs are not due to replicative stress to DNA; and 3) their numbers become reduced during prolonged propagation after reprogramming. It also shows that prolonged passaging of hiPSCs is associated with a decrease in their DNA repair capacity and that this is not the case for the hESC line CCTL-14. From a technical point of view, this study documents that solid accuracy in analysing numbers of DSBs requires discrimination between the cells in G1 and S/G2 phases of their cell cycle. Collectively, hiPSCs at low passage contain more DSBs than hiPSCs at high passage, but they can repair them more efficiently.
p53-Binding protein 1
One-way analysis of variance
Copy number variation
Human dermal fibroblast
Human embryonic stem cell
Human induced pluripotent stem cell
Mouse embryonic fibroblast
Pluripotent stem cell
The authors are very grateful to colleagues from the Department of Histology and Embryology, Faculty of Medicine, Masaryk University, Brno, Czech Republic, for assistance during the teratoma assay. We thank our colleagues from the National Tissue Centre, Brno, for the generous gift of fibroblast cells. We are grateful to Z. Racil from the Department of Internal Medicine, Haematology and Oncology, Masaryk University and University Hospital, Brno, for providing blood samples. We thank American Journal Experts for their editorial services.
This study was generously supported by the Grant Agency of the Czech Republic (302/12/G157) and by the Czech Health Research Council (16-31501A).
Availability of data and materials
All datasets supporting this article have been uploaded as part of the Additional files.
PS was responsible for conception and design, collection and assembly of data, data analysis and interpretation, and manuscript writing. LT was responsible for conception and design, data analysis and interpretation, and manuscript writing. DR was responsible for collection and assembly of data, and data analysis and interpretation. PM was responsible for collection and assembly of data, data analysis and interpretation, and manuscript writing. SS was responsible for collection and assembly of data, data analysis and interpretation, and manuscript writing. AH was responsible for conception and design, data analysis and interpretation, and manuscript writing. IK was responsible for conception and design, data analysis and interpretation, and manuscript writing. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Human tissue samples were collected under approved guidelines set by the Faculty of Medicine, Masaryk University, Brno, and by the National Tissue Centre, Brno. All patients signed an informed consent form.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Takahashi K, Yamanaka S. Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. Cell. 2006;126:663–76.View ArticlePubMedGoogle Scholar
- Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131:861–72.View ArticlePubMedGoogle Scholar
- Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, Nie J, Jonsdottir GA, Ruotti V, Stewart R, et al. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318:1917–20.View ArticlePubMedGoogle Scholar
- Simara P, Motl JA, Kaufman DS. Pluripotent stem cells and gene therapy. Transl Res. 2013;161:284–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Weissbein U, Benvenisty N, Ben-David U. Quality control: genome maintenance in pluripotent stem cells. J Cell Biol. 2014;204:153–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Gore A, Li Z, Fung HL, Young JE, Agarwal S, Antosiewicz-Bourget J, Canto I, Giorgetti A, Israel MA, Kiskinis E, et al. Somatic coding mutations in human induced pluripotent stem cells. Nature. 2011;471:63–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Laurent LC, Ulitsky I, Slavin I, Tran H, Schork A, Morey R, Lynch C, Harness JV, Lee S, Barrero MJ, et al. Dynamic changes in the copy number of pluripotency and cell proliferation genes in human ESCs and iPSCs during reprogramming and time in culture. Cell Stem Cell. 2011;8:106–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Cheng L, Hansen NF, Zhao L, Du Y, Zou C, Donovan FX, Chou BK, Zhou G, Li S, Dowey SN, et al. Low incidence of DNA sequence variation in human induced pluripotent stem cells generated by nonintegrating plasmid expression. Cell Stem Cell. 2012;10:337–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Taapken SM, Nisler BS, Newton MA, Sampsell-Barron TL, Leonhard KA, McIntire EM, Montgomery KD. Karotypic abnormalities in human induced pluripotent stem cells and embryonic stem cells. Nat Biotechnol. 2011;29:313–4.View ArticlePubMedGoogle Scholar
- Shrivastav M, De Haro LP, Nickoloff JA. Regulation of DNA double-strand break repair pathway choice. Cell Res. 2008;18:134–47.View ArticlePubMedGoogle Scholar
- Rothkamm K, Krüger I, Thompson LH, Löbrich M. Pathways of DNA double-strand break repair during the mammalian cell cycle. Mol Cell Biol. 2003;23:5706–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams BR, Golding SE, Rao RR, Valerie K. Dynamic dependence on ATR and ATM for double-strand break repair in human embryonic stem cells and neural descendants. PLoS One. 2010;5:e10001.View ArticlePubMedPubMed CentralGoogle Scholar
- Tichy ED, Pillai R, Deng L, Liang L, Tischfield J, Schwemberger SJ, Babcock GF, Stambrook PJ. Mouse embryonic stem cells, but not somatic cells, predominantly use homologous recombination to repair double-strand DNA breaks. Stem Cells Dev. 2010;19:1699–711.View ArticlePubMedPubMed CentralGoogle Scholar
- Maynard S, Swistowska AM, Lee JW, Liu Y, Liu ST, Da Cruz AB, Rao M, de Souza-Pinto NC, Zeng X, Bohr VA. Human embryonic stem cells have enhanced repair of multiple forms of DNA damage. Stem Cells. 2008;26:2266–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Rothkamm K, Barnard S, Moquet J, Ellender M, Rana Z, Burdak-Rothkamm S. DNA damage foci: meaning and significance. Environ Mol Mutagen. 2015;56:491–504.View ArticlePubMedGoogle Scholar
- Panier S, Boulton SJ. Double-strand break repair: 53BP1 comes into focus. Nat Rev Mol Cell Biol. 2014;15:7–18.View ArticlePubMedGoogle Scholar
- Rogakou EP, Pilch DR, Orr AH, Ivanova VS, Bonner WM. DNA double-stranded breaks induce histone H2AX phosphorylation on serine 139. J Biol Chem. 1998;273:5858–68.View ArticlePubMedGoogle Scholar
- Kleiner RE, Verma P, Molloy KR, Chait BT, Kapoor TM. Chemical proteomics reveals a γH2AX-53BP1 interaction in the DNA damage response. Nat Chem Biol. 2015;11:807–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Bunting SF, Callén E, Wong N, Chen HT, Polato F, Gunn A, Bothmer A, Feldhahn N, Fernandez-Capetillo O, Cao L, et al. 53BP1 inhibits homologous recombination in Brca1-deficient cells by blocking resection of DNA breaks. Cell. 2010;141:243–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Dimitrova N, Chen YC, Spector DL, de Lange T. 53BP1 promotes non-homologous end joining of telomeres by increasing chromatin mobility. Nature. 2008;456:524–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Zimmermann M, Lottersberger F, Buonomo SB, Sfeir A, de Lange T. 53BP1 regulates DSB repair using Rif1 to control 5′ end resection. Science. 2013;339:700–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang X, Tanaka T, Kurose A, Traganos F, Darzynkiewicz Z. Constitutive histone H2AX phosphorylation on Ser-139 in cells untreated by genotoxic agents is cell-cycle phase specific and attenuated by scavenging reactive oxygen species. Int J Oncol. 2006;29:495–501.PubMedGoogle Scholar
- MacPhail SH, Banáth JP, Yu Y, Chu E, Olive PL. Cell cycle-dependent expression of phosphorylated histone H2AX: reduced expression in unirradiated but not X-irradiated G1-phase cells. Radiat Res. 2003;159:759–67.View ArticlePubMedGoogle Scholar
- Dale Rein I, Stokke C, Jalal M, Myklebust JH, Patzke S, Stokke T. New distinct compartments in the G2 phase of the cell cycle defined by the levels of γH2AX. Cell Cycle. 2015;14:3261–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Suchánková J, Kozubek S, Legartová S, Sehnalová P, Küntziger T, Bártová E. Distinct kinetics of DNA repair protein accumulation at DNA lesions and cell cycle-dependent formation of γH2AX- and NBS1-positive repair foci. Biol Cell. 2015;107:440–54.View ArticlePubMedGoogle Scholar
- Šimara P, Tesařová L, Padourová S, Koutná I. Generation of human induced pluripotent stem cells using genome integrating or non-integrating methods. Folia Biol. 2014;60:85–9.Google Scholar
- Tesařová L, Simara P, Stejskal S, Koutná I. Haematopoietic developmental potential of human pluripotent stem cell lines. Folia Biol. 2014;60:90–4.Google Scholar
- Adewumi O, Aflatoonian B, Ahrlund-Richter L, Amit M, Andrews PW, Beighton G, Bello PA, Benvenisty N, Berry LS, Bevan S, et al. Characterization of human embryonic stem cell lines by the International Stem Cell Initiative. Nat Biotechnol. 2007;25:803–16.View ArticlePubMedGoogle Scholar
- Zhang M, Yang C, Liu H, Sun Y. Induced pluripotent stem cells are sensitive to DNA damage. Genomics Proteomics Bioinformatics. 2013;11:320–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Matula P, Maška M, Daněk O, Matula P, Kozubek M. Acquiarium: free software for acquisition and analysis of 3D images of cells in fluorescence microscopy. In: 6th IEEE International Symposium on Biomedical Imaging. Boston. 2009
- Štěpka K, Matula P, Wörz S, Rohr K, Kozubek M. Performance and sensitivity evaluation of 3D spot detection methods in confocal microscopy. Cytometry A. 2015;87:759–72.View ArticlePubMedGoogle Scholar
- Schultz LB, Chehab NH, Malikzay A, Halazonetis TD. p53 binding protein 1 (53BP1) is an early participant in the cellular response to DNA double-strand breaks. J Cell Biol. 2000;151:1381–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogakou EP, Boon C, Redon C, Bonner WM. Megabase chromatin domains involved in DNA double-strand breaks in vivo. J Cell Biol. 1999;146:905–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Redon CE, Dickey JS, Bonner WM, Sedelnikova OA. γ-H2AX as a biomarker of DNA damage induced by ionizing radiation in human peripheral blood lymphocytes and artificial skin. Adv Space Res. 2009;43:1171–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Salic A, Mitchison TJ. A chemical method for fast and sensitive detection of DNA synthesis in vivo. Proc Natl Acad Sci U S A. 2008;105:2415–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Stead E, White J, Faast R, Conn S, Goldstone S, Rathjen J, Dhingra U, Rathjen P, Walker D, Dalton S. Pluripotent cell division cycles are driven by ectopic Cdk2, cyclin A/E and E2F activities. Oncogene. 2002;21:8320–33.View ArticlePubMedGoogle Scholar
- Becker KA, Ghule PN, Therrien JA, Lian JB, Stein JL, van Wijnen AJ, Stein GS. Self-renewal of human embryonic stem cells is supported by a shortened G1 cell cycle phase. J Cell Physiol. 2006;209:883–93.View ArticlePubMedGoogle Scholar
- Calder A, Roth-Albin I, Bhatia S, Pilquil C, Lee JH, Bhatia M, Levadoux-Martin M, McNicol J, Russell J, Collins T, et al. Lengthened G1 phase indicates differentiation status in human embryonic stem cells. Stem Cells Dev. 2013;22:279–95.View ArticlePubMedGoogle Scholar
- Barta T, Dolezalova D, Holubcova Z, Hampl A. Cell cycle regulation in human embryonic stem cells: links to adaptation to cell culture. Exp Biol Med. 2013;238:271–5.View ArticleGoogle Scholar
- Suzuki K, Okada H, Yamauchi M, Oka Y, Kodama S, Watanabe M. Qualitative and quantitative analysis of phosphorylated ATM foci induced by low-dose ionizing radiation. Radiat Res. 2006;165:499–504.View ArticlePubMedGoogle Scholar
- Werbowetski-Ogilvie TE, Bossé M, Stewart M, Schnerch A, Ramos-Mejia V, Rouleau A, Wynder T, Smith MJ, Dingwall S, Carter T, et al. Characterization of human embryonic stem cells with features of neoplastic progression. Nat Biotechnol. 2009;27:91–7.View ArticlePubMedGoogle Scholar
- Yang S, Lin G, Tan YQ, Zhou D, Deng LY, Cheng DH, Luo SW, Liu TC, Zhou XY, Sun Z, et al. Tumor progression of culture-adapted human embryonic stem cells during long-term culture. Genes Chromosomes Cancer. 2008;47:665–79.View ArticlePubMedGoogle Scholar
- Hartwell LH, Weinert TA. Checkpoints: controls that ensure the order of cell cycle events. Science. 1989;246:629–34.View ArticlePubMedGoogle Scholar
- Kastan MB, Onyekwere O, Sidransky D, Vogelstein B, Craig RW. Participation of p53 protein in the cellular response to DNA damage. Cancer Res. 1991;51:6304–11.PubMedGoogle Scholar
- Kuerbitz SJ, Plunkett BS, Walsh WV, Kastan MB. Wild-type p53 is a cell cycle checkpoint determinant following irradiation. Proc Natl Acad Sci U S A. 1992;89:7491–5.View ArticlePubMedPubMed CentralGoogle Scholar
- González F, Georgieva D, Vanoli F, Shi ZD, Stadtfeld M, Ludwig T, Jasin M, Huangfu D. Homologous recombination DNA repair genes play a critical role in reprogramming to a pluripotent state. Cell Rep. 2013;3:651–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Raschke S, Spickermann S, Toncian T, Swantusch M, Boeker J, Giesen U, Iliakis G, Willi O, Boege F. Ultra-short laser-accelerated proton pulses have similar DNA-damaging effectiveness but produce less immediate nitroxidative stress than conventional proton beams. Sci Rep. 2016;6:32441.View ArticlePubMedPubMed CentralGoogle Scholar
- Woodward WA, Chen MS, Behbod F, Alfaro MP, Buchholz TA, Rosen JM. WNT/beta-catenin mediates radiation resistance of mouse mammary progenitor cells. Proc Natl Acad Sci U S A. 2007;104:618–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu T, Liu Y, Wen D, Tseng Z, Tahmasian M, Zhong M, Rafii S, Stadtfeld M, Hochedlinger K, Xiao A. Histone variant H2A.X deposition pattern serves as a functional epigenetic mark for distinguishing the developmental potentials of iPSCs. Cell Stem Cell. 2014;15:281–94.View ArticlePubMedGoogle Scholar
- Hussein SM, Batada NN, Vuoristo S, Ching RW, Autio R, Närvä E, Ng S, Sourour M, Hämäläinen R, Olsson C, et al. Copy number variation and selection during reprogramming to pluripotency. Nature. 2011;471:58–62.View ArticlePubMedGoogle Scholar
- Kulkarni A, Zschenker O, Reynolds G, Miller D, Murnane JP. Effect of telomere proximity on telomere position effect, chromosome healing, and sensitivity to DNA double-strand breaks in a human tumor cell line. Mol Cell Biol. 2010;30:578–89.View ArticlePubMedGoogle Scholar
- Zschenker O, Kulkarni A, Miller D, Reynolds GE, Granger-Locatelli M, Pottier G, Sabatier L, Murnane JP. Increased sensitivity of subtelomeric regions to DNA double-strand breaks in a human cancer cell line. DNA Repair. 2009;8:886–900.View ArticlePubMedPubMed CentralGoogle Scholar