The distribution of genomic variations in human iPSCs is related to replication-timing reorganization during reprogramming

Incomplete reprogramming of DNA replication timing in induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) are the foundation of cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which could affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing, a process linked to both genome regulation and genome stability, is efficiently reprogrammed to the embryonic state. To answer this, we compare genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicate their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibits delayed replication at heterochromatic regions containing genes downregulated in iPSCs with incompletely reprogrammed DNA methylation. DNA replication delays are not the result of gene expression or DNA methylation aberrations and persist after cells differentiate to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and influence the quality of iPSCs.

Incomplete Reprogramming of DNA Replication Timing in Induced Pluripotent Stem Cells

Induced pluripotent stem cells (iPSC) are a widely used cell system and a foundation for cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which have the potential to affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing - a process linked to both genome regulation and genome stability - is efficiently reprogrammed to the embryonic state. To answer this, we profiled and compared genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicated their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibit delayed replication at heterochromatic regions containing genes downregulated in iPSC with incompletely reprogrammed DNA methylation. DNA replication delays were not the result of gene expression and DNA methylation aberrations and persisted after differentiating cells to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and lead to undesirable phenotypes in iPSCs, establishing it as an important genomic feature to consider when evaluating iPSC lines.

Locus-specific transcription silencing at the FHIT gene suppresses replication stress-induced copy number variant formation and associated replication delay

Impaired replication progression leads to de novo copy number variant (CNV) formation at common fragile sites (CFSs). We previously showed that these hotspots for genome instability reside in late-replicating domains associated with large transcribed genes and provided indirect evidence that transcription is a factor in their instability. Here, we compared aphidicolin (APH)-induced CNV and CFS frequency between wild-type and isogenic cells in which FHIT gene transcription was ablated by promoter deletion. Two promoter-deletion cell lines showed reduced or absent CNV formation and CFS expression at FHIT despite continued instability at the NLGN1 control locus. APH treatment led to critical replication delays that remained unresolved in G2/M in the body of many, but not all, large transcribed genes, an effect that was reversed at FHIT by the promoter deletion. Altering RNase H1 expression did not change CNV induction frequency and DRIP-seq showed a paucity of R-loop formation in the central regions of large genes, suggesting that R-loops are not the primary mediator of the transcription effect. These results demonstrate that large gene transcription is a determining factor in replication stress-induced genomic instability and support models that CNV hotspots mainly result from the transcription-dependent passage of unreplicated DNA into mitosis.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

Detection of base analogs incorporated during DNA replication by nanopore sequencing

DNA synthesis is a fundamental requirement for cell proliferation and DNA repair, but no single method can identify the location, direction and speed of replication forks with high resolution. Mammalian cells have the ability to incorporate thymidine analogs along with the natural A, T, G and C bases during DNA synthesis, which allows for labeling of replicating or repaired DNA. Here, we demonstrate the use of the Oxford Nanopore Technologies MinION to detect 11 different thymidine analogs including CldU, BrdU, IdU as well as EdU alone or coupled to Biotin and other bulky adducts in synthetic DNA templates. We also show that the large adduct Biotin can be distinguished from the smaller analog IdU, which opens the possibility of using analog combinations to identify the location and direction of DNA synthesis. Furthermore, we detect IdU label on single DNA molecules in the genome of mouse pluripotent stem cells and using CRISPR/Cas9-mediated enrichment, determine replication rates using newly synthesized DNA strands in human mitochondrial DNA. We conclude that this novel method, termed Replipore sequencing, has the potential for on target examination of DNA replication in a wide range of biological contexts.

Copy number variant hotspots in Han Taiwanese population induced pluripotent stem cell lines - lessons from establishing the Taiwan human disease iPSC Consortium Bank

Background

The Taiwan Human Disease iPSC Service Consortium was established to accelerate Taiwan’s growing stem cell research initiatives and provide a platform for researchers interested in utilizing induced pluripotent stem cell (iPSC) technology. The consortium has generated and characterized 83 iPSC lines: 11 normal and 72 disease iPSC lines covering 21 different diseases, several of which are of high incidence in Taiwan. Whether there are any reprogramming-induced recurrent copy number variant (CNV) hotspots in iPSCs is still largely unknown.

Methods

We performed genome-wide copy number variant screening of 83 Han Taiwanese iPSC lines and compared them with 1093 control subjects using an Affymetrix genome-wide human SNP array.

Results

In the iPSCs, we identified ten specific CNV loci and seven “polymorphic” CNV regions that are associated with the reprogramming process. Additionally, we established several differentiation protocols for our iPSC lines. We demonstrated that our iPSC-derived cardiomyocytes respond to pharmacological agents and were successfully engrafted into the mouse myocardium demonstrating their potential application in cell therapy.

Conclusions

The CNV hotspots induced by cell reprogramming have successfully been identified in the current study. This finding may be used as a reference index for evaluating iPSC quality for future clinical applications. Our aim was to establish a national iPSC resource center generating iPSCs, made available to researchers, to benefit the stem cell community in Taiwan and throughout the world.

Control of DNA replication timing in the 3D genome

The 3D organization of mammalian chromatin was described more than 30 years ago by visualizing sites of DNA synthesis at different times during the S phase of the cell cycle. These early cytogenetic studies revealed structurally stable chromosome domains organized into subnuclear compartments. Active-gene-rich domains in the nuclear interior replicate early, whereas more condensed chromatin domains that are largely at the nuclear and nucleolar periphery replicate later. During the past decade, this spatiotemporal DNA replication programme has been mapped along the genome and found to correlate with epigenetic marks, transcriptional activity and features of 3D genome architecture such as chromosome compartments and topologically associated domains. But the causal relationship between these features and DNA replication timing and the regulatory mechanisms involved have remained an enigma. The recent identification of cis-acting elements regulating the replication time and 3D architecture of individual replication domains and of long non-coding RNAs that coordinate whole chromosome replication provide insights into such mechanisms.

Insights into the Link between the Organization of DNA Replication and the Mutational Landscape

The generation of a complete and accurate copy of the genetic material during each cell cycle is integral to cell growth and proliferation. However, genetic diversity is essential for adaptation and evolution, and the process of DNA replication is a fundamental source of mutations. Genome alterations do not accumulate randomly, with variations in the types and frequencies of mutations that arise in different genomic regions. Intriguingly, recent studies revealed a striking link between the mutational landscape of a genome and the spatial and temporal organization of DNA replication, referred to as the replication program. In our review, we discuss how this program may contribute to shaping the profile and spectrum of genetic alterations, with implications for genome dynamics and organismal evolution in natural and pathological contexts.

Linking the organization of DNA replication with genome maintenance

The spatial and temporal organization of genome duplication, also referred to as the replication program, is defined by the distribution and the activities of the sites of replication initiation across the genome. Alterations to the replication profile are associated with cell fate changes during development and in pathologies, but the importance of undergoing S phase with distinct and specific programs remains largely unexplored. We have recently addressed this question, focusing on the interplay between the replication program and genome maintenance. In particular, we demonstrated that when cells encounter challenges to DNA synthesis, the organization of DNA replication drives the response to replication stress that is mediated by the ATR/Rad3 checkpoint pathway, thus shaping the pattern of genome instability along the chromosomes. In this review, we present the major findings of our study and discuss how they may bring new perspectives to our understanding of the biological importance of the replication program.

New considerations for hiPSC-based models of neuropsychiatric disorders

The development of human-induced pluripotent stem cells (hiPSCs) has made possible patient-specific modeling across the spectrum of human disease. Here, we discuss recent advances in psychiatric genomics and post-mortem studies that provide critical insights concerning cell-type composition and sample size that should be considered when designing hiPSC-based studies of complex genetic disease. We review recent hiPSC-based models of SZ, in light of our new understanding of critical power limitations in the design of hiPSC-based studies of complex genetic disorders. Three possible solutions are a movement towards genetically stratified cohorts of rare variant patients, application of CRISPR technologies to engineer isogenic neural cells to study the impact of common variants, and integration of advanced genetics and hiPSC-based datasets in future studies. Overall, we emphasize that to advance the reproducibility and relevance of hiPSC-based studies, stem cell biologists must contemplate statistical and biological considerations that are already well accepted in the field of genetics. We conclude with a discussion of the hypothesis of biological convergence of disease-through molecular, cellular, circuit, and patient level phenotypes-and how this might emerge through hiPSC-based studies.

Perturbations in the Replication Program Contribute to Genomic Instability in Cancer

Cancer and genomic instability are highly impacted by the deoxyribonucleic acid (DNA) replication program. Inaccuracies in DNA replication lead to the increased acquisition of mutations and structural variations. These inaccuracies mainly stem from loss of DNA fidelity due to replication stress or due to aberrations in the temporal organization of the replication process. Here we review the mechanisms and impact of these major sources of error to the replication program.

Divergent Levels of Marker Chromosomes in an hiPSC-Based Model of Psychosis

In the process of generating presumably clonal human induced pluripotent stem cells (hiPSCs) from two carriers of a complex structural rearrangement, each having a psychotic disorder, we also serendipitously generated isogenic non-carrier control hiPSCs, finding that the rearrangement occurs as an extrachromosomal marker (mar) element. All confirmed carrier hiPSCs and differentiated neural progenitor cell lines were found to be mosaic. We caution that mar elements may be difficult to functionally evaluate in hiPSC cultures using currently available methods, as it is difficult to distinguish cells with and without mar elements in live mosaic cultures.

Kataegis Expression Signature in Breast Cancer Is Associated with Late Onset, Better Prognosis, and Higher HER2 Levels

Kataegis is a mutational process observed in ∼55% of breast tumors that results in hypermutation in localized genomic regions. Using whole-genome sequence data of 97 tumors, we examined the distribution of kataegis loci, showing that these somatic mutations are over-represented on chromosomes 8, 17, and 22 and enriched in genic regions and active chromatin elements. We show that tumors harboring kataegis are associated with transcriptome-wide expression changes consistent with low invasive potential. We exploit the kataegis expression signature to predict kataegis status in 412 breast cancers with transcriptome but not whole-genome sequence data and show that kataegis loci are enriched in high-grade, HER2(+) tumors in patients diagnosed with breast cancer at an older age and who have a later age at death. Our study demonstrates that kataegis loci are associated with important clinical features in breast cancer and may serve as a marker of good prognosis.

The inhibition of TDP-43 mitochondrial localization blocks its neuronal toxicity

Genetic mutations in TAR DNA-binding protein 43 (TARDBP, also known as TDP-43) cause amyotrophic lateral sclerosis (ALS), and an increase in the presence of TDP-43 (encoded by TARDBP) in the cytoplasm is a prominent histopathological feature of degenerating neurons in various neurodegenerative diseases. However, the molecular mechanisms by which TDP-43 contributes to ALS pathophysiology remain elusive. Here we have found that TDP-43 accumulates in the mitochondria of neurons in subjects with ALS or frontotemporal dementia (FTD). Disease-associated mutations increase TDP-43 mitochondrial localization. In mitochondria, wild-type (WT) and mutant TDP-43 preferentially bind mitochondria-transcribed messenger RNAs (mRNAs) encoding respiratory complex I subunits ND3 and ND6, impair their expression and specifically cause complex I disassembly. The suppression of TDP-43 mitochondrial localization abolishes WT and mutant TDP-43-induced mitochondrial dysfunction and neuronal loss, and improves phenotypes of transgenic mutant TDP-43 mice. Thus, our studies link TDP-43 toxicity directly to mitochondrial bioenergetics and propose the targeting of TDP-43 mitochondrial localization as a promising therapeutic approach for neurodegeneration.

Influence of ATM-Mediated DNA Damage Response on Genomic Variation in Human Induced Pluripotent Stem Cells

Genome instability is a potential limitation to the research and therapeutic application of induced pluripotent stem cells (iPSCs). Observed genomic variations reflect the combined activities of DNA damage, cellular DNA damage response (DDR), and selection pressure in culture. To understand the contribution of DDR on the distribution of copy number variations (CNVs) in iPSCs, we mapped CNVs of iPSCs with mutations in the central DDR gene ATM onto genome organization landscapes defined by genome-wide replication timing profiles. We show that following reprogramming the early and late replicating genome is differentially affected by CNVs in ATM-deficient iPSCs relative to wild-type iPSCs. Specifically, the early replicating regions had increased CNV losses during retroviral (RV) reprogramming. This differential CNV distribution was not present after later passage or after episomal reprogramming. Comparison of different reprogramming methods in the setting of defective DDR reveals unique vulnerability of early replicating open chromatin to RV vectors.

A genomic view on epilepsy and autism candidate genes

Epilepsy is a common complex disorder most frequently associated with psychiatric and neurological diseases. Massive parallel sequencing of individual or cohort genomes and exomes led the identification of several disease associated genes. We review here the candidate genes in epilepsy genetics with focus on exome and gene panel data. Together with the examination of brain expressed genes and post synaptic proteome the results show that: (1) Non-metabolic epilepsies and autism candidate genes tend to be AT-rich and (2) large transcript size and local AT-richness are characteristic features of genes involved in developmental brain disorders and synaptic functions. These results point to the preferential location of core epilepsy and autism candidate genes in late replicating, GC-poor chromosomal regions (isochores). These results indicate that the genomic alterations leading to some brain disorders are confined to responsive chromatin areas harboring brain critical genes.

Passage number is a major contributor to genomic structural variations in mouse iPSCs

Emergence of genomic instability is a practical issue in preparing neural stem cells (NSCs) and induced pluripotent stem cells (iPSCs). However, it is still not fully understood what the origins and mechanisms for formation are for the genomic alternations observed. Here, we studied the extent of genomic variation on the scale of individual cells originating from the same animal. We used mouse NSCs grown from embryonic cells and iPSCs generated from embryonic brain cells, B cells or fibroblasts, and performed comparative analysis with cultures of fibroblasts from the same mouse. In the first passage of these cell lines, aneuploidies were only observed for chromosomes 6, 11, 12, 19, and Y, which is overall at a rate lower than previously reported; de novo copy number variations (CNVs) were observed in 4.3% of neural iPSCs, 29% of B cell iPSCs, 10% of fibroblast iPSCs, and 1.3% of neurospheres. In contrast, propagation of these first passage cells to a later passage induced additional aneuploidies and CNVs. Breakpoint sequencing analysis suggested that the majority of the detected CNVs arose by replicative mechanisms. Interestingly, we detected identical de novo CNVs in different single cell colonies that appeared to have arisen independently from each other, which suggests a novel CNV formation mechanism in these cells. Our findings provide insights into mechanisms of CNV formation during reprogramming and suggest that replicative mechanisms for CNV formation accompany mitotic divisions.

From "directed differentiation" to "neuronal induction": modeling neuropsychiatric disease

Aberrant behavior and function of neurons are believed to be the primary causes of most neurological diseases and psychiatric disorders. Human postmortem samples have limited availability and, while they provide clues to the state of the brain after a prolonged illness, they offer limited insight into the factors contributing to disease onset. Conversely, animal models cannot recapitulate the polygenic origins of neuropsychiatric disease. Novel methods, such as somatic cell reprogramming, deliver nearly limitless numbers of pathogenic human neurons for the study of the mechanism of neuropsychiatric disease initiation and progression. First, this article reviews the advent of human induced pluripotent stem cell (hiPSC) technology and introduces two major methods, “directed differentiation” and “neuronal induction,” by which it is now possible to generate neurons for modeling neuropsychiatric disease. Second, it discusses the recent applications, and the limitations, of these technologies to in vitro studies of psychiatric disorders.

Unique features of mutations revealed by sequentially reprogrammed induced pluripotent stem cells

Although viable mice can be generated from induced pluripotent stem cells (iPSCs), the impact of accumulated mutations on the developmental potential of the resulting iPSCs remains to be determined. Here, we demonstrate that all-iPSC mice generated through tetraploid blastocysts complementation can tolerate the accumulation of somatic mutations for up to six generations using a Tet-on inducible reprogramming system. But, the viability of the all-iPS mice decreased with increasing generations. A whole-genome sequencing survey revealed that thousands of single-nucleotide variations (SNVs), including 44 non-synonymous ones, accumulated throughout the sequential reprogramming process. Subsequent analysis provides evidence that these accumulated SNVs account for the gradual reduction in viability of the resultant all-iPSC mice. Unexpectedly, our present reprogramming system revealed that pluripotent stem cells are heterogeneous in terms of possessing a set of copy-number alterations (CNAs). These CNAs are unique for pluripotent cells and subsequently disappear in the differentiating progenies.

Genome architecture and its roles in human copy number variation

Besides single-nucleotide variants in the human genome, large-scale genomic variants, such as copy number variations (CNVs), are being increasingly discovered as a genetic source of human diversity and the pathogenic factors of diseases. Recent experimental findings have shed light on the links between different genome architectures and CNV mutagenesis. In this review, we summarize various genomic features and discuss their contributions to CNV formation. Genomic repeats, including both low-copy and high-copy repeats, play important roles in CNV instability, which was initially known as DNA recombination events. Furthermore, it has been found that human genomic repeats can also induce DNA replication errors and consequently result in CNV mutations. Some recent studies showed that DNA replication timing, which reflects the high-order information of genomic organization, is involved in human CNV mutations. Our review highlights that genome architecture, from DNA sequence to high-order genomic organization, is an important molecular factor in CNV mutagenesis and human genomic instability.

Abnormal dosage of ultraconserved elements is highly disfavored in healthy cells but not cancer cells

Ultraconserved elements (UCEs) are strongly depleted from segmental duplications and copy number variations (CNVs) in the human genome, suggesting that deletion or duplication of a UCE can be deleterious to the mammalian cell. Here we address the process by which CNVs become depleted of UCEs. We begin by showing that depletion for UCEs characterizes the most recent large-scale human CNV datasets and then find that even newly formed de novo CNVs, which have passed through meiosis at most once, are significantly depleted for UCEs. In striking contrast, CNVs arising specifically in cancer cells are, as a rule, not depleted for UCEs and can even become significantly enriched. This observation raises the possibility that CNVs that arise somatically and are relatively newly formed are less likely to have established a CNV profile that is depleted for UCEs. Alternatively, lack of depletion for UCEs from cancer CNVs may reflect the diseased state. In support of this latter explanation, somatic CNVs that are not associated with disease are depleted for UCEs. Finally, we show that it is possible to observe the CNVs of induced pluripotent stem (iPS) cells become depleted of UCEs over time, suggesting that depletion may be established through selection against UCE-disrupting CNVs without the requirement for meiotic divisions.