Stem Cell

iPS cell generation-associated point mutations include many C > T substitutions via different cytosine modification mechanisms – Nature Communications

Time Stamp: June 10, 2024 8:00 PM
Source Node: 611455

Origins of the mutations identified in iPSC genomes

Whether or not genomic mutations occur during the cell lineage conversion from somatic cells to iPSCs is a critically important question to resolve from both a medical and a scientific point of view. There are reports that mutations detected in iPSCs are pre-existing single nucleotide variants (SNVs) and were therefore already present in the parent somatic cells, and this debate continues¹⁵. In this present study, we aimed to obtain a more definitive answer on this issue using three approaches.

We first established iPSCs from the same somatic cells using different methods and the mutations detected in each cell line were then compared. Human dermal fibroblasts derived from a single individual (18year-old female, FC-0024, Lot No. 00967, Lifeline Cell Technology, Frederick, MD, USA) were used as the parent cells and iPSCs were generated using retroviral vectors or integration-free plasmid vectors to introduce reprogramming factors^10,11. Mutations were identified by WGS analysis of three and five iPSC lines, respectively. The results showed a statistically significant difference in the number of mutations between the methods, i.e. 461.3 ± 23.1 for retro virus-mediated vs 710.0 ± 135.6 for episomal vector-mediated (Fig. 1a, b and Supplementary Data 1). This result, that the number of mutations is dependent on the gene delivery method, indicates many mutations in iPSCs cannot be explained by pre-existing SNVs.

Normally, WGS results for the parental somatic cell populations used for iPSC generation are employed as a control, or reference sequence (ref seq), for iPSC mutation analysis (Supplementary Fig. 1). However, this ref seq may not include mutations that appeared during the culture of the parental cell population, i.e., mutations that exist in only a few cells in the population, due to the detection limit problem. On the other hand, since iPSCs are derived from a single cell, it has been pointed out that mutations that are missed in the analysis of such parent somatic cell populations can be detected and recognized as mutations unique to iPSC lines¹⁵. To address this issue herein, we focused on sister iPSC lines established from the same parental somatic cells. Since mutation identification for each sister line was performed independently using the parental somatic ref seq as usual, pre-existing SNVs that were missed may be detected in sister iPSC lines and, depending on the developmental relationship between them, may be detected as common mutations between sister cell lines (i.e., shared SNVs), and this possibility increases as the number of sister lines increases.

We therefore examined the 35 sister lines we had previously established using the same parental cell population for the presence of shared SNVs. The iPS cell lines were: 33 iPS cell lines and 2 ntESC lines independently established from one type of MEF (MEF5 prepared from a single embryo) (Fig. 1c, d). To our knowledge, no comparative study of as many as 35 lines of iPS cells established from the same parental somatic cells has been reported, providing crucial information in the discussion of pre-existing SNVs.

The number of shared SNVs detected among these sister strains is highlighted in Fig. 1c in blue (e.g. there was one shared SNV between R3F-60 and NAC1-12). WGS analysis revealed a total of 25,966 SNVs (mean, 741.9 ± 320.0 per cell line) identified across 35 lines, of which only 24 (0.28%) were shared. There were six cell lines (539, 19004, R3F-70, 1170, NAC1-7, NAC1-11) in which no shared SNVs were found, several lines in which only one was detected, one pair in which four were detected (594 vs. 866), and two pairs in which five were detected (Thy-48 vs. 420 or 427). Even with the extreme cases, there were only one pair in which 11 shared SNVs was found (420 vs. 427; Fig. 1c).

Among the mutations that arise in somatic cell populations during culture, those that occur at an earlier stage should be excluded by using the ref seq data of the parental somatic cell population, since they are expected to form larger clusters within the population. Notably in this regard, even the PCR-based ultradeep sequencing test of parental cells, which can detect mutations in 1/2500–1/5000 cells within a parental somatic cell population^10,16, has not detected the mutations detected in each iPSC by our WGS mutation analysis using the parental cell ref seq data. Hence, the lack of shared SNVs even in the 35 sister clones we analyzed by WGS indicates that most of the mutations detected in these 35 iPSC clones, averaging 741.9 ± 320.0 cells/line, were not pre-existing but arose intensively over a short period of time in the final phase of the culture period of parent cell population, including cell lineage conversion step into iPSCs.

As mentioned above, standard mutation analysis of iPSCs uses the WGS data of the parent somatic cell genome as the reference sequences. However, this may cause a sensitivity problem. Hence, we further performed mutation analysis in our present study using sister iPSC lines as well as parental somatic cell for ref seq data and compared the mutations between the two. If there were many different SNVs in each cell of the parent somatic cell population, and if these variations were detected as iPS clone-specific mutations, then the ref seq information of one of the sister lines instead of the parent somatic cell should miss significantly more mutations than the ref seq data of the parental cells. This is because the sister line is of single cell origin and is therefore only a subset of the parental somatic cell population. On the other hand, the sister line WGS data should detect SNVs at very low frequencies that are not detected in the parental cell ref seq due to sensitivity issues. The number and content of mutations detected will be expected therefore to differ between the parental cell ref seq and the sister line WGS data. Furthermore, if the data from the different sister cell lines are used as the ref seq, the results should also differ.

In our specific analyses, five independent sister clones established using somatic cells from the same single individual were used¹¹. The mutations identified in each sister clone were compared with those in the parent somatic cells and with each of the other four different sister clones as references (Fig. 1e). This analysis resulted in a total of 25 (5 clones × 5 different ref seq = 25) sets of mutation analyses (Fig. 1e, upper right). The results indicated almost no differences in the outcomes of using the parent somatic cell and any of the sister clones as the reference sequences (99.0–100% concordance rate of SNVs detected; Fig. 1e). This finding is very surprising in two ways. First, the sister iPSC ref seq data, all of which are derived from a single cell, contained all the information contained in the ref seq derived from the parent somatic cell, which is a population. Second, the ref seq data for the parental cell population, in which each cell constituting the population should have many different mutations, contained the same information as the single-cell-derived iPSC-derived ref seq information, i.e., no sensitivity limitation issues were observed. These results thus indicated that the parent somatic cell population and its five independently established sister strains harbored nearly identical information as a ref seq, such that there was no heterogeneity in the parent somatic cell population genome that could explain the hundreds of variations found in iPSCs⁵.

The results of these present analyses indicate that pre-existing SNVs are present in some cells in the iPSC parent somatic cell population. Even with this apparent heterogeneity in the parental cell population in terms of mutations, its effect on correct mutation identification in iPSCs will be very low as long as bioinformatics processing using ref seq information of the parental cell population is rigorously performed.

Variations of up to 100-fold exist in the SNV numbers among human iPS cell lines

Since we confirmed that the mutations observed in iPSCs are de novo and that sister lines can also be used to generate ref seq information, we next employed our analysis system to conduct mutation analysis of 75 human iPSC lines, 53 for which WGS information on sister clones is also available in public databases, along with the WGS results of their parent somatic cells (https://www.hipsci.org/)¹⁷, and 22 of which were generated in our laboratory (Fig. 2a, Supplementary Table 1 and Supplementary Data 1).

**Fig. 2: C to T mutations in iPS cells.**

We found that 38 (50.7%) of the 75 iPSC lines analyzed had more than 1000 mutations, including 5 lines with more than 5000 mutations and 2 lines with more than 10,000 mutations, the greatest amount being 12,857 (Fig. 2a). The cell lines with more than 10,000 SNVs would be expected to have more than 100 coding sequence (CDS) mutations and had 111 and 106 such mutations indeed (Supplementary Data 1). Notably also, the 53 clones established outside of our laboratory using the Sendai virus vector were found to have particularly large numbers of mutations, while there were also large differences in the number of mutations between lines (Fig. 2a). Not only did the frequency of mutations vary by generation procedure (102.1 ± 21.2, 710 ± 135.6, 461.3 ± 23.1, and 2382.8 ± 2432.2 mutations for CB-epi, E-HDFa, R-HDFa, and HipSci, respectively; Supplementary Table 1)¹⁰, there were also large differences in the number of mutations among the sister cell lines (Supplementary Fig. 2a). On the other hand, there was no positive correlation between mutation frequency and the age of the individual from whom the parent somatic cells were derived (Supplementary Fig. 2b). Hence, our mutations results for the 75 lines tested in this analysis also supported that iPSC mutations are not pre-existing SNVs. It should be noted that the results of this analysis do not necessarily indicate that the use of Sendai virus vectors for iPSC generation causes more mutations. It is necessary to consider other factors in these lines since large differences were observed between clones established by this method.

To further explore the characteristics of mutations detected in iPSC genomes, we conducted mutation signature analysis and normalized the results by the number of each sequence present in the genome. The results revealed that CpG is a prominent mutation-prone site, with a resulting C to T transition efficiency that is more than 10-fold higher than other low-frequency sites, with the same trend observed in both human and mouse iPSC lines (Fig. 2b and Supplementary Fig. 3). Although mutations in CpG sequences are well known polymorphisms among individuals (Fig. 2c), the frequency in iPSCs was found to be higher than that among individuals (approximately 2.0- to 14.0-fold; Fig. 2d). Furthermore, in each iPSC line, C to T transitions at a CpG site (hereafter denoted as CpG C > T) correlated with the total number of mutations detected (Supplementary Fig. 4). The intragenomic distribution of mutations also presented a unique profile in iPSCs. CpG C > T mutations in the germline neatly reflects the distribution of CpG sequences within the genome, whereas those in iPSCs showed a bias toward intergenic regions (Fig. 3a). More interestingly, the amount of CpG C > T mutations, which was increased in iPSCs, showed elevated numbers arose in retrotransposons, especially in short interspersed nuclear elements (SINEs) (Fig. 3b and Supplementary Table 2). On the other hand, no CpG specific elevation of C > T mutations was observed in long interspersed nuclear elements (LINEs) and long terminal repeat (LTR)-type retrotransposons. Most interestingly, among the SINE subfamilies, one type of subfamily, AluY that is the most active¹⁸, showed a large increase of CpG C > T mutations in iPS cells (Fig. 3c, d). Considering that the C to T mutations in non-CpG sites do not show an increased mutation frequency in SINEs in iPS cells, our results seem to suggest that CpG-specific mutations in SINEs in iPS cells are caused by a different molecular mechanism than the previously known deamination of cytosine (Fig. 3b).

T transitions in retrotransposons in iPSCs.”>

The C to T transition mutation profiles of iPS cell genomes differ from those accumulated in germlines as inter-individual differences: relationship between DNA demethylation and mutation

The existence of iPSCs with more than 10,000 de novo SNVs was demonstrated in this present study, and a large number of C > T transitions were revealed among them (Fig. 2a). The number of C > T mutations showed a positive correlation with the total number of mutations. However, when the cytosine-containing sequences in which mutations occurred were divided into CpG and non-CpG, a contrasting correlation emerged i.e., the percentage of C > T transitions at non-CpG sites correlated positively with the total number of SNVs, whereas the percentage of mutations occurring at CpG sites was found to be inversely correlated (Fig. 4a, b). Given that the C > T transition at non-CpG sites is thought to be due to deamination, this finding also indicated that a mechanism other than this was involved in generating the C > T transition at CpG sites in the iPSCs with low total mutation numbers and with a small percentage of these transitions at non-CpG sites. On the other hand, the significant increase of C > T mutations at non-CpG sites in iPSCs was also particularly noteworthy. The ratio of C > T changes at non-CpG to CpG sites appeared to be different in iPSCs from that in the germline (Supplementary Fig. 5). Although the mutation profile in the germline was consistent with that in previous reports, i.e. the rate of 5mC to T in the germline was reported previously to be about four times higher than the rate of C to U deamination^19,20,21, non-CpG C > T mutations were significantly increased in our iPSCs, depending on the derivation method. The aforementioned ratio was thereby completely reversed in HipSci iPSCs compared to the germline. This observation indicates a C > T mutagenic system that is unique to iPSCs and radically different from the germline.

**Fig. 4: Relationship between total number of mutations and C to T transition (%) (at CpG vs. non-CpG sites).**

A possible CpG-specific reaction other than demamination is the methylation of DNA. Gene expression analysis during the cell lineage conversion from somatic cell to iPSC has confirmed the induction of Tet1 and Tet2, enzymes thought to be responsible for DNA demethylation, during iPSC formation^22,23,24,25. We then conducted forced expression experiments with Tet1 in iPSC generation (OSKC + Tet1)(Fig. 5a–c). As in previous studies, an increase in iPSC generation frequency was observed upon introduction of Tet1²⁵(Fig. 5d). Subsequent WGS analysis of six Tet1-iPSC lines and three control lines revealed an increase in CpG sequence-specific C to T transition frequency upon Tet1 transducion (Fig. 5e, f). In addition, we also observed an increase in mutations in SINEs (Fig. 5g). This observation seemed to reproduce the phenomenon observed in the mutation analysis of human iPSCs genome, i.e. increased mutations in AluY retrotransposons (Fig. 3c).

**Fig. 5: Effects of Tet1 on C to T mutation in iPSCs.**

Hence, although Tet1 gene transduction experiments have suggested that DNA demethylation during iPSC generation has mutagenic potential, further analysis is needed to confirm that the C > T mutations detected in iPSC genomes are indeed due to DNA demethylation. Since DNA demethylation mediated by the Tet products is a key reaction in genome reprogramming and occurs on a large scale, it is crucial to determine whether this process is mutagenic or not. One direct way to confirm this idea would be through Tet suppression. However, it has been reported that suppression of the Tet family inhibits genome reprogramming itself^24,25. Here, we also performed suppression experiments using shTet1 and shTet2, but iPSC generation itself was inhibited (Supplementary Fig. 6a). Moreover, the iPSCs that were generated under these conditions showed an incomplete morphology, especially at the beginning of their generation and SNV analysis of these imperfect iPSCs showed no difference from controls (Supplementary Fig. 6b).

This situation required an approach other than reverse genetics. We therefore investigated the findings that have strongly suggested a close relationship between DNA demethylation and CpG C > T mutations. Further analysis of the results shown in Fig. 4 confirmed the contrasting correlation i.e., the percentage of C to T transitions at non-CpG sites showed a strong positive correlation (r = 0.82) with the total number of SNVs, whereas the percentage of mutations occurring at CpG sites showed a weak negative correlation (r = −0.21) (Fig. 6a, b). Interestingly, it was further noted that in iPS lines with significantly lower total SNV counts, the ratio of CpG C > T to all mutations was significantly higher (Fig. 6a).

T in Erythroblasts-derived iPSCs.”>

Given that non-CpG C > T transitions are thought to be due to deamination, we considered that a mechanism other than deamination might be operating in iPSCs with such a low total number of mutations and in which C > T transitions at CpG sites are common, but at non-CpG sites are rare. We focused our analysis therefore on iPSCs with few mutations, especially those derived from erythroblasts and HDFs. We compared the mutation rates at DMR (differentially methylated region) sites²⁶ and other genomic regions and found that CpG C > T mutation at DMRs occurred more frequently in erythroblast-derived iPSCs (Fig. 6c). A similar trend was clearly observed in HDF-derived iPSCs, which are the second least mutated iPSCs after erythroblast-derived iPSCs, although it was weaker than in the erythroblast-derived iPSCs. In addition, the expression of the TET1 gene in their parental cells showed a similar trend, i.e., the higher its expression level, the more likely the cells were to have CpG C > T mutations in the DMR (Fig. 6d). By contrast, no similar trend was observed in mature iPSCs i.e., between erythroblast-derived and HDF-derived iPSCs, or between TET1 gene expression and the number of CpG C > Ts in HipSci iPSCs (Supplementary Fig. 7a, b). To further assess this relationship, analyses of the very beginning of the cell lineage conversion process (day 2-3), when a large-scale reorganization of genome methylation pattern is expected to take place, would also be necessary.

Further to this, the marked elevation of the C > T mutation frequency at CpG sites in SINEs (most notably in evolutionarily young families i.e. AluY families) in iPSCs, which is not seen at non-CpG sites, also appears to suggest a close relationship between DNA demethylation and this mutation (Fig. 3b–d). Indeed, upregulation of the expression of AluY families in iPSCs compared to their parental cells was also clearly demonstrated (Supplementary Fig. 8). Furthermore, the degree of upregulation differed between iPSC clones for almost exclusively AluY families (shown in orange), but not for the non-AluY retrotransposons (gray). This observation that the degree of expression induction varies with each iPSC generation event suggested a profound relationship between iPSC generation and transcriptional induction of AluY (Supplementary Fig. 8). Increased retrotransposition activity in iPSCs has recently been reported²⁷, and DNA demethylation, followed by transcription, is the first step in this reaction.

These observations suggest that during genome reprogramming, DNA demethylation reactions are closely related to CpG C > T mutations. Because of the biological significance of this, further analysis and more direct evidence is needed in the future.

Furthermore, elucidating the relationship between DNA demethylation and CpG C > T mutations in biological phenomena other than genome reprogramming will be an important issue that may provide insights into the mechanisms of carcinogenesis and other diseases. Therefore, we also analyzed the effect of forced expression of Tet1 on mutagenesis in somatic cells that are not undergoing genome reprogramming. However, WGS analysis of a Tet1-transduced somatic cell population cannot provide valid data because this cell population is very heterogeneous in terms of mutations. To accurately study the effect of the Tet1 gene on mutagenesis therefore, it will be necessary to isolate a single transgenic cell from a post-transfection cell population, and this requires a system that can culture a single cell to a sufficient population size.

To date single cell expansion from primary cell populations remain not possible in mice and very few human cell types are capable of expansion from single cells to WGS-enabling levels. We therefore used p53KO-MEFs (5 and 16) to compare Tet1-transduced and untransduced (control) cells (Supplementary Fig. 9a, b). Briefly, the cells were infected with a retroviral vector containing Tet1 gene, and single cell isolation was performed three days later. The cells were then allowed to grow for one month (and some were set aside for WGS). After this, single cell isolation was performed again, and the number of cells was allowed to grow to 500–700 cell populations and subjected to WGS. In this system therefore, Tet1 was expressed intracellularly for one month to observe its effect. As a result, a trend toward increased CpG C > T mutations was observed in Tet1-transduced cells, although this was not statistically significant (Supplementary Fig. 9c). These observations suggested that mutations may also occur during biological events other than genome reprogramming, such as carcinogenesis following deletion of the p53 gene, and regeneration events in amphibians, where a transient decrease in p53 activity occurs²⁸. This sheds new light on the molecular mechanisms underlying these important biological phenomena. In our present experiments, we could only use p53−/− cells due to technical limitations, but addressing the question of whether DNA demethylation produces mutations even under normal p53 conditions is definitely critical, as it may be strongly associated with development and differentiation.

Biological relevance of the C to T transition at CpG sites

Although we observed dramatically less mutated iPSCs than reported previously, i.e., erythroblast-derived iPSCs (about 1/5 of conventional iPSCs), even these cell lines contained several-fold more mutations than ES cells^7,10. Detailed analysis revealed that the C to T transition occurred at a significantly higher frequency in the iPSCs, especially at CpG sites. The frequency of this mutation was about 10-fold higher than similar mutations observed in germlines (Fig. 2d). We therefore discussed the biological effects of a C > T single nucleotide substitution occurring at a CpG site and examined each of the C > T mutations detected in the iPSC genomes in this present study. Four possible mechanisms of biological effects were considered as follows: 1) mutations in the coding sequences (CDS); 2) mutations affecting splicing; 3) mutations in genomic regions other than the CDS, such as intronic and intergenic regions; and 4) mutations in the CpG island (CGI). In the case of 2), no similar instances were found among the mutations detected in the iPSCs, but in the other three types of mutations were confirmed in this study.

With regards to CDS mutations described in 1) above, a total of 168 mutations were detected in the CDS of iPSC genomes in a nonsense: missense: silent ration of 8: 97: 63, of which 58 (34.5%) were substitutions of the basic amino acid arginine (Arg)(Supplementary Fig. 10a–c and Supplementary Data 2). Furthermore, eight nonsense mutations were detected that, without exception, involved an Arg (CGA) to stop codon (TGA) transition (Fig. 7a). Charged amino acid substitutions were also characteristic. The acidic amino acid glutamine (Glu) was substituted with a basic amino acid lysine (Lys) in all 11 locations where it was detected, which can cause dramatic changes in its function, such as complete inactivation. Mutations of proline (Pro) residues, which fix the protein structure, were detected at 14 sites, 11 of which (78.6%) were silent, but 3 of which involved a leucine (Leu) substitution. Pro to Leu mutations of INS gene have been reported to cause diabetes²⁹. The genes containing the CpG mutations detected in this study included three that have been previously reported to be associated with cancer: BRIP1, SETD1B, and TERT (Supplementary Fig. 10d) (https://cancer.sanger.ac.uk/census)³⁰.

T mutations.”>

Second, the C > T transitions at CpG sites to splicing in the peri-exon region is also known to have many effects, causing dramatic conformational changes in proteins such as deletion and exon skipping³¹. For the situation in 3), the critical question is whether CpG C > T mutations in genomic regions other than coding regions are associated with biological aberrations. A large number of such mutations have been reported as SNPs and are powerful tools for identifying genes or regions highly associated with disease (i.e. disease-associated SNPs). In recent years, attempts have also been made to search for causal SNPs in disease. Alsheikh et al. conducted a comprehensive literature review of non-coding SNPs and reported that 309 causal SNPs have been experimentally validated³². These included 52 individual CpG C > T mutations, many of which are in transcriptional regulatory regions and are associated with reduced gene expression. Some have even been shown to have altered transcription factor binding. In addition, CpG C > T mutations have been shown to be risk factors for type 2 diabetes, Parkinson’s disease, ischemic heart disease, etc. (Supplementary Table 3). Indeed, among the mutations detected in iPSCs in this present study, we detected a close relationship between the mutation and methylation profile. We used sister clones to compare cells with similar methylation backgrounds to study the effect of the CpG C > T mutation on methylation. As a result, an interesting pair of sister iPSC clones (nosn_1 and nosn_6) was detected; in this pair, only nosn_6, which has the CpG C > T mutation, showed decreased methylation levels around this mutation. However, this result does not necessarily indicate that this mutation affects methylation. Therefore, we further examined the methylation levels in this region in the remaining 24 sister clone pairs and found that even in a small number of pairs (2/24 pairs), there were large differences in methylation levels between the sisters (all five hypomethylated CpG sites around the mutation detected in nosn_6 had a β value of 0.2 or greater). This result indicates that this region can be methylation variable even without the identified CpG C > T mutation, suggesting that the mutation is a consequence rather than a cause of the hypomethylation. Notably, the difference in methylation between the sisters is greatest between nosn_1 and nosn_6, so the effect of CpG C > T in such a region prone to methylation fluctuations is a question worth clarifying in the future (Fig. 7b, c). Thus, a C to T transition occurring at a single CpG site, even outside the coding region, has the potential to cause abnormalities. Finally, with regard to 4), many CpG C > T mutations detected in iPSCs were located around a TSS region, where mutations are usually rare. There is also a tendency for mutations to occur less frequently closer to the TSS region, which was not observed in iPSCs (Fig. 7d, e).

SEO Powered Content & PR Distribution. Get Amplified Today.
PlatoData.Network Vertical Generative Ai. Empower Yourself. Access Here.
PlatoAiStream. Web3 Intelligence. Knowledge Amplified. Access Here.
PlatoESG. Carbon, CleanTech, Energy, Environment, Solar, Waste Management. Access Here.
PlatoHealth. Biotech and Clinical Trials Intelligence. Access Here.
Source: https://www.nature.com/articles/s41467-024-49335-5

iPS cell generation-associated point mutations include many C > T substitutions via different cytosine modification mechanisms – Nature Communications

Origins of the mutations identified in iPSC genomes

Variations of up to 100-fold exist in the SNV numbers among human iPS cell lines

The C to T transition mutation profiles of iPS cell genomes differ from those accumulated in germlines as inter-individual differences: relationship between DNA demethylation and mutation

Biological relevance of the C to T transition at CpG sites

Latest Intelligence

Systemic and local regulation of hematopoietic homeostasis in health and disease – Nature Cardiovascular Research

Mechanics of morphogenesis – Nature Methods

Mostly discouraging trial of MSC-NPs or purported neural cells from MSCs for MS – The Niche

Mostly discouraging trial of MSC-NPs or purported neural cells from MSCs for MS – The Niche

Red2Flpe-SCON: a versatile, multicolor strategy for generating mosaic conditional knockout mice – Nature Communications

Red2Flpe-SCON: a versatile, multicolor strategy for generating mosaic conditional knockout mice – Nature Communications

Embryo models need consistent ethical oversight

A histone code controls positional identity in limb regeneration – Nature Reviews Molecular Cell Biology

Recommended reads: FDA warning on stem cells for pets, chromatin, aging & HSCs – The Niche

An efficient low cost means of biophysical gene transfection in primary cells – Scientific Reports