Conservation of replication timing reveals global and local regulation of replication origin activity

Integrating high-throughput analysis to create an atlas of replication origins in Trypanosoma cruzi in the context of genome structure and variability

Trypanosoma cruzi is the etiologic agent of the most prevalent human parasitic disease in Latin America, Chagas disease. Its genome is rich in multigenic families that code for virulent antigens and are present in the rapidly evolving genomic compartment named Disruptive. DNA replication is a meticulous biological process in which flaws can generate mutations and changes in chromosomal and gene copy numbers. Here, integrating high-throughput and single-molecule analyses, we were able to identify Predominant, Flexible, and Dormant Orc1Cdc6-dependent origins as well as Orc1Cdc6-independent origins. Orc1Cdc6-dependent origins were found in multigenic family loci, while independent origins were found in the Core compartment that contains conserved and hypothetical protein-coding genes, in addition to multigenic families. In addition, we found that Orc1Cdc6 density is related to the firing of origins and that Orc1Cdc6-binding sites within fired origins are depleted of a specific class of nucleosomes that we previously categorized as dynamic. Together, these data suggest that Orc1Cdc6-dependent origins may contribute to the rapid evolution of the Disruptive compartment and, therefore, to the success of T. cruzi infection and that the local epigenome landscape is also involved in this process.IMPORTANCETrypanosoma cruzi, responsible for Chagas disease, affects millions globally, particularly in Latin America. Lack of vaccine or treatment underscores the need for research. Parasite's genome, with virulent antigen-coding multigenic families, resides in the rapidly evolving Disruptive compartment. Study sheds light on the parasite's dynamic DNA replication, discussing the evolution of the Disruptive compartment. Therefore, the findings represent a significant stride in comprehending T. cruzi's biology and the molecular bases that contribute to the success of infection caused by this parasite.

Design, construction, and functional characterization of a tRNA neochromosome in yeast

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Cohesin maintains replication timing to suppress DNA damage on cancer genes

Cohesin loss-of-function mutations are frequently observed in tumors, but the mechanism underlying its role in tumorigenesis is unclear. Here, we found that depletion of RAD21, a core subunit of cohesin, leads to massive genome-wide DNA breaks and 147 translocation hotspot genes, co-mutated with cohesin in multiple cancers. Increased DNA damages are independent of RAD21-loss-induced transcription alteration and loop anchor disruption. However, damage-induced chromosomal translocations coincide with the asymmetrically distributed Okazaki fragments of DNA replication, suggesting that RAD21 depletion causes replication stresses evidenced by the slower replication speed and increased stalled forks. Mechanistically, approximately 30% of the human genome exhibits an earlier replication timing after RAD21 depletion, caused by the early initiation of >900 extra dormant origins. Correspondingly, most translocation hotspot genes lie in timing-altered regions. Therefore, we conclude that cohesin dysfunction causes replication stresses induced by excessive DNA replication initiation, resulting in gross DNA damages that may promote tumorigenesis.

Contributions of replicative and translesion DNA polymerases to mutagenic bypass of canonical and atypical UV photoproducts

UV exposure induces a mutation signature of C > T substitutions at dipyrimidines in skin cancers. We recently identified additional UV-induced AC > TT and A > T substitutions that could respectively cause BRAF V600K and V600E oncogenic mutations. The mutagenic bypass mechanism past these atypical lesions, however, is unknown. Here, we whole genome sequenced UV-irradiated yeast and used reversion reporters to delineate the roles of replicative and translesion DNA polymerases in mutagenic bypass of UV-lesions. Our data indicates that yeast DNA polymerase eta (pol η) has varied impact on UV-induced mutations: protecting against C > T substitutions, promoting T > C and AC > TT substitutions, and not impacting A > T substitutions. Surprisingly, deletion rad30Δ increased novel UV-induced C > A substitutions at CA dinucleotides. In contrast, DNA polymerases zeta (pol ζ) and epsilon (pol ε) participated in AC > TT and A > T mutations. These results uncover lesion-specific accurate and mutagenic bypass of UV lesions, which likely contribute to key driver mutations in melanoma.

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.

Construction of a novel plasmid for an industrial yeast Candida glycerinogenes by dual-autonomously replicating sequence strategy

Due to the lack of available episomal plasmid, the improvement of many industrial strains, especially exogenous gene expression, is severely restricted. The failure of autonomous replication or low copy number of episomal plasmids is the main reason for the failure of many episomal plasmids construction. In this paper, Candida glycerinogenes, an industrial strain lacking episomal plasmids, was employed as the topic. A series of GFP-based plasmids containing autonomously replicating sequence (ARS) from different strain sources were constructed and analyzed for performance, and it was found that only the panARS from Kluyveromyces lactis compared with other nine low capacity ARSs proved to have the best performance and could be used to construct episomal plasmid. Further, the dual-ARS strategy was used to optimize the episomal plasmid, and the results indicated that only the dual-ARS plasmid +PPARS2 with double different ARSs, not the dual-ARS plasmid +panARS with double same ARSs, showed an improvement in all properties, with an increase in transformation efficiency of about 36% and a synchronous trend of fluorescence intensity and copy number, both by about 40%. In addition, constructed episomal plasmids were used to express the exogenous gene CrGES, and the fact that geraniol was found proved the versatility of the plasmids. The successful construction of episomal plasmids will also substantially facilitate genetic engineering research and industrial use of C. glycerinogenes in the future, as well as providing a feasible approach to create episomal plasmids for industrial strains.

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.

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.

Is the mutation rate lower in genomic regions of stronger selective constraints?

A study of the plant Arabidopsis thaliana detected lower mutation rates in genomic regions where mutations are more likely to be deleterious, challenging the principle that mutagenesis is blind to its consequence. To examine the generality of this finding, we analyze large mutational data from baker's yeast and humans. The yeast data do not exhibit this trend, whereas the human data show an opposite trend that disappears upon the control of potential confounders. We find that the Arabidopsis study identified substantially more mutations than reported in the original data-generating studies and expected from Arabidopsis' mutation rate. These extra mutations are enriched in polynucleotide tracts and have relatively low sequencing qualities so are likely sequencing errors. Furthermore, the polynucleotide “mutations” can produce the purported mutational trend in Arabidopsis. Together, our results do not support lower mutagenesis of genomic regions of stronger selective constraints in the plant, fungal, and animal models examined.

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.

Genomic analysis finds no evidence of canonical eukaryotic DNA processing complexes in a free-living protist

Cells replicate and segregate their DNA with precision. Previous studies showed that these regulated cell-cycle processes were present in the last eukaryotic common ancestor and that their core molecular parts are conserved across eukaryotes. However, some metamonad parasites have secondarily lost components of the DNA processing and segregation apparatuses. To clarify the evolutionary history of these systems in these unusual eukaryotes, we generated a genome assembly for the free-living metamonad Carpediemonas membranifera and carried out a comparative genomics analysis. Here, we show that parasitic and free-living metamonads harbor an incomplete set of proteins for processing and segregating DNA. Unexpectedly, Carpediemonas species are further streamlined, lacking the origin recognition complex, Cdc6 and most structural kinetochore subunits. Carpediemonas species are thus the first known eukaryotes that appear to lack this suite of conserved complexes, suggesting that they likely rely on yet-to-be-discovered or alternative mechanisms to carry out these fundamental processes.

Conservation of a DNA Replication Motif among Phylogenetically Distant Budding Yeast Species

Eukaryotic DNA replication begins at genomic loci termed origins, which are bound by the origin recognition complex (ORC). Although ORC is conserved across species, the sequence composition of origins is more varied. In the budding yeast Saccharomyces cerevisiae, the ORC-binding motif consists of an A/T-rich 17 bp “extended ACS” sequence adjacent to a B1 element composed of two 3-bp motifs. Similar sequences occur at origins in closely related species, but it is not clear when this type of replication origin arose and whether it predated a whole-genome duplication that occurred around 100 Ma in the budding yeast lineage. To address these questions, we identified the ORC-binding sequences in the nonduplicated species Torulaspora delbrueckii. We used chromatin immunoprecipitation followed by sequencing and identified 190 ORC-binding sites distributed across the eight T. delbrueckii chromosomes. Using these sites, we identified an ORC-binding motif that is nearly identical to the known motif in S. cerevisiae. We also found that the T. delbrueckii ORC-binding sites function as origins in T. delbrueckii when cloned onto a plasmid and that the motif is required for plasmid replication. Finally, we compared an S. cerevisiae origin with two T. delbrueckii ORC-binding sites and found that they conferred similar stabilities to a plasmid. These results reveal that the ORC-binding motif arose prior to the whole-genome duplication and has been maintained for over 100 Myr.

An evolutionary model identifies the main evolutionary biases for the evolution of genome-replication profiles

Recent results comparing the temporal program of genome replication of yeast species belonging to the Lachancea clade support the scenario that the evolution of the replication timing program could be mainly driven by correlated acquisition and loss events of active replication origins. Using these results as a benchmark, we develop an evolutionary model defined as birth-death process for replication origins and use it to identify the evolutionary biases that shape the replication timing profiles. Comparing different evolutionary models with data, we find that replication origin birth and death events are mainly driven by two evolutionary pressures, the first imposes that events leading to higher double-stall probability of replication forks are penalized, while the second makes less efficient origins more prone to evolutionary loss. This analysis provides an empirically grounded predictive framework for quantitative evolutionary studies of the replication timing program.

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

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

The capacity of origins to load MCM establishes replication timing patterns

Loading of the MCM replicative helicase at origins of replication is a highly regulated process that precedes DNA replication in all eukaryotes. The stoichiometry of MCM loaded at origins has been proposed to be a key determinant of when those origins initiate replication during S phase. Nevertheless, the genome-wide regulation of MCM loading stoichiometry and its direct effect on replication timing remain unclear. In order to investigate why some origins load more MCM than others, we perturbed MCM levels in budding yeast cells and, for the first time, directly measured MCM levels and replication timing in the same experiment. Reduction of MCM levels through degradation of Mcm4, one of the six obligate components of the MCM complex, slowed progression through S phase and increased sensitivity to replication stress. Reduction of MCM levels also led to differential loading at origins during G1, revealing origins that are sensitive to reductions in MCM and others that are not. Sensitive origins loaded less MCM under normal conditions and correlated with a weak ability to recruit the origin recognition complex (ORC). Moreover, reduction of MCM loading at specific origins of replication led to a delay in their replication during S phase. In contrast, overexpression of MCM had no effects on cell cycle progression, relative MCM levels at origins, or replication timing, suggesting that, under optimal growth conditions, cellular MCM levels are not limiting for MCM loading. Our results support a model in which the loading capacity of origins is the primary determinant of MCM stoichiometry in wild-type cells, but that stoichiometry is controlled by origins' ability to recruit ORC and compete for MCM when MCM becomes limiting.

Prospect of reprogramming replication licensing for cancer drug development

Eukaryotic chromosomal DNA replication is preceded by replication licensing which involves the identification of the origin of replication by origin recognition complex (ORC). The ORC loads pre-replication complexes (pre-RCs) through a series of tightly regulated mechanisms where the ORC interacts with Cdc6 to recruit cdt1-MCM2-7 complexes to the origin of replication. In eukaryotes, adherence to regulatory mechanisms of the replication program is required to ensure that all daughter cells carry the exact copy of genetic material as the parent cell. Failure of which leads to the development of genome instability syndromes like cancer, diabetes, etc. In an event of such occurrence, preventing cells from carrying the defaulted genetic material and passing it to other cells hinges on the regulation of chromosomal DNA replication. Thus, understanding the mechanisms underpinning chromosomal DNA replication and particularly replication licensing can expose druggable enzymes, effector molecules, and secondary messengers that can be targeted for diagnosis and therapeutic purposes. Effectively drugging these molecular markers to reprogram pre-replication events can be used to control the fate of chromosomal DNA replication for the treatment of genome instability disorders and in this case, cancer. This review discusses available knowledge of replication licensing in the contest of molecular drug discovery for the treatment of cancer.

Tos4 mediates gene expression homeostasis through interaction with HDAC complexes independently of H3K56 acetylation

Saccharomyces cerevisiae exhibits gene expression homeostasis, which is defined as the buffering of transcription levels against changes in DNA copy number during the S phase of the cell cycle. It has been suggested that S. cerevisiae employs an active mechanism to maintain gene expression homeostasis through Rtt109-Asf1-dependent acetylation of histone H3 on lysine 56 (H3K56). Here, we show that gene expression homeostasis can be achieved independently of H3K56 acetylation by Tos4 (Target of Swi6-4). Using Nanostring technology, we establish that Tos4-dependent gene expression homeostasis depends on its forkhead-associated (FHA) domain, which is a phosphopeptide recognition domain required to bind histone deacetylases (HDACs). We demonstrate that the mechanism of Tos4-dependent gene expression homeostasis requires its interaction with the Rpd3L HDAC complex. However, this is independent of Rpd3's well-established roles in both histone deacetylation and controlling the DNA replication timing program, as established by deep sequencing of Fluorescence-Activated Cell Sorted (FACS) S and G2 phase populations. Overall, our data reveals that Tos4 mediates gene expression homeostasis through its FHA domain-dependent interaction with the Rpd3L complex, which is independent of H3K56ac.

Comparing DNA replication programs reveals large timing shifts at centromeres of endocycling cells in maize roots

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed.

In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

Read, Write, Adapt: Challenges and Opportunities during Kinetoplastid Genome Replication

The genomes of all organisms are read throughout their growth and development, generating new copies during cell division and encoding the cellular activities dictated by the genome’s content. However, genomes are not invariant information stores but are purposefully altered in minor and major ways, adapting cellular behaviour and driving evolution. Kinetoplastids are eukaryotic microbes that display a wide range of such read–write genome activities, in many cases affecting critical aspects of their biology, such as host adaptation. Here we discuss the range of read–write genome changes found in two well-studied kinetoplastid parasites, Trypanosoma brucei and Leishmania, focusing on recent work that suggests such adaptive genome variation is linked to novel strategies the parasites use to replicate their unconventional genomes.

Polycistronic transcription dominates and shapes kinetoplastid genomes, inevitably leading to clashes with DNA replication. By harnessing the resultant DNA damage for adaptation, kinetoplastids have huge potential for dynamic read–write genome variation.

Major origins of DNA replication are confined to the boundaries of polycistronic transcription units in the Trypanosoma brucei and Leishmania genomes, putatively limiting DNA damage. Subtelomeres may lack this arrangement, generating read–write hotspots.

In T. brucei, early replication of the highly transcribed subtelomeric variant surface glycoprotein (VSG) expression site may ensure replication-transcription clashes within this site to trigger DNA recombination, an event critical for antigenic variation.

Leishmania genomes show extensive aneuploidy and copy number variation. Notably, DNA replication requires recombination factors and relies on post-S phase replication of subtelomeres.

Evolution of compartmentalised DNA replication programmes underpin important aspects of genome biology in kinetoplastids, illustrating the consolidation of genome maintenance strategies to promote genome plasticity.

Genome duplication in Leishmania major relies on persistent subtelomeric DNA replication

DNA replication is needed to duplicate a cell’s genome in S phase and segregate it during cell division. Previous work in Leishmania detected DNA replication initiation at just a single region in each chromosome, an organisation predicted to be insufficient for complete genome duplication within S phase. Here, we show that acetylated histone H3 (AcH3), base J and a kinetochore factor co-localise in each chromosome at only a single locus, which corresponds with previously mapped DNA replication initiation regions and is demarcated by localised G/T skew and G4 patterns. In addition, we describe previously undetected subtelomeric DNA replication in G2/M and G1-phase-enriched cells. Finally, we show that subtelomeric DNA replication, unlike chromosome-internal DNA replication, is sensitive to hydroxyurea and dependent on 9-1-1 activity. These findings indicate that Leishmania’s genome duplication programme employs subtelomeric DNA replication initiation, possibly extending beyond S phase, to support predominantly chromosome-internal DNA replication initiation within S phase.

Isolation and analysis of rereplicated DNA by Rerep-Seq

Changes in gene copy number contribute to genomic instability, the onset and progression of cancer, developmental abnormalities and adaptive potential. The origins of gene amplifications have remained elusive; however, DNA rereplication has been implicated as a source of gene amplifications. The inability to determine which sequences are rereplicated and under what conditions have made it difficult to determine the validity of the proposed models. Here we present Rerep-Seq, a technique that selectively enriches for rereplicated DNA in preparation for analysis by DNA sequencing that can be applied to any species. We validated Rerep-Seq by simulating DNA rereplication in yeast and human cells. Using Rerep-Seq, we demonstrate that rereplication induced in Saccharomyces cerevisiae by deregulated origin licensing is non-random and defined by broad domains that span multiple replication origins and topological boundaries.

Conditional knockout of RAD51-related genes in Leishmania major reveals a critical role for homologous recombination during genome replication

Homologous recombination (HR) has an intimate relationship with genome replication, both during repair of DNA lesions that might prevent DNA synthesis and in tackling stalls to the replication fork. Recent studies led us to ask if HR might have a more central role in replicating the genome of Leishmania, a eukaryotic parasite. Conflicting evidence has emerged regarding whether or not HR genes are essential, and genome-wide mapping has provided evidence for an unorthodox organisation of DNA replication initiation sites, termed origins. To answer this question, we have employed a combined CRISPR/Cas9 and DiCre approach to rapidly generate and assess the effect of conditional ablation of RAD51 and three RAD51-related proteins in Leishmania major. Using this approach, we demonstrate that loss of any of these HR factors is not immediately lethal but in each case growth slows with time and leads to DNA damage and accumulation of cells with aberrant DNA content. Despite these similarities, we show that only loss of RAD51 or RAD51-3 impairs DNA synthesis and causes elevated levels of genome-wide mutation. Furthermore, we show that these two HR factors act in distinct ways, since ablation of RAD51, but not RAD51-3, has a profound effect on DNA replication, causing loss of initiation at the major origins and increased DNA synthesis at subtelomeres. Our work clarifies questions regarding the importance of HR to survival of Leishmania and reveals an unanticipated, central role for RAD51 in the programme of genome replication in a microbial eukaryote.

Replication origin location might contribute to genetic variability in Trypanosoma cruzi

Background

DNA replication in trypanosomatids operates in a uniquely challenging environment, since most of their genomes are constitutively transcribed. Trypanosoma cruzi, the etiological agent of Chagas disease, presents high variability in both chromosomes size and copy number among strains, though the underlying mechanisms are unknown.

Results

Here we have mapped sites of DNA replication initiation across the T. cruzi genome using Marker Frequency Analysis, which has previously only been deployed in two related trypanosomatids. The putative origins identified in T. cruzi show a notable enrichment of GC content, a preferential position at subtelomeric regions, coinciding with genes transcribed towards the telomeres, and a pronounced enrichment within coding DNA sequences, most notably in genes from the Dispersed Gene Family 1 (DGF-1).

Conclusions

These findings suggest a scenario where collisions between DNA replication and transcription are frequent, leading to increased genetic variability, as seen by the increase SNP levels at chromosome subtelomeres and in DGF-1 genes containing putative origins.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

DNA copy-number measurement of genome replication dynamics by high-throughput sequencing: the sort-seq, sync-seq and MFA-seq family

Genome replication follows a defined temporal programme that can change during cellular differentiation and disease onset. DNA replication results in an increase in DNA copy number that can be measured by high-throughput sequencing. Here we present a protocol to determine genome replication dynamics using DNA copy-number measurements. Cell populations can be obtained in three variants of the method. First, sort-seq reveals the average replication dynamics across S phase in an unperturbed cell population; FACS is used to isolate replicating and non-replicating subpopulations from asynchronous cells. Second, sync-seq measures absolute replication time at specific points during S phase using a synchronized cell population. Third, marker frequency analysis can be used to reveal the average replication dynamics using copy-number analysis in any proliferating asynchronous cell culture. These approaches have been used to reveal genome replication dynamics in prokaryotes, archaea and a wide range of eukaryotes, including yeasts and mammalian cells. We have found this approach straightforward to apply to other organisms and highlight example studies from across the three domains of life. Here we present a Saccharomyces cerevisiae version of the protocol that can be performed in 7-10 d. It requires basic molecular and cellular biology skills, as well as a basic understanding of Unix and R.

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

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

Linking the organization of DNA replication with genome maintenance

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

Rapid high-resolution measurement of DNA replication timing by droplet digital PCR

Genomes are replicated in a reproducible temporal pattern. Current methods for assaying allele replication timing are time consuming and/or expensive. These include high-throughput sequencing which can be used to measure DNA copy number as a proxy for allele replication timing. Here, we use droplet digital PCR to study DNA replication timing at multiple loci in budding yeast and human cells. We establish that the method has temporal and spatial resolutions comparable to the high-throughput sequencing approaches, while being faster than alternative locus-specific methods. Furthermore, the approach is capable of allele discrimination. We apply this method to determine relative replication timing across timing transition zones in cultured human cells. Finally, multiple samples can be analysed in parallel, allowing us to rapidly screen kinetochore mutants for perturbation to centromere replication timing. Therefore, this approach is well suited to the study of locus-specific replication and the screening of cis- and trans-acting mutants to identify mechanisms that regulate local genome replication timing.

Genomic Analysis of DNA Double-Strand Break Repair in Escherichia coli

Counting DNA whole genome sequencing reads is providing new insight into DNA double-strand break repair (DSBR) in the model organism Escherichia coli. We describe the application of RecA chromatin immunoprecipitation coupled to genomic DNA sequencing (RecA-ChIP-seq) and marker frequency analysis (MFA) to analyze the genomic consequences of DSBR. We provide detailed procedures for the preparation of DNA and the analysis of data. We compare different ways of visualizing ChIP data and show that alternative protocols for the preparation of DNA for MFA differentially affect the recovery of branched DNA molecules containing Holliday junctions.

The organization of genome duplication is a critical determinant of the landscape of genome maintenance

Genome duplication is essential for cell proliferation, and the mechanisms regulating its execution are highly conserved. These processes give rise to a spatiotemporal organization of replication initiation across the genome, referred to as the replication program. Despite the identification of such programs in diverse eukaryotic organisms, their biological importance for cellular physiology remains largely unexplored. We address this fundamental question in the context of genome maintenance, taking advantage of the inappropriate origin firing that occurs when fission yeast cells lacking the Rad3/ATR checkpoint kinase are subjected to replication stress. Using this model, we demonstrate that the replication program quantitatively dictates the extent of origin de-regulation and the clustered localization of these events. Furthermore, our results uncover an accumulation of abnormal levels of single-stranded DNA (ssDNA) and the Rad52 repair protein at de-regulated origins. We show that these loci constitute a defining source of the overall ssDNA and Rad52 hotspots in the genome, generating a signature pattern of instability along the chromosomes. We then induce a genome-wide reprogramming of origin usage and evaluate its consequences in our experimental system. This leads to a complete redistribution of the sites of both inappropriate initiation and associated Rad52 recruitment. We therefore conclude that the organization of genome duplication governs the checkpoint control of origin-associated hotspots of instability and plays an integral role in shaping the landscape of genome maintenance.

Regulation of the program of DNA replication by CDK: new findings and perspectives

Progression through the cell cycle is driven by the activities of the cyclin-dependent kinase (CDK) family of enzymes, which establish an ordered passage through the cell cycle phases. CDK activity is crucial for the cellular transitions from G1 to S and G2 to M, which are highly controlled to promote the faithful duplication of the genetic material and the transmission of the genome into daughter cells, respectively. While oscillations in CDK activity are essential for cell division, how its specific dynamics may shape cellular processes remains an open question. Recently, we have investigated the potential role of CDK in establishing the profile of replication initiation along the chromosomes, also referred to as the replication program. Our results demonstrated that the timing and level of CDK activity at G1/S provide two critical and independent inputs that modulate the pattern of origin usage. In this review, we will present the conclusions of our study and discuss the implications of our findings for cellular function and physiology.

The evolution of the temporal program of genome replication

Genome replication is highly regulated in time and space, but the rules governing the remodeling of these programs during evolution remain largely unknown. We generated genome-wide replication timing profiles for ten Lachancea yeasts, covering a continuous evolutionary range from closely related to more divergent species. We show that replication programs primarily evolve through a highly dynamic evolutionary renewal of the cohort of active replication origins. We found that gained origins appear with low activity yet become more efficient and fire earlier as they evolutionarily age. By contrast, origins that are lost comprise the complete range of firing strength. Additionally, they preferentially occur in close vicinity to strong origins. Interestingly, despite high evolutionary turnover, active replication origins remain regularly spaced along chromosomes in all species, suggesting that origin distribution is optimized to limit large inter-origin intervals. We propose a model on the evolutionary birth, death, and conservation of active replication origins.

Temporal programs of genome replication show different levels of conservation between closely or distantly related species. Here, the authors generate genome-wide replication timing profiles for ten yeast species, and analyze their evolutionary dynamics.

G2-seq: A High Throughput Sequencing-based Technique for Identifying Late Replicating Regions of the Genome

Numerous techniques have been developed to follow the progress of DNA replication through the S phase of the cell cycle. Most of these techniques have been directed toward elucidation of the location and timing of initiation of genome duplication rather than its completion. However, it is critical that we understand regions of the genome that are last to complete replication, because these regions suffer elevated levels of chromosomal breakage and mutation, and they have been associated with both disease and aging. Here we describe how we have extended a technique that has been used to monitor replication initiation to instead identify those regions of the genome last to complete replication. This approach is based on a combination of flow cytometry and high throughput sequencing. Although this report focuses on the application of this technique to yeast, the approach can be used with any cells that can be sorted in a flow cytometer according to DNA content.

The effects of transcription and recombination on mutational dynamics of short tandem repeats

Short tandem repeats (STR) are ubiquitous components of the genomic architecture of most living organisms. Recent work has highlighted the widespread functional significance of such repeats, particularly around gene regulation, but the mutational processes underlying the evolution of these highly abundant and highly variable sequences are not fully understood. Traditional models assume that strand misalignment during replication is the predominant mechanism, but empirical data suggest the involvement of other processes including recombination and transcription. Despite this evidence, the relative influences of these processes have not previously been tested experimentally on a genome-wide scale. Using deep sequencing, we identify mutations at >200 microsatellites, across 700 generations in replicated populations of two otherwise identical sexual and asexual Saccharomyces cerevisiae strains. Using generalized linear models, we investigate correlates of STR mutability including the nature of the mutation, STR composition and contextual factors including recombination, transcription and replication origins. Sexual capability was not a significant predictor of microsatellite mutability, but, intriguingly, we identify transcription as a significant positive predictor. We also find that STR density is substantially increased in regions neighboring, but not within, recombination hotspots.

Conservation and Variation in Strategies for DNA Replication of Kinetoplastid Nuclear Genomes

Introduction:

Understanding how the nuclear genome of kinetoplastid parasites is replicated received experimental stimulus from sequencing of the Leishmania major, Trypanosoma brucei and Trypanosoma cruzi genomes around 10 years ago. Gene annotations suggested key players in DNA replication initiation could not be found in these organisms, despite considerable conservation amongst characterised eukaryotes. Initial studies that indicated trypanosomatids might possess an archaeal-like Origin Recognition Complex (ORC), composed of only a single factor termed ORC1/CDC6, have been supplanted by the more recent identification of an ORC in T. brucei. However, the constituent subunits of T. brucei ORC are highly diverged relative to other eukaryotic ORCs and the activity of the complex appears subject to novel, positive regulation. The availability of whole genome sequences has also allowed the deployment of genome-wide strategies to map DNA replication dynamics, to date in T. brucei and Leishmania. ORC1/CDC6 binding and function in T. brucei displays pronounced overlap with the unconventional organisation of gene expression in the genome. Moreover, mapping of sites of replication initiation suggests pronounced differences in replication dynamics in Leishmania relative to T. brucei.

Conclusion:

Here we discuss what implications these emerging data may have for parasite and eukaryotic biology of DNA replication.

eRP arrangement: a strategy for assembled genomic contig rearrangement based on replication profiling in bacteria

Background

The reduced cost of sequencing has made de novo sequencing and the assembly of draft microbial genomes feasible in any ordinary biology lab. However, the process of finishing and completing the genome remains labor-intensive and computationally challenging in some cases, such as in the study of complete genome sequences, genomic rearrangements, long-range syntenic relationships, and structural variations.

Methods

Here, we show a contig reordering strategy based on experimental replication profiling (eRP) to recapitulate the bacterial genome structure within draft genomes. During the exponential growth phase, the majority of bacteria show a global genomic copy number gradient that is enriched near the replication origin and gradually declines toward the terminus. Therefore, if genome sequencing is performed with appropriate timing, the short-read coverage reflects this copy number gradient, providing information about the contig positions relative to the replication origin and terminus.

Results

We therefore investigated the appropriate timing for genomic DNA sampling and developed an algorithm for the reordering of the contigs based on eRP. As a result, this strategy successfully recapitulates the genomic structure of various structural mutants with draft genome sequencing.

Conclusions

Our strategy was successful for contig rearrangement with intracellular DNA replication behavior mechanisms and can be applied to almost all bacteria because the DNA replication system is highly conserved. Therefore, eRP makes it possible to understand genomic structural information and long-range syntenic relationships using a draft genome that is based on short reads.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4162-z) contains supplementary material, which is available to authorized users.

A stitch in time: Replicate early and escape dosage compensation to express more

María Gómez previews work by Müller and Nieduszynski that describes a biological role for conserved replication timing differences in eukaryotic cells.

The biological significance of conserved replication timing patterns in eukaryotic genomes remains a mystery. In this issue, Müller and Nieduszynski (2017. J. Cell Biol.https://doi.org/10.1083/jcb.201701061) find that early replication is a requirement for the highest expression levels of certain genes.

DNA replication timing influences gene expression level

Eukaryotic genomes are replicated in a reproducible temporal order whose physiological significance is poorly understood. Müller and Nieduszynski compare the temporal order of genome replication in phylogenetically diverse yeast species and identify genes for which conserved replication timing contributes to maximal expression.

Eukaryotic genomes are replicated in a reproducible temporal order; however, the physiological significance is poorly understood. We compared replication timing in divergent yeast species and identified genomic features with conserved replication times. Histone genes were among the earliest replicating loci in all species. We specifically delayed the replication of HTA1-HTB1 and discovered that this halved the expression of these histone genes. Finally, we showed that histone and cell cycle genes in general are exempt from Rtt109-dependent dosage compensation, suggesting the existence of pathways excluding specific loci from dosage compensation mechanisms. Thus, we have uncovered one of the first physiological requirements for regulated replication time and demonstrated a direct link between replication timing and gene expression.

Differential paralog divergence modulates genome evolution across yeast species

Evolutionary outcomes depend not only on the selective forces acting upon a species, but also on the genetic background. However, large timescales and uncertain historical selection pressures can make it difficult to discern such important background differences between species. Experimental evolution is one tool to compare evolutionary potential of known genotypes in a controlled environment. Here we utilized a highly reproducible evolutionary adaptation in Saccharomyces cerevisiae to investigate whether experimental evolution of other yeast species would select for similar adaptive mutations. We evolved populations of S. cerevisiae, S. paradoxus, S. mikatae, S. uvarum, and interspecific hybrids between S. uvarum and S. cerevisiae for ~200–500 generations in sulfate-limited continuous culture. Wild-type S. cerevisiae cultures invariably amplify the high affinity sulfate transporter gene, SUL1. However, while amplification of the SUL1 locus was detected in S. paradoxus and S. mikatae populations, S. uvarum cultures instead selected for amplification of the paralog, SUL2. We measured the relative fitness of strains bearing deletions and amplifications of both SUL genes from different species, confirming that, converse to S. cerevisiae, S. uvarum SUL2 contributes more to fitness in sulfate limitation than S. uvarum SUL1. By measuring the fitness and gene expression of chimeric promoter-ORF constructs, we were able to delineate the cause of this differential fitness effect primarily to the promoter of S. uvarum SUL1. Our data show evidence of differential sub-functionalization among the sulfate transporters across Saccharomyces species through recent changes in noncoding sequence. Furthermore, these results show a clear example of how such background differences due to paralog divergence can drive changes in genome evolution.

Both comparative genomics and experimental evolution are powerful tools that can be used to make inferences about evolutionary processes. Together, these approaches provide the opportunity to observe evolutionary adaptation over millions of years where selective history is largely unknown, and over short timescales under controlled selective pressures in the laboratory. We have used comparative experimental evolution to observe the evolutionary fate of an adaptive mutation, and determined to what degree the outcome is conditional on the genetic background. We evolved several populations of different yeast species for over 200 generations in sulfate-limited conditions to determine how the differences in genomic context can alter evolutionary routes when challenged with a nutrient limitation selection pressure. We find that the gene encoding a high affinity sulfur transporter becomes amplified in most species of Saccharomyces, except in S. uvarum, in which the amplification of the paralogous sulfate transporter gene SUL2 is recovered. We attribute this change in amplification preference to mutations in the non-coding region of SUL1, likely due to reduced expression of this gene in S. uvarum. We conclude that the adaptive mutations selected for in each organism depend on the genomic context, even when faced with the same environmental condition.

SIR2 suppresses replication gaps and genome instability by balancing replication between repetitive and unique sequences

Replication gaps that persist into mitosis likely represent important threats to genome stability, but experimental identification of these gaps has proved challenging. We have developed a technique that allows us to explore the dynamics by which genome replication is completed before mitosis. Using this approach, we demonstrate that excessive allocation of replication resources to origins within repetitive regions, induced by SIR2 deletion, leads to persistent replication gaps and genome instability. Conversely, the weakening of replication origins in repetitive regions suppresses these gaps. Given known age- and cancer-associated changes in chromatin accessibility at repetitive sequences, we suggest that replication gaps resulting from misallocation of replication resources underlie age- and disease-associated genome instability.

Model-based analysis of DNA replication profiles: predicting replication fork velocity and initiation rate by profiling free-cycling cells

Eukaryotic cells initiate DNA synthesis by sequential firing of hundreds of origins. This ordered replication is described by replication profiles, which measure the DNA content within a cell population. Here, we show that replication dynamics can be deduced from replication profiles of free-cycling cells. While such profiles lack explicit temporal information, they are sensitive to fork velocity and initiation capacity through the passive replication pattern, namely the replication of origins by forks emanating elsewhere. We apply our model-based approach to a compendium of profiles that include most viable budding yeast mutants implicated in replication. Predicted changes in fork velocity or initiation capacity are verified by profiling synchronously replicating cells. Notably, most mutants implicated in late (or early) origin effects are explained by global modulation of fork velocity or initiation capacity. Our approach provides a rigorous framework for analyzing DNA replication profiles of free-cycling cells.

Mechanisms Governing DDK Regulation of the Initiation of DNA Replication

The budding yeast Dbf4-dependent kinase (DDK) complex—comprised of cell division cycle (Cdc7) kinase and its regulatory subunit dumbbell former 4 (Dbf4)—is required to trigger the initiation of DNA replication through the phosphorylation of multiple minichromosome maintenance complex subunits 2-7 (Mcm2-7). DDK is also a target of the radiation sensitive 53 (Rad53) checkpoint kinase in response to replication stress. Numerous investigations have determined mechanistic details, including the regions of Mcm2, Mcm4, and Mcm6 phosphorylated by DDK, and a number of DDK docking sites. Similarly, the way in which the Rad53 forkhead-associated 1 (FHA1) domain binds to DDK—involving both canonical and non-canonical interactions—has been elucidated. Recent work has revealed mutual promotion of DDK and synthetic lethal with dpb11-1 3 (Sld3) roles. While DDK phosphorylation of Mcm2-7 subunits facilitates their interaction with Sld3 at origins, Sld3 in turn stimulates DDK phosphorylation of Mcm2. Details of a mutually antagonistic relationship between DDK and Rap1-interacting factor 1 (Rif1) have also recently come to light. While Rif1 is able to reverse DDK-mediated Mcm2-7 complex phosphorylation by targeting the protein phosphatase glycogen 7 (Glc7) to origins, there is evidence to suggest that DDK can counteract this activity by binding to and phosphorylating Rif1.

Chromosomal Replication Complexity: A Novel DNA Metrics and Genome Instability Factor

As the ratio of the copy number of the most replicated to the unreplicated regions in the same chromosome, the definition of chromosomal replication complexity (CRC) appears to leave little room for variation, being either two during S-phase or one otherwise. However, bacteria dividing faster than they replicate their chromosome spike CRC to four and even eight. A recent experimental inquiry about the limits of CRC in Escherichia coli revealed two major reasons to avoid elevating it further: (i) increased chromosomal fragmentation and (ii) complications with subsequent double-strand break repair. Remarkably, examples of stable elevated CRC in eukaryotic chromosomes are well known under various terms like "differential replication," "underreplication," "DNA puffs," "onion-skin replication," or "re-replication" and highlight the phenomenon of static replication fork (sRF). To accurately describe the resulting "amplification by overinitiation," I propose a new term: "replification" (subchromosomal overreplication). In both prokaryotes and eukaryotes, replification, via sRF processing, causes double-strand DNA breaks and, with their repair elevating chromosomal rearrangements, represents a novel genome instability factor. I suggest how static replication bubbles could be stabilized and speculate that some tandem duplications represent such persistent static bubbles. Moreover, I propose how static replication bubbles could be transformed into tandem duplications, double minutes, or inverted triplications. Possible experimental tests of these models are discussed.

High Throughput Analyses of Budding Yeast ARSs Reveal New DNA Elements Capable of Conferring Centromere-Independent Plasmid Propagation

The ability of plasmids to propagate in Saccharomyces cerevisiae has been instrumental in defining eukaryotic chromosomal control elements. Stable propagation demands both plasmid replication, which requires a chromosomal replication origin (i.e., an ARS), and plasmid distribution to dividing cells, which requires either a chromosomal centromere for segregation or a plasmid-partitioning element. While our knowledge of yeast ARSs and centromeres is relatively advanced, we know less about chromosomal regions that can function as plasmid partitioning elements. The Rap1 protein-binding site (RAP1) present in transcriptional silencers and telomeres of budding yeast is a known plasmid-partitioning element that functions to anchor a plasmid to the inner nuclear membrane (INM), which in turn facilitates plasmid distribution to daughter cells. This Rap1-dependent INM-anchoring also has an important chromosomal role in higher-order chromosomal structures that enhance transcriptional silencing and telomere stability. Thus, plasmid partitioning can reflect fundamental features of chromosome structure and biology, yet a systematic screen for plasmid partitioning elements has not been reported. Here, we couple deep sequencing with competitive growth experiments of a plasmid library containing thousands of short ARS fragments to identify new plasmid partitioning elements. Competitive growth experiments were performed with libraries that differed only in terms of the presence or absence of a centromere. Comparisons of the behavior of ARS fragments in the two experiments allowed us to identify sequences that were likely to drive plasmid partitioning. In addition to the silencer RAP1 site, we identified 74 new putative plasmid-partitioning motifs predicted to act as binding sites for DNA binding proteins enriched for roles in negative regulation of gene expression and G2/M-phase associated biology. These data expand our knowledge of chromosomal elements that may function in plasmid partitioning and suggest underlying biological roles shared by such elements.

A Versatile Procedure to Generate Genome-Wide Spatiotemporal Program of Replication in Yeast Species

Here, we describe a complete protocol, comprising both the experimental and the analytical procedures, that allows to generate genome-wide spatiotemporal program of replication and to find the location of chromosomally active replication origins in yeast. The first step consists on synchronizing a cell population by physical discrimination of G1 cells according to their sedimentation coefficient. G1 cells are then synchronously released into S-phase and time-point samples are regularly taken until they reach the G2 phase. Progression through the cell cycle is monitored by measuring DNA content variation by flow cytometry. DNA samples, covering the entire S-phase, are then extracted and analyzed using deep sequencing. The gradual change of DNA copy number is measured to determine the mean replication time along the genome. A simple method of peak calling allows to infer from the replication profile the location of replication origins along the chromosomes. Our protocol is versatile enough to be applied to virtually any yeast species of interest and generate its replication profile.

Large-Scale Chromatin Structure-Function Relationships during the Cell Cycle and Development: Insights from Replication Timing

Chromosome architecture has received a lot of attention since the recent development of genome-scale methods to measure chromatin interactions (Hi-C), enabling the first sequence-based models of chromosome tertiary structure. A view has emerged of chromosomes as a string of structural units (topologically associating domains; TADs) whose boundaries persist through the cell cycle and development. TADs with similar chromatin states tend to aggregate, forming spatially segregated chromatin compartments. However, high-resolution Hi-C has revealed substructure within TADs (subTADs) that poses a challenge for models that attribute significance to structural units at any given scale. More than 20 years ago, the DNA replication field independently identified stable structural (and functional) units of chromosomes (replication foci) as well as spatially segregated chromatin compartments (early and late foci), but lacked the means to link these units to genomic map units. Genome-wide studies of replication timing (RT) have now merged these two disciplines by identifying individual units of replication regulation (replication domains; RDs) that correspond to TADs and are arranged in 3D to form spatiotemporally segregated subnuclear compartments. Furthermore, classifying RDs/TADs by their constitutive versus developmentally regulated RT has revealed distinct classes of chromatin organization, providing unexpected insight into the relationship between large-scale chromosome structure and function.

Genome-wide mapping reveals single-origin chromosome replication in Leishmania, a eukaryotic microbe

BACKGROUND

DNA replication initiates on defined genome sites, termed origins. Origin usage appears to follow common rules in the eukaryotic organisms examined to date: all chromosomes are replicated from multiple origins, which display variations in firing efficiency and are selected from a larger pool of potential origins. To ask if these features of DNA replication are true of all eukaryotes, we describe genome-wide origin mapping in the parasite Leishmania.

RESULTS

Origin mapping in Leishmania suggests a striking divergence in origin usage relative to characterized eukaryotes, since each chromosome appears to be replicated from a single origin. By comparing two species of Leishmania, we find evidence that such origin singularity is maintained in the face of chromosome fusion or fission events during evolution. Mapping Leishmania origins suggests that all origins fire with equal efficiency, and that the genomic sites occupied by origins differ from related non-origins sites. Finally, we provide evidence that origin location in Leishmania displays striking conservation with Trypanosoma brucei, despite the latter parasite replicating its chromosomes from multiple, variable strength origins.

CONCLUSIONS

The demonstration of chromosome replication for a single origin in Leishmania, a microbial eukaryote, has implications for the evolution of origin multiplicity and associated controls, and may explain the pervasive aneuploidy that characterizes Leishmania chromosome architecture.

Genome-wide replication landscape of Candida glabrata

Background

The opportunistic pathogen Candida glabrata is a member of the Saccharomycetaceae yeasts. Like its close relative Saccharomyces cerevisiae, it underwent a whole-genome duplication followed by an extensive loss of genes. Its genome contains a large number of very long tandem repeats, called megasatellites. In order to determine the whole replication program of the C. glabrata genome and its general chromosomal organization, we used deep-sequencing and chromosome conformation capture experiments.

Results

We identified 253 replication fork origins, genome wide. Centromeres, HML and HMR loci, and most histone genes are replicated early, whereas natural chromosomal breakpoints are located in late-replicating regions. In addition, 275 autonomously replicating sequences (ARS) were identified during ARS-capture experiments, and their relative fitness was determined during growth competition. Analysis of ARSs allowed us to identify a 17-bp consensus, similar to the S. cerevisiae ARS consensus sequence but slightly more constrained. Megasatellites are not in close proximity to replication origins or termini. Using chromosome conformation capture, we also show that early origins tend to cluster whereas non-subtelomeric megasatellites do not cluster in the yeast nucleus.

Conclusions

Despite a shorter cell cycle, the C. glabrata replication program shares unexpected striking similarities to S. cerevisiae, in spite of their large evolutionary distance and the presence of highly repetitive large tandem repeats in C. glabrata. No correlation could be found between the replication program and megasatellites, suggesting that their formation and propagation might not be directly caused by replication fork initiation or termination.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0177-6) contains supplementary material, which is available to authorized users.