Genomic methods for measuring DNA replication dynamics

Genome replication in asynchronously growing microbial populations

Biological cells replicate their genomes in a well-planned manner. The DNA replication program of an organism determines the timing at which different genomic regions are replicated, with fundamental consequences for cell homeostasis and genome stability. In a growing cell culture, genomic regions that are replicated early should be more abundant than regions that are replicated late. This abundance pattern can be experimentally measured using deep sequencing. However, a general quantitative theory linking this pattern to the replication program is still lacking. In this paper, we predict the abundance of DNA fragments in asynchronously growing cultures from any given stochastic model of the DNA replication program. As key examples, we present stochastic models of the DNA replication programs in budding yeast and Escherichia coli. In both cases, our model results are in excellent agreement with experimental data and permit to infer key information about the replication program. In particular, our method is able to infer the locations of known replication origins in budding yeast with high accuracy. These examples demonstrate that our method can provide insight into a broad range of organisms, from bacteria to eukaryotes.

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.

H4S47 O-GlcNAcylation regulates the activation of mammalian replication origins

The transmission and maintenance of genetic information in eukaryotic cells relies on the faithful duplication of the entire genome. In each round of division, excessive replication origins are licensed, with only a fraction activated to give rise to bi-directional replication forks in the context of chromatin. However, it remains elusive how eukaryotic replication origins are selectively activated. Here we demonstrate that O-GlcNAc transferase (OGT) enhances replication initiation by catalyzing H4S47 O-GlcNAcylation. Mutation of H4S47 impairs DBF4-dependent protein kinase (DDK) recruitment on chromatin, causing reduced phosphorylation of the replicative helicase mini-chromosome maintenance (MCM) complex and compromised DNA unwinding. Our short nascent-strand sequencing results further confirm the importance of H4S47 O-GlcNAcylation in origin activation. We propose that H4S47 O-GlcNAcylation directs origin activation through facilitating MCM phosphorylation, and this may shed light on the control of replication efficiency by chromatin environment.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

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.

Telomere-to-telomere human DNA replication timing profiles

The spatiotemporal organization of DNA replication produces a highly robust and reproducible replication timing profile. Sequencing-based methods for assaying replication timing genome-wide have become commonplace, but regions of high repeat content in the human genome have remained refractory to analysis. Here, we report the first nearly-gapless telomere-to-telomere replication timing profiles in human, using the T2T-CHM13 genome assembly and sequencing data for five cell lines. We find that replication timing can be successfully assayed in centromeres and large blocks of heterochromatin. Centromeric regions replicate in mid-to-late S-phase and contain replication-timing peaks at a similar density to other genomic regions, while distinct families of heterochromatic satellite DNA differ in their bias for replicating in late S-phase. The high degree of consistency in centromeric replication timing across chromosomes within each cell line prompts further investigation into the mechanisms dictating that some cell lines replicate their centromeres earlier than others, and what the consequences of this variation are.

High-throughput analysis of single human cells reveals the complex nature of DNA replication timing control

DNA replication initiates from replication origins firing throughout S phase. Debate remains about whether origins are a fixed set of loci, or a loose agglomeration of potential sites used stochastically in individual cells, and about how consistent their firing time is. We develop an approach to profile DNA replication from whole-genome sequencing of thousands of single cells, which includes in silico flow cytometry, a method for discriminating replicating and non-replicating cells. Using two microfluidic platforms, we analyze up to 2437 replicating cells from a single sample. The resolution and scale of the data allow focused analysis of replication initiation sites, demonstrating that most occur in confined genomic regions. While initiation order is remarkably similar across cells, we unexpectedly identify several subtypes of initiation regions in late-replicating regions. Taken together, high throughput, high resolution sequencing of individual cells reveals previously underappreciated variability in replication initiation and progression.

Mapping Replication Timing in Single Mammalian Cells

Replication timing (RT) is the temporal order in which genomic DNA is replicated during S phase. Early and late replication correlate with transcriptionally active and inactive chromatin compartments, but mechanistic links between large-scale chromosome structure, transcription, and replication are still enigmatic. A proper RT program is necessary to maintain the global epigenome that defines cell identity, suggesting that RT is critical for epigenome integrity by facilitating the assembly of different types of chromatin at different times during S phase. RT is regulated during development and has been found to be altered in disease. Thus, RT can identify stable epigenetic differences distinguishing cell types, and can be used to help stratify patient outcomes and identify markers of disease. Most methods to profile RT require thousands of S-phase cells. In cases where cells are rare (e.g., early-stage embryos or rare primary cell types) or consist of a heterogeneous mixture of cell states (e.g., differentiation intermediates), or when the interest is in determining the degree of stable epigenetic heterogeneity within a population of cells, single-cell measurements of RT are necessary. We have previously developed single cell Repli-seq, a method to measure replication timing in single cells using DNA copy number quantification. To date, however, single-cell Repli-seq suffers from relatively low throughput and high costs. Here, we describe an improved single-cell Repli-seq protocol that uses degenerate oligonucleotide-primed PCR (DOP-PCR) for uniform whole-genome amplification and uniquely barcoded primers that permit early pooling of single-cell samples into a single library preparation. We also provide a bioinformatics platform for analysis of the data. The improved throughput and decreased costs of this method relative to previously published single-cell Repli-seq protocols should make it considerably more accessible to a broad range of investigators. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Whole Genome Amplification (WGA) of single cells and sequence library construction. Basic Protocol 2: Deriving and displaying single-cell replication timing data from whole genome sequencing.

The Clinical Significance and Potential Molecular Mechanism of Upregulated CDC28 Protein Kinase Regulatory Subunit 1B in Osteosarcoma

BACKGROUND

CDC28 Protein Kinase Regulatory Subunit 1B (CKS1B) is a member of cyclin-dependent kinase subfamily and the relationship between CKS1B and osteosarcoma (OS) remains to be explored.

METHODS

80 OS and 41 nontumor tissue samples were arranged to conduct immunohistochemistry (IHC) to evaluate CKS1B expression between OS and nontumor samples. The standard mean deviation (SMD) was calculated based on in-house IHC and tissue microarrays and exterior high-throughput datasets for further verification of CKS1B expression in OS. The effect of CKS1B expression on clinicopathological and overall survival of OS patients was measured through public high-throughput datasets, and analysis of immune infiltration and single-cell RNA-seq was applied to ascertain molecular mechanism of CKS1B in OS.

RESULTS

A total of 197 OS samples and 83 nontumor samples (including tissue and cell line) were obtained from in-house IHC, microarrays, and exterior high-throughput datasets. The analysis of integrated expression status demonstrated upregulation of CKS1B in OS (SMD = 1.38, 95% CI [0.52-2.25]) and the significant power of CKS1B expression in distinguishing OS samples from nontumor samples (Area under the Curve (AUC) = 0.89, 95% CI [0.86-0.91]). Clinicopathological and prognosis analysis indicated no remarkable significance but inference of immune infiltration and single-cell RNA-seq prompted that OS patients with overexpressed CKS1B were more likely to suffer OS metastasis while MYC Protooncogene may be the upstream regulon of CKS1B in proliferating osteoblastic OS cells.

CONCLUSIONS

In this study, sufficient evidence was provided for upregulation of CKS1B in OS. The advanced effect of CKS1B on OS progression indicates a foreground of CKS1B as a biomarker for OS.

Delayed DNA replication in haploid human embryonic stem cells

Haploid human embryonic stem cells (ESCs) provide a powerful genetic system but diploidize at high rates. We hypothesized that diploidization results from aberrant DNA replication. To test this, we profiled DNA replication timing in isogenic haploid and diploid ESCs. The greatest difference was the earlier replication of the X Chromosome in haploids, consistent with the lack of X-Chromosome inactivation. We also identified 21 autosomal regions that had delayed replication in haploids, extending beyond the normal S phase and into G2/M. Haploid-delays comprised a unique set of quiescent genomic regions that are also underreplicated in polyploid placental cells. The same delays were observed in female ESCs with two active X Chromosomes, suggesting that increased X-Chromosome dosage may cause delayed autosomal replication. We propose that incomplete replication at the onset of mitosis could prevent cell division and result in re-entry into the cell cycle and whole genome duplication.

Replication timing analysis in polyploid cells reveals Rif1 uses multiple mechanisms to promote underreplication in Drosophila

Regulation of DNA replication and copy number is necessary to promote genome stability and maintain cell and tissue function. DNA replication is regulated temporally in a process known as replication timing (RT). Rap1-interacting factor 1 (Rif1) is a key regulator of RT and has a critical function in copy number control in polyploid cells. Previously, we demonstrated that Rif1 functions with SUUR to inhibit replication fork progression and promote underreplication (UR) of specific genomic regions. How Rif1-dependent control of RT factors into its ability to promote UR is unknown. By applying a computational approach to measure RT in Drosophila polyploid cells, we show that SUUR and Rif1 have differential roles in controlling UR and RT. Our findings reveal that Rif1 acts to promote late replication, which is necessary for SUUR-dependent underreplication. Our work provides new insight into the process of UR and its links to RT.

Efficiency and equity in origin licensing to ensure complete DNA replication

The cell division cycle must be strictly regulated during both development and adult maintenance, and efficient and well-controlled DNA replication is a key event in the cell cycle. DNA replication origins are prepared in G1 phase of the cell cycle in a process known as origin licensing which is essential for DNA replication initiation in the subsequent S phase. Appropriate origin licensing includes: (1) Licensing enough origins at adequate origin licensing speed to complete licensing before G1 phase ends; (2) Licensing origins such that they are well-distributed on all chromosomes. Both aspects of licensing are critical for replication efficiency and accuracy. In this minireview, we will discuss recent advances in defining how origin licensing speed and distribution are critical to ensure DNA replication completion and genome stability.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

TIGER: inferring DNA replication timing from whole-genome sequence data

MOTIVATION

Genomic DNA replicates according to a reproducible spatiotemporal program, with some loci replicating early in S phase while others replicate late. Despite being a central cellular process, DNA replication timing studies have been limited in scale due to technical challenges.

RESULTS

We present TIGER (Timing Inferred from Genome Replication), a computational approach for extracting DNA replication timing information from whole genome sequence data obtained from proliferating cell samples. The presence of replicating cells in a biological specimen leads to non-uniform representation of genomic DNA that depends on the timing of replication of different genomic loci. Replication dynamics can hence be observed in genome sequence data by analyzing DNA copy number along chromosomes while accounting for other sources of sequence coverage variation. TIGER is applicable to any species with a contiguous genome assembly and rivals the quality of experimental measurements of DNA replication timing. It provides a straightforward approach for measuring replication timing and can readily be applied at scale.

AVAILABILITY AND IMPLEMENTATION

TIGER is available at https://github.com/TheKorenLab/TIGER.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

Positive and Negative Regulation of DNA Replication Initiation

Genomic DNA is replicated every cell cycle by the programmed activation of replication origins at specific times and chromosomal locations. The factors that define the locations of replication origins and their typical activation times in eukaryotic cells are poorly understood. Previous studies highlighted the role of activating factors and epigenetic modifications in regulating replication initiation. Here, we review the role that repressive pathways - and their alleviation - play in establishing the genomic landscape of replication initiation. Several factors mediate this repression, in particular, factors associated with inactive chromatin. Repression can support organized, yet stochastic, replication initiation, and its absence could explain instances of rapid and random replication or re-replication.

PRIM1 promotes the proliferation of hepatocellular carcinoma cells in vitro and in vivo

PRIM1 plays an important role during oncogenesis, however it has never been reported in liver cancer, and thus our objective is to explore the role of PRIM1 in liver cancer. We selected RNAseq data of 50 paired liver cancer samples from the Cancer Gene Atlas (TCGA), and then bioinformatics methods and Mann-Whitney U test were used to analyze the correlation between PRIM1 and the clinical pathological stage of liver cancer. Quantitative polymerase chain reaction (QPCR) was used to detect mRNA expression of PRIM1 in BEL-7404, BCL-7402, HepG2 and SMMC-7721 cell lines. LV-PRIM1-RNAi was transfected into BEL-7404 and SMMC-7721 cells by lentivirus, and then Celigo imaging cytometer, Caspase3/7 Assay, flow cytometry and MTT assay were used to detect the proliferation and apoptosis of BEL-7404 and SMMC-7721 cells with ≥50% gene reduction rate after lentivirus transfection detected by QPCR. BEL-7404 and SMMC-7721 carrying PRIM1 gene were used for oncogenesis in vitro to observe the weight and fluorescence intensity of the tumor. Bioinformatics method was used to obtain the information about PRIM1 gene, and the correlation between PRIM1 and clinical pathological stage of liver cancer was analyzed by Mann-Whitney U test. QPCR results showed that PRIM1 was expressed in BEL-7404, BCL-7402, HepG2 and SMMC-7721 cell lines, which was highest in BCL-7404 cell line. Celigo imaging cytometer, Caspase3/7 Assay, flow cytometry and MTT assay showed that the proliferative ability of BEL-7404 and SMMC-7721 were decreased after LV-PRIM1-RNAi transfection. Furthermore, the weight and the fluorescence intensity of the tumors in vitro formed by LV-PRIM1-RNAi cells on SCID mice were decreased. So, interference of PRIM1 expression can inhibit the proliferation of BEL-7404 and SMMC-7721 cells, as well as induce the apoptosis of liver cancer cells.

A sampling of methods to study chromosome and genome structure and function

As a scientist, one’s perspective of the human genome is informed by the way it is studied – at the level of single nucleotides, a single gene, a specific genomic region, an entire chromosome, the complete karyotype, or the nucleus that encompasses both the genome and the nuclear components that support genome structure, function, stability, and inheritance. Experimentally investigating aspects of genome structure and chromosome number and higher order packaging requires different technical approaches that offer varying levels of resolution. This special issue of Chromosome Research provides overviews of a few current methodologies to study chromosome and genome organization and function, with a particular focus on contemporary sequencing-based approaches.