DNA replication timing

Sex-specific DNA-replication in the early mammalian embryo

The timing of DNA replication in mammals is crucial for minimizing errors and influenced by genome usage and chromatin states. Replication timing in the newly formed mammalian embryo remains poorly understood. Here, we have investigated replication timing in mouse zygotes and 2-cell embryos, revealing that zygotes lack a conventional replication timing program, which then emerges in 2-cell embryos. This program differs from embryonic stem cells and generally correlates with transcription and genome compartmentalization of both parental genomes. However, consistent and systematic differences existed between the replication timing of the two parental genomes, including considerably later replication of maternal pericentromeric regions compared to paternal counterparts. Moreover, maternal chromatin modified by Polycomb Repressive Complexes in the oocyte, undergoes early replication, despite belonging to the typically late-replicating B-compartment of the genome. This atypical and asynchronous replication of the two parental genomes may advance our understanding of replication stress in early human embryos and trigger strategies to reduce errors and aneuploidies.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%-18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

Activity of DNA polymerase κ across the genome in human fibroblasts

DNA polymerase κ (Polκ) is a specialized polymerase that has multiple cellular roles such as translesion DNA synthesis, replication of repetitive sequences, and nucleotide excision repair. We have developed a method for capturing DNA synthesized by Polκ utilizing a Polκ-specific substrate, N2-(4-ethynylbenzyl)-2'-deoxyguanosine (EBndG). After shearing of the DNA into 200 to 500 bp lengths, the EBndG-containing DNA was covalently bound to biotin using the Cu(I)-catalyzed alkyne-azide cycloaddition reaction and isolated with streptavidin beads. Isolated DNA was then ligated to adaptors, followed by PCR amplification and next-generation sequencing to generate genome-wide repair maps. We have termed this method polymerase κ sequencing. Here, we present the human genome maps for Polκ activity in an undamaged cell line. We found that Polκ activity was enhanced in GC-rich regions, euchromatin regions, the promoter of genes, and in DNA that is replicated early in the S phase.

Pioneer factors for DNA replication initiation in metazoans

Precise DNA replication is fundamental for genetic inheritance. In eukaryotes, replication initiates at multiple origins that are first "licensed" and subsequently "fired" to activate DNA synthesis. Despite the success in identifying origins with specific DNA motifs in Saccharomyces cerevisiae, no consensus sequence or sequences with a predictive value of replication origins have been recognized in metazoan genomes. Rather, epigenetic rules and chromatin structures are believed to play important roles in governing the selection and activation of replication origins. We propose that replication initiation is facilitated by a group of sequence-specific "replication pioneer factors," which function to increase chromatin accessibility and foster a chromatin environment that is conducive to the loading of the prereplication complex. Dysregulation of the function of these factors may lead to gene duplication, genomic instability, and ultimately the occurrence of pathological conditions such as cancer.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Acute multi-level response to defective de novo chromatin assembly in S-phase

Long-term perturbation of de novo chromatin assembly during DNA replication has profound effects on epigenome maintenance and cell fate. The early mechanistic origin of these defects is unknown. Here, we combine acute degradation of Chromatin Assembly Factor 1 (CAF-1), a key player in de novo chromatin assembly, with single-cell genomics, quantitative proteomics, and live-microscopy to uncover these initiating mechanisms in human cells. CAF-1 loss immediately slows down DNA replication speed and renders nascent DNA hyperaccessible. A rapid cellular response, distinct from canonical DNA damage signaling, is triggered and lowers histone mRNAs. As a result, histone variants usage and their modifications are altered, limiting transcriptional fidelity and delaying chromatin maturation within a single S-phase. This multi-level response induces a cell-cycle arrest after mitosis. Our work reveals the immediate consequences of defective de novo chromatin assembly during DNA replication, explaining how at later times the epigenome and cell fate can be altered.

Sir2 and Fun30 regulate ribosomal DNA replication timing via Mcm helicase positioning and nucleosome occupancy

The association between late replication timing and low transcription rates in eukaryotic heterochromatin is well-known, yet the specific mechanisms underlying this link remain uncertain. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). We have previously reported that in the absence of SIR2, a derepressed RNA PolII repositions MCM replicative helicases from their loading site at the ribosomal origin, where they abut well-positioned, high-occupancy nucleosomes, to an adjacent region with lower nucleosome occupancy. By developing a method that can distinguish activation of closely spaced MCM complexes, here we show that the displaced MCMs at rDNA origins have increased firing propensity compared to the non-displaced MCMs. Furthermore, we found that both, activation of the repositioned MCMs and low occupancy of the adjacent nucleosomes critically depend on the chromatin remodeling activity of FUN30. Our study elucidates the mechanism by which Sir2 delays replication timing, and it demonstrates, for the first time, that activation of a specific replication origin in vivo relies on the nucleosome context shaped by a single chromatin remodeler.

TopBP1 utilises a bipartite GINS binding mode to support genome replication

Activation of the replicative Mcm2-7 helicase by loading GINS and Cdc45 is crucial for replication origin firing, and as such for faithful genetic inheritance. Our biochemical and structural studies demonstrate that the helicase activator GINS interacts with TopBP1 through two separate binding surfaces, the first involving a stretch of highly conserved amino acids in the TopBP1-GINI region, the second a surface on TopBP1-BRCT4. The two surfaces bind to opposite ends of the A domain of the GINS subunit Psf1. Mutation analysis reveals that either surface is individually able to support TopBP1-GINS interaction, albeit with reduced affinity. Consistently, either surface is sufficient for replication origin firing in Xenopus egg extracts and becomes essential in the absence of the other. The TopBP1-GINS interaction appears sterically incompatible with simultaneous binding of DNA polymerase epsilon (Polε) to GINS when bound to Mcm2-7-Cdc45, although TopBP1-BRCT4 and the Polε subunit PolE2 show only partial competitivity in binding to Psf1. Our TopBP1-GINS model improves the understanding of the recently characterised metazoan pre-loading complex. It further predicts the coordination of three molecular origin firing processes, DNA polymerase epsilon arrival, TopBP1 ejection and GINS integration into Mcm2-7-Cdc45.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

Tissue-Predisposition to Cancer Driver Mutations

Driver mutations are considered the cornerstone of cancer initiation. They are defined as mutations that convey a competitive fitness advantage, and hence, their mutation frequency in premalignant tissue is expected to exceed the basal mutation rate. In old terms, that translates to "the survival of the fittest" and implies that a selective process underlies the frequency of cancer driver mutations. In that sense, each tissue is its own niche that creates a molecular selective pressure that may favor the propagation of a mutation or not. At the heart of this stands one of the biggest riddles in cancer biology: the tissue-predisposition to cancer driver mutations. The frequency of cancer driver mutations among tissues is non-uniform: for instance, mutations in APC are particularly frequent in colorectal cancer, and 99% of chronic myeloid leukemia patients harbor the driver BCR-ABL1 fusion mutation, which is rarely found in solid tumors. Here, we provide a mechanistic framework that aims to explain how tissue-specific features, ranging from epigenetic underpinnings to the expression of viral transposable elements, establish a molecular basis for selecting cancer driver mutations in a tissue-specific manner.

Diffusion controls local versus dispersed inheritance of histones during replication and shapes epigenomic architecture

The dynamics of inheritance of histones and their associated modifications across cell divisions can have major consequences on maintenance of the cellular epigenomic state. Recent experiments contradict the long-held notion that histone inheritance during replication is always local, suggesting that active and repressed regions of the genome exhibit fundamentally different histone dynamics independent of transcription-coupled turnover. Here we develop a stochastic model of histone dynamics at the replication fork and demonstrate that differential diffusivity of histones in active versus repressed chromatin is sufficient to quantitatively explain these recent experiments. Further, we use the model to predict patterns in histone mark similarity between pairs of genomic loci that should be developed as a result of diffusion, but cannot originate from either PRC2 mediated mark spreading or transcriptional processes. Interestingly, using a combination of CHIP-seq, replication timing and Hi-C datasets we demonstrate that all the computationally predicted patterns are consistently observed for both active and repressive histone marks in two different cell lines. While direct evidence for histone diffusion remains controversial, our results suggest that dislodged histones in euchromatin and facultative heterochromatin may exhibit some level of diffusion within "Diffusion-Accessible-Domains" (DADs), leading to redistribution of epigenetic marks within and across chromosomes. Preservation of the epigenomic state across cell divisions therefore might be achieved not by passing on strict positional information of histone marks, but by maintaining the marks in somewhat larger DADs of the genome.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, "a site to which MCM is bound in G1" might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, a technique referred to as "Chromatin Endogenous Cleavage", we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least 3 orders of magnitude. Published data quantifying the production of ssDNA during S phase showed clear evidence of replication initiation among the most abundant 1600 of these sites, with replication activity decreasing in concert with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5,500 sites. Specifically, these sites (1) appeared in intergenic nucleosome-free regions that were flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. Furthermore, the high resolution of this technique allowed us to demonstrate a strong bias for detecting Mcm double-hexamers downstream rather than upstream of the ACS, which is consistent with the directionality of Mcm loading by Orc that has been observed in vitro. We conclude that, if sites at which Mcm double-hexamers are loaded can function as replication origins, then DNA replication origins are at least 3-fold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than is widely assumed.

Genome-material costs and functional trade-offs in the autopolyploid Solidago gigantea (giant goldenrod) series

PREMISE

Increased genome-material costs of N and P atoms inherent to organisms with larger genomes have been proposed to limit growth under nutrient scarcities and to promote growth under nutrient enrichments. Such responsiveness may reflect a nutrient-dependent diploid versus polyploid advantage that could have vast ecological and evolutionary implications, but direct evidence that material costs increase with ploidy level and/or influence cytotype-dependent growth, metabolic, and/or resource-use trade-offs is limited.

METHODS

We grew diploid, autotetraploid, and autohexaploid Solidago gigantea plants with one of four ambient or enriched N:P ratios and measured traits related to material costs, primary and secondary metabolism, and resource-use.

RESULTS

Relative to diploids, polyploids invested more N and P into cells, and tetraploids grew more with N enrichments, suggesting that material costs increase with ploidy level. Polyploids also generally exhibited strategies that could minimize material-cost constraints over both long (reduced monoploid genome size) and short (more extreme transcriptome downsizing, reduced photosynthesis rates and terpene concentrations, enhanced N-use efficiencies) evolutionary time periods. Furthermore, polyploids had lower transpiration rates but higher water-use efficiencies than diploids, both of which were more pronounced under nutrient-limiting conditions.

CONCLUSIONS

N and P material costs increase with ploidy level, but material-cost constraints might be lessened by resource allocation/investment mechanisms that can also alter ecological dynamics and selection. Our results enhance mechanistic understanding of how global increases in nutrients might provide a release from material-cost constraints in polyploids that could impact ploidy (or genome-size)-specific performances, cytogeographic patterning, and multispecies community structuring.

ATXN3 controls DNA replication and transcription by regulating chromatin structure

The deubiquitinating enzyme Ataxin-3 (ATXN3) contains a polyglutamine (PolyQ) region, the expansion of which causes spinocerebellar ataxia type-3 (SCA3). ATXN3 has multiple functions, such as regulating transcription or controlling genomic stability after DNA damage. Here we report the role of ATXN3 in chromatin organization during unperturbed conditions, in a catalytic-independent manner. The lack of ATXN3 leads to abnormalities in nuclear and nucleolar morphology, alters DNA replication timing and increases transcription. Additionally, indicators of more open chromatin, such as increased mobility of histone H1, changes in epigenetic marks and higher sensitivity to micrococcal nuclease digestion were detected in the absence of ATXN3. Interestingly, the effects observed in cells lacking ATXN3 are epistatic to the inhibition or lack of the histone deacetylase 3 (HDAC3), an interaction partner of ATXN3. The absence of ATXN3 decreases the recruitment of endogenous HDAC3 to the chromatin, as well as the HDAC3 nuclear/cytoplasm ratio after HDAC3 overexpression, suggesting that ATXN3 controls the subcellular localization of HDAC3. Importantly, the overexpression of a PolyQ-expanded version of ATXN3 behaves as a null mutant, altering DNA replication parameters, epigenetic marks and the subcellular distribution of HDAC3, giving new insights into the molecular basis of the disease.

Cell-type specificity of the human mutation landscape with respect to DNA replication dynamics

The patterns of genomic mutations are associated with various genomic features, most notably late replication timing, yet it remains contested which mutation types and signatures relate to DNA replication dynamics and to what extent. Here, we perform high-resolution comparisons of mutational landscapes between lymphoblastoid cell lines, chronic lymphocytic leukemia tumors, and three colon adenocarcinoma cell lines, including two with mismatch repair deficiency. Using cell-type-matched replication timing profiles, we demonstrate that mutation rates exhibit heterogeneous replication timing associations among cell types. This cell-type heterogeneity extends to the underlying mutational pathways, as mutational signatures show inconsistent replication timing bias between cell types. Moreover, replicative strand asymmetries exhibit similar cell-type specificity, albeit with different relationships to replication timing than mutation rates. Overall, we reveal an underappreciated complexity and cell-type specificity of mutational pathways and their relationship to replication timing.

The dynamical organization of the core pluripotency transcription factors responds to differentiation cues in early S-phase

DNA replication in stem cells is a major challenge for pluripotency preservation and cell fate decisions. This process involves massive changes in the chromatin architecture and the reorganization of many transcription-related molecules in different spatial and temporal scales. Pluripotency is controlled by the master transcription factors (TFs) OCT4, SOX2 and NANOG that partition into condensates in the nucleus of embryonic stem cells. These condensates are proposed to play relevant roles in the regulation of gene expression and the maintenance of pluripotency. Here, we asked whether the dynamical distribution of the pluripotency TFs changes during the cell cycle, particularly during DNA replication. Since the S phase is considered to be a window of opportunity for cell fate decisions, we explored if differentiation cues in G1 phase trigger changes in the distribution of these TFs during the subsequent S phase. Our results show a spatial redistribution of TFs condensates during DNA replication which was not directly related to chromatin compaction. Additionally, fluorescence fluctuation spectroscopy revealed TF-specific, subtle changes in the landscape of TF-chromatin interactions, consistent with their particularities as key players of the pluripotency network. Moreover, we found that differentiation stimuli in the preceding G1 phase triggered a relatively fast and massive reorganization of pluripotency TFs in early-S phase. Particularly, OCT4 and SOX2 condensates dissolved whereas the lifetimes of TF-chromatin interactions increased suggesting that the reorganization of condensates is accompanied with a change in the landscape of TF-chromatin interactions. Notably, NANOG showed impaired interactions with chromatin in stimulated early-S cells in line with its role as naïve pluripotency TF. Together, these findings provide new insights into the regulation of the core pluripotency TFs during DNA replication of embryonic stem cells and highlight their different roles at early differentiation stages.

Drug-induced loss of imprinting revealed using bioluminescent reporters of Cdkn1c

Genomic imprinting is an epigenetically mediated mechanism that regulates allelic expression of genes based upon parent-of-origin and provides a paradigm for studying epigenetic silencing and release. Here, bioluminescent reporters for the maternally-expressed imprinted gene Cdkn1c are used to examine the capacity of chromatin-modifying drugs to reverse paternal Cdkn1c silencing. Exposure of reporter mouse embryonic stem cells (mESCs) to 5-Azacytidine, HDAC inhibitors, BET inhibitors or GSK-J4 (KDM6A/B inhibitor) relieved repression of paternal Cdkn1c, either selectively or by inducing biallelic effects. Treatment of reporter fibroblasts with HDAC inhibitors or GSK-J4 resulted in similar paternal Cdkn1c activation, whereas BET inhibitor-induced loss of imprinting was specific to mESCs. Changes in allelic expression were generally not sustained in dividing cultures upon drug removal, indicating that the underlying epigenetic memory of silencing was maintained. In contrast, Cdkn1c de-repression by GSK-J4 was retained in both mESCs and fibroblasts following inhibitor removal, although this impact may be linked to cellular stress and DNA damage. Taken together, these data introduce bioluminescent reporter cells as tools for studying epigenetic silencing and disruption, and demonstrate that Cdkn1c imprinting requires distinct and cell-type specific chromatin features and modifying enzymes to enact and propagate a memory of silencing.

The evolution of the human DNA replication timing program

DNA is replicated according to a defined spatiotemporal program that is linked to both gene regulation and genome stability. The evolutionary forces that have shaped replication timing programs in eukaryotic species are largely unknown. Here, we studied the molecular causes and consequences of replication timing evolution across 94 humans, 95 chimpanzees, and 23 rhesus macaques. Replication timing differences recapitulated the species' phylogenetic tree, suggesting continuous evolution of the DNA replication timing program in primates. Hundreds of genomic regions had significant replication timing variation between humans and chimpanzees, of which 66 showed advances in replication origin firing in humans, while 57 were delayed. Genes overlapping these regions displayed correlated changes in expression levels and chromatin structure. Many human-chimpanzee variants also exhibited interindividual replication timing variation, pointing to ongoing evolution of replication timing at these loci. Association of replication timing variation with genetic variation revealed that DNA sequence evolution can explain replication timing variation between species. Taken together, DNA replication timing shows substantial and ongoing evolution in the human lineage that is driven by sequence alterations and could impact regulatory evolution at specific genomic sites.

Mutation Rate and Spectrum of the Silkworm in Normal and Temperature Stress Conditions

Mutation rate is a crucial parameter in evolutionary genetics. However, the mutation rate of most species as well as the extent to which the environment can alter the genome of multicellular organisms remain poorly understood. Here, we used parents–progeny sequencing to investigate the mutation rate and spectrum of the domestic silkworm (Bombyx mori) among normal and two temperature stress conditions (32 °C and 0 °C). The rate of single-nucleotide mutations in the normal temperature rearing condition was 0.41 × 10−8 (95% confidence interval, 0.33 × 10−8–0.49 × 10−8) per site per generation, which was up to 1.5-fold higher than in four previously studied insects. Moreover, the mutation rates of the silkworm under the stresses are significantly higher than in normal conditions. Furthermore, the mutation rate varies less in gene regions under normal and temperature stresses. Together, these findings expand the known diversity of the mutation rate among eukaryotes but also have implications for evolutionary analysis that assumes a constant mutation rate among species and environments.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

Origins of DNA replication in eukaryotes

Errors occurring during DNA replication can result in inaccurate replication, incomplete replication, or re-replication, resulting in genome instability that can lead to diseases such as cancer or disorders such as autism. A great deal of progress has been made toward understanding the entire process of DNA replication in eukaryotes, including the mechanism of initiation and its control. This review focuses on the current understanding of how the origin recognition complex (ORC) contributes to determining the location of replication initiation in the multiple chromosomes within eukaryotic cells, as well as methods for mapping the location and temporal patterning of DNA replication. Origin specification and configuration vary substantially between eukaryotic species and in some cases co-evolved with gene-silencing mechanisms. We discuss the possibility that centromeres and origins of DNA replication were originally derived from a common element and later separated during evolution.

The human pre-replication complex is an open complex

In eukaryotes, DNA replication initiation requires assembly and activation of the minichromosome maintenance (MCM) 2-7 double hexamer (DH) to melt origin DNA strands. However, the mechanism for this initial melting is unknown. Here, we report a 2.59-Å cryo-electron microscopy structure of the human MCM-DH (hMCM-DH), also known as the pre-replication complex. In this structure, the hMCM-DH with a constricted central channel untwists and stretches the DNA strands such that almost a half turn of the bound duplex DNA is distorted with 1 base pair completely separated, generating an initial open structure (IOS) at the hexamer junction. Disturbing the IOS inhibits DH formation and replication initiation. Mapping of hMCM-DH footprints indicates that IOSs are distributed across the genome in large clusters aligning well with initiation zones designed for stochastic origin firing. This work unravels an intrinsic mechanism that couples DH formation with initial DNA melting to license replication initiation in human cells.

Double-strand break repair and mis-repair in 3D

DNA double-strand breaks (DSBs) are lesions that arise frequently from exposure to damaging agents as well as from ongoing physiological DNA transactions. Mis-repair of DSBs leads to rearrangements and structural variations in chromosomes, including insertions, deletions, and translocations implicated in disease. The DNA damage response (DDR) limits pathologic mutations and large-scale chromosome rearrangements. DSB repair initiates in 2D at DNA lesions with the stepwise recruitment of repair proteins and local chromatin remodeling which facilitates break accessibility. More complex structures are then formed via protein assembly into nanodomains and via genome folding into chromatin loops. Subsequently, 3D reorganization of DSBs is guided by clustering forces which drive the assembly of repair domains harboring multiple lesions. These domains are further stabilized and insulated into condensates via liquid-liquid phase-separation. Here, we discuss the benefits and risks associated with this 3D reorganization of the broken genome.

Integrative chromatin domain annotation through graph embedding of Hi-C data

MOTIVATION

The organization of the genome into domains plays a central role in gene expression and other cellular activities. Researchers identify genomic domains mainly through two views: 1D functional assays such as ChIP-seq, and chromatin conformation assays such as Hi-C. Fully understanding domains requires integrative modeling that combines these two views. However, the predominant form of integrative modeling uses segmentation and genome annotation (SAGA) along with the rigid assumption that loci in contact are more likely to share the same domain type, which is not necessarily true for epigenomic domain types and genome-wide chromatin interactions.

RESULTS

Here, we present an integrative approach that annotates domains using both 1D functional genomic signals and Hi-C measurements of genome-wide 3D interactions without the use of a pairwise prior. We do so by using a graph embedding to learn structural features corresponding to each genomic region, then inputting learned structural features along with functional genomic signals to a SAGA algorithm. We show that our domain types recapitulate well-known subcompartments with an additional granularity that distinguishes a combination of the spatial and functional states of the genomic regions. In particular, we identified a division of the previously identified A2 subcompartment such that the divided domain types have significantly varying expression levels.

AVAILABILITY AND IMPLEMENTATION

https://github.com/nedashokraneh/IChDA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Drosophila SUMM4 complex couples insulator function and DNA replication control

Asynchronous replication of chromosome domains during S phase is essential for eukaryotic genome function, but the mechanisms establishing which domains replicate early versus late in different cell types remain incompletely understood. Intercalary heterochromatin domains replicate very late in both diploid chromosomes of dividing cells and in endoreplicating polytene chromosomes where they are also underreplicated. Drosophila SNF2-related factor SUUR imparts locus-specific underreplication of polytene chromosomes. SUUR negatively regulates DNA replication fork progression; however, its mechanism of action remains obscure. Here, we developed a novel method termed MS-Enabled Rapid protein Complex Identification (MERCI) to isolate a stable stoichiometric native complex SUMM4 that comprises SUUR and a chromatin boundary protein Mod(Mdg4)-67.2. Mod(Mdg4) stimulates SUUR ATPase activity and is required for a normal spatiotemporal distribution of SUUR in vivo. SUUR and Mod(Mdg4)-67.2 together mediate the activities of gypsy insulator that prevent certain enhancer-promoter interactions and establish euchromatin-heterochromatin barriers in the genome. Furthermore, SuUR or mod(mdg4) mutations reverse underreplication of intercalary heterochromatin. Thus, SUMM4 can impart late replication of intercalary heterochromatin by attenuating the progression of replication forks through euchromatin/heterochromatin boundaries. Our findings implicate a SNF2 family ATP-dependent motor protein SUUR in the insulator function, reveal that DNA replication can be delayed by a chromatin barrier, and uncover a critical role for architectural proteins in replication control. They suggest a mechanism for the establishment of late replication that does not depend on an asynchronous firing of late replication origins.

DNA replication timing: Biochemical mechanisms and biological significance

The regulation of DNA replication is a fascinating biological problem both from a mechanistic angle-How is replication timing regulated?-and from an evolutionary one-Why is replication timing regulated? Recent work has provided significant insight into the first question. Detailed biochemical understanding of the mechanism and regulation of replication initiation has made possible robust hypotheses for how replication timing is regulated. Moreover, technical progress, including high-throughput, single-molecule mapping of replication initiation and single-cell assays of replication timing, has allowed for direct testing of these hypotheses in mammalian cells. This work has consolidated the conclusion that differential replication timing is a consequence of the varying probability of replication origin initiation. The second question is more difficult to directly address experimentally. Nonetheless, plausible hypotheses can be made and one-that replication timing contributes to the regulation of chromatin structure-has received new experimental support.

Accurate Identification of DNA Replication Origin by Fusing Epigenomics and Chromatin Interaction Information

DNA replication initiation is a complex process involving various genetic and epigenomic signatures. The correct identification of replication origins (ORIs) could provide important clues for the study of a variety of diseases caused by replication. Here, we design a computational approach named iORI-Epi to recognize ORIs by incorporating epigenome-based features, sequence-based features, and 3D genome-based features. The iORI-Epi displays excellent robustness and generalization ability on both training datasets and independent datasets of K562 cell line. Further experiments confirm that iORI-Epi is highly scalable in other cell lines (MCF7 and HCT116). We also analyze and clarify the regulatory role of epigenomic marks, DNA motifs, and chromatin interaction in DNA replication initiation of eukaryotic genomes. Finally, we discuss gene enrichment pathways from the perspective of ORIs in different replication timing states and heuristically dissect the effect of promoters on replication initiation. Our computational methodology is worth extending to ORI identification in other eukaryotic species.

Exploration of Tools for the Interpretation of Human Non-Coding Variants

The advent of Whole Genome Sequencing (WGS) broadened the genetic variation detection range, revealing the presence of variants even in non-coding regions of the genome, which would have been missed using targeted approaches. One of the most challenging issues in WGS analysis regards the interpretation of annotated variants. This review focuses on tools suitable for the functional annotation of variants falling into non-coding regions. It couples the description of non-coding genomic areas with the results and performance of existing tools for a functional interpretation of the effect of variants in these regions. Tools were tested in a controlled genomic scenario, representing the ground-truth and allowing us to determine software performance.

Global early replication disrupts gene expression and chromatin conformation in a single cell cycle

Background

The early embryonic divisions of many organisms, including fish, flies, and frogs, are characterized by a very rapid S-phase caused by high rates of replication initiation. In somatic cells, S-phase is much longer due to both a reduction in the total number of initiation events and the imposition of a temporal order of origin activation. The physiological importance of changes in the rate and timing of replication initiation in S-phase remains unclear.

Results

Here we assess the importance of the temporal control of replication initiation using a conditional system in budding yeast to drive the early replication of the majority of origins in a single cell cycle. We show that global early replication disrupts the expression of over a quarter of all genes. By deleting individual origins, we show that delaying replication is sufficient to restore normal gene expression, directly implicating origin firing control in this regulation. Global early replication disrupts nucleosome positioning and transcription factor binding during S-phase, suggesting that the rate of S-phase is important to regulate the chromatin landscape.

Conclusions

Together, these data provide new insight into the role of the temporal control of origin firing during S-phase for coordinating replication, gene expression, and chromatin establishment as occurs in the early embryo.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02788-7.

Rpd3 regulates single-copy origins independently of the rDNA array by opposing Fkh1-mediated origin stimulation

The faithful replication of eukaryotic genomes requires balancing the replication capacities of different genomic regions, such as repetitive versus single-copy genetic elements, which may compete for limiting replication resources, possibly leading to replication stress and genome instability. We examined the function of histone deacetylases Rpd3 and Sir2 in balancing replication between unique genome sequences and the multicopy ribosomal DNA genes. Our findings support prior conclusions that Sir2 directly suppresses early firing of rDNA origins, thereby enabling balanced replication of the genome. We further show that Rpd3’s function in delaying firing of later-firing, single-copy origins is independent of Sir2 and rDNA load. Instead, Rpd3 appears to oppose the Fkh1/2 origin activation pathway by regulating binding of the origin-stimulator Fkh1.

Eukaryotic chromosomes are organized into structural and functional domains with characteristic replication timings, which are thought to contribute to epigenetic programming and genome stability. Differential replication timing results from epigenetic mechanisms that positively and negatively regulate the competition for limiting replication initiation factors. Histone deacetylase Sir2 negatively regulates initiation of the multicopy (∼150) rDNA origins, while Rpd3 histone deacetylase negatively regulates firing of single-copy origins. However, Rpd3’s effect on single-copy origins might derive indirectly from a positive function for Rpd3 in rDNA origin firing shifting the competitive balance. Our quantitative experiments support the idea that origins compete for limiting factors; however, our results show that Rpd3’s effect on single-copy origin is independent of rDNA copy-number and of Sir2’s effects on rDNA origin firing. Whereas RPD3 deletion and SIR2 deletion alter the early S phase dynamics of single-copy and rDNA origin firings in opposite fashion, unexpectedly only RPD3 deletion suppresses overall rDNA origin efficiency across S phase. Increased origin activation in rpd3Δ requires Fkh1/2, suggesting that Rpd3 opposes Fkh1/2-origin stimulation, which involves recruitment of Dbf4-dependent kinase (DDK). Indeed, Fkh1 binding increases at Rpd3-regulated origins in rpd3Δ cells in G1, supporting a mechanism whereby Rpd3 influences initiation timing of single-copy origins directly through modulation of Fkh1-origin binding. Genetic suppression of a DBF4 hypomorphic mutation by RPD3 deletion further supports the conclusion that Rpd3 impedes DDK recruitment by Fkh1, revealing a mechanism of Rpd3 in origin regulation.

Nucleosome-directed replication origin licensing independent of a consensus DNA sequence

The numerous enzymes and cofactors involved in eukaryotic DNA replication are conserved from yeast to human, and the budding yeast Saccharomyces cerevisiae (S.c.) has been a useful model organism for these studies. However, there is a gap in our knowledge of why replication origins in higher eukaryotes do not use a consensus DNA sequence as found in S.c. Using in vitro reconstitution and single-molecule visualization, we show here that S.c. origin recognition complex (ORC) stably binds nucleosomes and that ORC-nucleosome complexes have the intrinsic ability to load the replicative helicase MCM double hexamers onto adjacent nucleosome-free DNA regardless of sequence. Furthermore, we find that Xenopus laevis nucleosomes can substitute for yeast ones in engaging with ORC. Combined with re-analyses of genome-wide ORC binding data, our results lead us to propose that the yeast origin recognition machinery contains the cryptic capacity to bind nucleosomes near a nucleosome-free region and license origins, and that this nucleosome-directed origin licensing paradigm generalizes to all eukaryotes.

Most eukaryotes do not use a consensus DNA sequence as binding sites for the origin recognition complex (ORC) to initiate DNA replication, however budding yeast do. Here the authors show S. cerevisiae ORC can bind nucleosomes near nucleosome-free regions and recruit replicative helicases to form a pre-replication complex independent of the DNA sequence.

DNA repair as a shared hallmark in cancer and ageing

Increasing evidence demonstrates that DNA damage and genome instability play a crucial role in ageing. Mammalian cells have developed a wide range of complex and well‐orchestrated DNA repair pathways to respond to and resolve many different types of DNA lesions that occur from exogenous and endogenous sources. Defects in these repair pathways lead to accelerated or premature ageing syndromes and increase the likelihood of cancer development. Understanding the fundamental mechanisms of DNA repair will help develop novel strategies to treat ageing‐related diseases. Here, we revisit the processes involved in DNA damage repair and how these can contribute to diseases, including ageing and cancer. We also review recent mechanistic insights into DNA repair and discuss how these insights are being used to develop novel therapeutic strategies for treating human disease. We discuss the use of PARP inhibitors in the clinic for the treatment of breast and ovarian cancer and the challenges associated with acquired drug resistance. Finally, we discuss how DNA repair pathway‐targeted therapeutics are moving beyond PARP inhibition in the search for ever more innovative and efficacious cancer therapies.

Characterization of a Novel Heterochromatin Protein 1 Homolog “HP1c” in the Silkworm, Bombyx mori

Simple Summary

Heterochromatin protein 1 (HP1) plays a major role in the formation and maintenance of heterochromatin and in the regulation of gene expression. Five HP1 genes have been found in Drosophila melanogaster and three HP1 genes in Homo sapiens, while in Bombyx mori, two HP1 genes (BmHP1a and BmHP1b) have been reported. In the present study, we analyzed the function of the novel Bombyx mori HP1 gene (BmHP1c), the third HP1 gene in silkworm. BmHP1c has different characteristics from BmHP1a and BmHP1b in terms of transcriptional repression activity, dimer formation, subcellular localization, and effects of RNAi on cell cycle progression. These findings indicate that BmHP1c plays a different role than BmHP1a and BmHP1b.

Abstract

Heterochromatin protein 1 plays an important role in chromatin structure and gene expression regulation. Three HP1 genes have been found in Homo sapiens, and five HP1 genes have been reported in Drosophila melanogaster. On the other hand, in Bombyx mori, only two HP1 genes, BmHP1a and BmHP1b, were reported. In this research, we have reported the molecular and functional characterization of a novel Bombyx mori HP1 gene (BmHP1c), which had stronger transcriptional repression activity than BmHP1a. BmHP1a and BmHP1b is reported to form homo- and heterodimers, but in co-immunoprecipitation experiments, no homo- or hetero-dimer formation of BmHP1c with the other silkworm HP1s is detected. The intracellular localization of BmHP1c is not only in the nucleus but also in the cytoplasm like mammalian HP1γ. In contrast to human HP1a and b, all three BmHP1s were localized preferentially in the regions poorly stained with DAPI. Interestingly, the double knockdown of BmHP1a and b, but not BmHP1c with a or b, arrested the cell cycle at the G2/M phase. These results suggest that BmHP1c is not essential for cell progression and plays a different role than BmHP1a and BmHP1b.

Single-Cell Mononucleotide Microsatellite Analysis Reveals Differential Insertion-Deletion Dynamics in Mouse T Cells

Microsatellite sequences are particularly prone to slippage during DNA replication, forming insertion-deletion loops that, if left unrepaired, result in de novo mutations (expansions or contractions of the repeat array). Mismatch repair (MMR) is a critical DNA repair mechanism that corrects these insertion-deletion loops, thereby maintaining microsatellite stability. MMR deficiency gives rise to the molecular phenotype known as microsatellite instability (MSI). By sequencing MMR-proficient and -deficient (Mlh1+/+ and Mlh1−/−) single-cell exomes from mouse T cells, we reveal here several previously unrecognized features of in vivo MSI. Specifically, mutational dynamics of insertions and deletions were different on multiple levels. Factors that associated with propensity of mononucleotide microsatellites to insertions versus deletions were: microsatellite length, nucleotide composition of the mononucleotide tract, gene length and transcriptional status, as well replication timing. Here, we show on a single-cell level that deletions — the predominant MSI type in MMR-deficient cells — are preferentially associated with longer A/T tracts, long or transcribed genes and later-replicating genes.

Latent Representation of the Human Pan-Celltype Epigenome Through a Deep Recurrent Neural Network

The availability of thousands of assays of epigenetic activity necessitates compressed representations of these data sets that summarize the epigenetic landscape of the genome. Until recently, most such representations were cell type-specific, applying to a single tissue or cell state. Recently, neural networks have made it possible to summarize data across tissues to produce a pan-cell type representation. In this work, we propose Epi-LSTM, a deep long short-term memory (LSTM) recurrent neural network autoencoder to capture the long-term dependencies in the epigenomic data. The latent representations from Epi-LSTM capture a variety of genomic phenomena, including gene-expression, promoter-enhancer interactions, replication timing, frequently interacting regions, and evolutionary conservation. These representations outperform existing methods in a majority of cell types while yielding smoother representations along the genomic axis due to their sequential nature.

Learning representations of chromatin contacts using a recurrent neural network identifies genomic drivers of conformation

Despite the availability of chromatin conformation capture experiments, discerning the relationship between the 1D genome and 3D conformation remains a challenge, which limits our understanding of their affect on gene expression and disease. We propose Hi-C-LSTM, a method that produces low-dimensional latent representations that summarize intra-chromosomal Hi-C contacts via a recurrent long short-term memory neural network model. We find that these representations contain all the information needed to recreate the observed Hi-C matrix with high accuracy, outperforming existing methods. These representations enable the identification of a variety of conformation-defining genomic elements, including nuclear compartments and conformation-related transcription factors. They furthermore enable in-silico perturbation experiments that measure the influence of cis-regulatory elements on conformation.

Dimerization of Firing Factors for Replication Origin Activation in Eukaryotes: A Crucial Process for Simultaneous Assembly of Bidirectional Replication Forks?

Simple Summary

Chromosomal DNA must be faithfully duplicated and segregated into two daughter cells when cells divide. DNA synthesis initiates from specific regions known as the origins of replication. When it starts, a pair of the replication fork is established, and each replication fork moves away from replication origins. In each replication fork, replicative helicase unwinds DNA from a double to a single strand. This implies that two sets of active helicase are generated from each replication origin. To make this possible, two sets of replicative helicases are loaded onto replication origins as inactive dimers first. When S-phase specific cyclin-dependent kinases, S-CDKs, are activated, the inactive helicase is converted into the active helicase with the aid of other factors called firing factors. Although two sets of firing factors seem to be required to activate two sets of helicase, it is largely unknown whether two sets of firing factors function simultaneously to establish bidirectional replication forks in a coordinated way. We introduce our current understanding of firing factor dimerization and discuss its potential contribution to bidirectional replication fork formation in this review.

Abstract

Controlling the activity of the heterohexameric Mcm2–7 replicative helicase is crucial for regulation of replication origin activity in eukaryotes. Because bidirectional replication forks are generated from every replication origin, when origins are licensed for replication in the first step of DNA replication, two inactive Mcm2–7 heterohexiameric complexes are loaded around double stranded DNA as a head-to-head double hexamer. The helicases are subsequently activated via a ‘firing’ reaction, in which the Mcm2–7 double hexamer is converted into two active helicase units, the CMG complex, by firing factors. Dimerization of firing factors may contribute to this process by allowing simultaneous activation of two sets of helicases and thus efficient assembly of bidirectional replication forks. An example of this is dimerization of the firing factor Sld3/Treslin/Ticrr via its binding partner, Sld7/MTBP. In organisms in which no Sld7 ortholog has been identified, such as the fission yeast Schizosaccharomyces pombe, Sld3 itself has a dimerization domain, and it has been suggested that this self-interaction is crucial for the firing reaction in this organism. Dimerization induces a conformational change in Sdl3 that appears to be critical for the firing reaction. Moreover, Mcm10 also seems to be regulated by self-interaction in yeasts. Although it is not yet clear to what extent dimerization of firing factors contributes to the firing reaction in eukaryotes, we discuss the possible roles of firing factor dimerization in simultaneous helicase activation.

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.

The Role of MTBP as a Replication Origin Firing Factor

The initiation step of replication at replication origins determines when and where in the genome replication machines, replisomes, are generated. Tight control of replication initiation helps facilitate the two main tasks of genome replication, to duplicate the genome accurately and exactly once each cell division cycle. The regulation of replication initiation must ensure that initiation occurs during the S phase specifically, that no origin fires more than once per cell cycle, that enough origins fire to avoid non-replicated gaps, and that the right origins fire at the right time but only in favorable circumstances. Despite its importance for genetic homeostasis only the main molecular processes of eukaryotic replication initiation and its cellular regulation are understood. The MTBP protein (Mdm2-binding protein) is so far the last core replication initiation factor identified in metazoan cells. MTBP is the orthologue of yeast Sld7. It is essential for origin firing, the maturation of pre-replicative complexes (pre-RCs) into replisomes, and is emerging as a regulation focus targeted by kinases and by regulated degradation. We present recent insight into the structure and cellular function of the MTBP protein in light of recent structural and biochemical studies revealing critical molecular details of the eukaryotic origin firing reaction. How the roles of MTBP in replication and other cellular processes are mutually connected and are related to MTBP's contribution to tumorigenesis remains largely unclear.

Differential condensation of sister chromatids acts with Cdc6 to ensure asynchronous S-phase entry in Drosophila male germline stem cell lineage

During Drosophila melanogaster male germline stem cell (GSC) asymmetric division, preexisting old versus newly synthesized histones H3 and H4 are asymmetrically inherited. However, the biological outcomes of this phenomenon have remained unclear. Here, we tracked old and new histones throughout the GSC cell cycle through the use of high spatial and temporal resolution microscopy. We found unique features that differ between old and new histone-enriched sister chromatids, including differences in nucleosome density, chromosomal condensation, and H3 Ser10 phosphorylation. These distinct chromosomal features lead to their differential association with Cdc6, a pre-replication complex component, and subsequent asynchronous DNA replication initiation in the resulting daughter cells. Disruption of asymmetric histone inheritance abolishes differential Cdc6 association and asynchronous S-phase entry, demonstrating that histone asymmetry acts upstream of these critical cell-cycle progression events. Furthermore, disruption of these GSC-specific chromatin features leads to GSC defects, indicating a connection between histone inheritance, cell-cycle progression, and cell fate determination.

S Phase Duration Is Determined by Local Rate and Global Organization of Replication

Simple Summary

In order for a cell to divide into two cells, it must first copy its DNA. Although the time required for this process tends not to vary much, many examples of the importance of variability have been reported. In this review, we discuss the methods used to study this question, present some of the examples of variation, and attempt to explain the factors that determine the time required in simple terms. We will show that the overall time depends on the rate of DNA replication within a region, and on the temporal organization of the regions relative to each other.

Abstract

The duration of the cell cycle has been extensively studied and a wide degree of variability exists between cells, tissues and organisms. However, the duration of S phase has often been neglected, due to the false assumption that S phase duration is relatively constant. In this paper, we describe the methodologies to measure S phase duration, summarize the existing knowledge about its variability and discuss the key factors that control it. The local rate of replication (LRR), which is a combination of fork rate (FR) and inter-origin distance (IOD), has a limited influence on S phase duration, partially due to the compensation between FR and IOD. On the other hand, the organization of the replication program, specifically the amount of replication domains that fire simultaneously and the degree of overlap between the firing of distinct replication timing domains, is the main determinant of S phase duration. We use these principles to explain the variation in S phase length in different tissues and conditions.

Inheritance of repressed chromatin domains during S phase requires the histone chaperone NPM1

The epigenetic process safeguards cell identity during cell division through the inheritance of appropriate gene expression profiles. We demonstrated previously that parental nucleosomes are inherited by the same chromatin domains during DNA replication only in the case of repressed chromatin. We now show that this specificity is conveyed by NPM1, a histone H3/H4 chaperone. Proteomic analyses of late S-phase chromatin revealed NPM1 in association with both H3K27me3, an integral component of facultative heterochromatin, and MCM2, an integral component of the DNA replication machinery; moreover, NPM1 interacts directly with PRC2 and with MCM2. Given that NPM1 is essential, the inheritance of repressed chromatin domains was examined anew using mESCs expressing an auxin-degradable version of endogenous NPM1. Upon NPM1 degradation, cells accumulated in the G₁-S phase of the cell cycle and parental nucleosome inheritance from repressed chromatin domains was markedly compromised. NPM1 chaperone activity may contribute to the integrity of this process as appropriate inheritance required the NPM1 acidic patches.

In Silico Analysis Identified Putative Pathogenic Missense nsSNPs in Human SLITRK1 Gene

Human DNA contains several variations, which can affect the structure and normal functioning of a protein. These variations could be single nucleotide polymorphisms (SNPs) or insertion-deletions (InDels). SNPs, as opposed to InDels, are more commonly present in DNA and may cause genetic disorders. In the current study, several bioinformatic tools were used to prioritize the pathogenic variants in the SLITRK1 gene. Out of all of the variants, 16 were commonly predicted to be pathogenic by these tools. All the variants had very low frequency, i.e., <0.0001 in the global population. The secondary structure of all filtered variants was predicted, but no structural change was observed at the site of variation in any variant. Protein stability analysis of these variants was then performed, which determined a decrease in protein stability of 10 of the variants. Amino acid conservation analysis revealed that all the amino acids were highly conserved, indicating their structural and functional importance. Protein 3D structure of wildtype SLITRK1 and all of its variants was predicted using I-TASSER, and the effect of variation on 3D structure of the protein was observed using the Missense3D tool, which presented the probable structural loss in three variants, i.e., Asn529Lys, Leu496Pro and Leu94Phe. The wildtype SLITRK1 protein and these three variants were independently docked with their close interactor protein PTPRD, and remarkable differences were observed in the docking sites of normal and variants, which will ultimately affect the functional activity of the SLITRK1 protein. Previous studies have shown that mutations in SLITRK1 are involved in Tourette syndrome. The present study may assist a molecular geneticist in interpreting the variant pathogenicity in research as well as diagnostic setup.

Rif1-Dependent Control of Replication Timing

Successful duplication of the genome requires the accurate replication of billions of base pairs of DNA within a relatively short time frame. Failure to accurately replicate the genome results in genomic instability and a host of diseases. To faithfully and rapidly replicate the genome, DNA replication must be tightly regulated and coordinated with many other nuclear processes. These regulations, however, must also be flexible as replication kinetics can change through development and differentiation. Exactly how DNA replication is regulated and how this regulation changes through development is an active field of research. One aspect of genome duplication where much remains to be discovered is replication timing (RT), which dictates when each segment of the genome is replicated during S phase. All organisms display some level of RT, yet the precise mechanisms that govern RT remain are not fully understood. The study of Rif1, a protein that actively regulates RT from yeast to humans, provides a key to unlock the underlying molecular mechanisms controlling RT. The paradigm for Rif1 function is to delay helicase activation within certain regions of the genome, causing these regions to replicate late in S phase. Many questions, however, remain about the intricacies of Rif1 function. Here, we review the current models for the activity of Rif1 with the goal of trying to understand how Rif1 functions to establish the RT program.

Chemically Induced Chromosomal Interaction (CICI) method to study chromosome dynamics and its biological roles

Numerous intra- and inter-chromosomal contacts have been mapped in eukaryotic genomes, but it remains challenging to link these 3D structures to their regulatory functions. To establish the causal relationships between chromosome conformation and genome functions, we  develop a method, Chemically Induced Chromosomal Interaction (CICI), to selectively perturb the chromosome conformation at targeted loci. Using this method, long-distance chromosomal interactions can be induced dynamically between intra- or inter-chromosomal loci pairs, including the ones with very low Hi-C contact frequencies. Measurement of CICI formation time allows us to probe chromosome encounter dynamics between different loci pairs across the cell cycle. We also conduct two functional tests of CICI. We perturb the chromosome conformation near a DNA double-strand break and observe altered donor preference in homologous recombination; we force interactions between early and late-firing DNA replication origins and find no significant changes in replication timing. These results suggest that chromosome conformation plays a deterministic role in homology-directed DNA repair, but not in the establishment of replication timing. Overall, our study demonstrates that CICI is a powerful tool to study chromosome dynamics and 3D genome function.

Methods to selectively manipulate specific long-distance chromosomal interactions are limited. Here the authors develop a method called Chemically Induced Chromosomal Interaction (CICI) to engineer interactions and demonstrate that 3D conformation plays a causal role in establishing donor DNA preference during DNA repair.

f = m*a: A Framework for Investigating the Regulation of Replication Timing

Stochastic models of replication timing posit that origin firing timing is regulated by origin firing probability, with early-firing origins having a high probability of firing and late-firing origins having a lower probability. However, they offer no insight into why one origin should have a higher firing probability than another. Here, a simple framework is suggested for how to approach the question by noting that the firing probability (f) must be the product of the stoichiometry of the MCM replicative helicase loaded at the origin (m) and the probability with which that MCM is activated (a). This framework emphasizes that mechanistic understanding of replication timing must focus on MCM loading and activation and can be simplified to the equation f = m*a.

The consequences of differential origin licensing dynamics in distinct chromatin environments

Eukaryotic chromosomes contain regions of varying accessibility, yet DNA replication factors must access all regions. The first replication step is loading MCM complexes to license replication origins during the G1 cell cycle phase. It is not yet known how mammalian MCM complexes are adequately distributed to both accessible euchromatin regions and less accessible heterochromatin regions. To address this question, we combined time-lapse live-cell imaging with immunofluorescence imaging of single human cells to quantify the relative rates of MCM loading in euchromatin and heterochromatin throughout G1. We report here that MCM loading in euchromatin is faster than that in heterochromatin in early G1, but surprisingly, heterochromatin loading accelerates relative to euchromatin loading in middle and late G1. This differential acceleration allows both chromatin types to begin S phase with similar concentrations of loaded MCM. The different loading dynamics require ORCA-dependent differences in origin recognition complex distribution. A consequence of heterochromatin licensing dynamics is that cells experiencing a truncated G1 phase from premature cyclin E expression enter S phase with underlicensed heterochromatin, and DNA damage accumulates preferentially in heterochromatin in the subsequent S/G2 phase. Thus, G1 length is critical for sufficient MCM loading, particularly in heterochromatin, to ensure complete genome duplication and to maintain genome stability.

CCPE: cell cycle pseudotime estimation for single cell RNA-seq data

Pseudotime analysis from scRNA-seq data enables to characterize the continuous progression of various biological processes, such as the cell cycle. Cell cycle plays an important role in cell fate decisions and differentiation and is often regarded as a confounder in scRNA-seq data analysis when analyzing the role of other factors. Therefore, accurate prediction of cell cycle pseudotime and identification of cell cycle stages are important steps for characterizing the development-related biological processes. Here, we develop CCPE, a novel cell cycle pseudotime estimation method to characterize cell cycle timing and identify cell cycle phases from scRNA-seq data. CCPE uses a discriminative helix to characterize the circular process of the cell cycle and estimates each cell's pseudotime along the cell cycle. We evaluated the performance of CCPE based on a variety of simulated and real scRNA-seq datasets. Our results indicate that CCPE is an effective method for cell cycle estimation and competitive in various applications compared with other existing methods. CCPE successfully identified cell cycle marker genes and is robust to dropout events in scRNA-seq data. Accurate prediction of the cell cycle using CCPE can also effectively facilitate the removal of cell cycle effects across cell types or conditions.