Delayed DNA replication in haploid human embryonic stem cells

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in  >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Parent-of-origin-specific DNA replication timing is confined to large imprinted regions

Genomic imprinting involves differential DNA methylation and gene expression between homologous paternal and maternal loci. It remains unclear, however, whether DNA replication also shows parent-of-origin-specific patterns at imprinted or other genomic regions. Here, we investigate genome-wide asynchronous DNA replication utilizing uniparental human embryonic stem cells containing either maternal-only (parthenogenetic) or paternal-only (androgenetic) DNA. Four clusters of imprinted genes exhibited differential replication timing based on parent of origin, while the remainder of the genome, 99.82%, showed no significant replication asynchrony between parental origins. Active alleles in imprinted gene clusters replicated earlier than their inactive counterparts. At the Prader-Willi syndrome locus, replication asynchrony spanned virtually the entirety of S phase. Replication asynchrony was carried through differentiation to neuronal precursor cells in a manner consistent with gene expression. This study establishes asynchronous DNA replication as a hallmark of large imprinted gene clusters.

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.

Single-cell DNA replication dynamics in genomically unstable cancers

Dysregulated DNA replication is both a cause and a consequence of aneuploidy, yet the dynamics of DNA replication in aneuploid cell populations remains understudied. We developed a new method, PERT, for inferring cell-specific DNA replication states from single-cell whole genome sequencing, and investigated clone-specific DNA replication dynamics in >50,000 cells obtained from a collection of aneuploid and clonally heterogeneous cell lines, xenografts and primary cancer tissues. Clone replication timing (RT) profiles correlated with future copy number changes in serially passaged cell lines. Cell type was the strongest determinant of RT heterogeneity, while whole genome doubling and mutational process were associated with accumulation of late S-phase cells and weaker RT associations. Copy number changes affecting chromosome X had striking impact on RT, with loss of the inactive X allele shifting replication earlier, and loss of inactive Xq resulting in reactivation of Xp. Finally, analysis of time series xenografts illustrate how cell cycle distributions approximate clone proliferation, recapitulating expected relationships between proliferation and fitness in treatment-naive and chemotherapeutic contexts.

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.

Generation, genomic characterization, and differentiation of triploid human embryonic stem cells

Humans are diploid organisms, and triploidy in human embryos is responsible for ∼10% of spontaneous miscarriages. Surprisingly, some pregnancies proceed to triploid newborns that suffer from many neuro-developmental disorders. To investigate the impact of triploidy on human development, we generate triploid human embryonic stem cells (hESCs) by fusing isogenic haploid and diploid hESCs. Comparison of the transcriptome, methylome, and genome-wide replication timing shows general similarity between diploid and triploid hESCs. However, triploid cells have a larger volume than diploid cells, demonstrating decreased surface-area-to-volume ratio. This leads to a significant downregulation of cell surface ion channel genes, which are more essential in neural progenitors than in undifferentiated cells, leading to inhibition of differentiation, and it affects the neuronal differentiation ability of triploid hESCs, both in vitro and in vivo. Notably, our research establishes a platform to study triploidy in humans and points to their pathology as observed in triploid embryos.

Comprehensive analysis of DNA replication timing across 184 cell lines suggests a role for MCM10 in replication timing regulation

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.