A unifying model that explains the origins of human inverted copy number variants.

With the release of the telomere-to-telomere human genome sequence and the availability of both long-read sequencing and optical genome mapping techniques, the identification of copy number variants (CNVs) and other structural variants is providing new insights into human genetic disease. Different mechanisms have been proposed to account for the novel junctions in these complex architectures, including aberrant forms of DNA replication, non-allelic homologous recombination, and various pathways that repair DNA breaks. Here, we have focused on a set of structural variants that include an inverted segment and propose that they share a common initiating event: an inverted triplication with long, unstable palindromic junctions. The secondary rearrangement of these palindromes gives rise to the various forms of inverted structural variants. We postulate that this same mechanism (ODIRA: origin-dependent inverted-repeat amplification) that creates the inverted CNVs in inherited syndromes also generates the palindromes found in cancers.

Template switching between the leading and lagging strands at replication forks generates inverted copy number variants through hairpin-capped extrachromosomal DNA.

Inherited and germ-line de novo copy number variants (CNVs) are increasingly found to be correlated with human developmental and cancerous phenotypes. Several models for template switching during replication have been proposed to explain the generation of these gross chromosomal rearrangements. We proposed a model of template switching (ODIRA-origin dependent inverted repeat amplification) in which simultaneous ligation of the leading and lagging strands at diverging replication forks could generate segmental inverted triplications through an extrachromosomal inverted circular intermediate. Here, we created a genetic assay using split-ura3 cassettes to trap the proposed inverted intermediate. However, instead of recovering circular inverted intermediates, we found inverted linear chromosomal fragments ending in native telomeres-suggesting that a template switch had occurred at the centromere-proximal fork of a replication bubble. As telomeric inverted hairpin fragments can also be created through double strand breaks we tested whether replication errors or repair of double stranded DNA breaks were the most likely initiating event. The results from CRISPR/Cas9 cleavage experiments and growth in the replication inhibitor hydroxyurea indicate that it is a replication error, not a double stranded break that creates the inverted junctions. Since inverted amplicons of the SUL1 gene occur during long-term growth in sulfate-limited chemostats, we sequenced evolved populations to look for evidence of linear intermediates formed by an error in replication. All of the data are compatible with a two-step version of the ODIRA model in which sequential template switching at short inverted repeats between the leading and lagging strands at a replication fork, followed by integration via homologous recombination, generates inverted interstitial triplications.

The effects of manipulating levels of replication initiation factors on origin firing efficiency in yeast.

Chromosome replication in Saccharomyces cerevisiae is initiated from ~300 origins that are regulated by DNA sequence and by the limited abundance of six trans-acting initiation proteins (Sld2, Sld3, Dpb11, Dbf4, Sld7 and Cdc45). We set out to determine how the levels of individual factors contribute to time of origin activation and/or origin efficiency using induced depletion of single factors and overexpression of sets of multiple factors. Depletion of Sld2 or Sld3 slows growth and S phase progression, decreases origin efficiency across the genome and impairs viability as a result of incomplete replication of the rDNA. We find that the most efficient early origins are relatively unaffected by depletion of either Sld2 or Sld3. However, Sld3 levels, and to a lesser extent Sld2 levels, are critical for firing of the less efficient early origins. Overexpression of Sld3 simultaneously with Sld2, Dpb11 and Dbf4 preserves the relative efficiency of origins. Only when Cdc45 and Sld7 are also overexpressed is origin efficiency equalized between early- and late-firing origins. Our data support a model in which Sld3 together with Cdc45 (and/or Sld7) is responsible for the differential efficiencies of origins across the yeast genome.

Rif1 controls DNA replication by directing Protein Phosphatase 1 to reverse Cdc7-mediated phosphorylation of the MCM complex.

Initiation of eukaryotic DNA replication requires phosphorylation of the MCM complex by Dbf4-dependent kinase (DDK), composed of Cdc7 kinase and its activator, Dbf4. We report here that budding yeast Rif1 (Rap1-interacting factor 1) controls DNA replication genome-wide and describe how Rif1 opposes DDK function by directing Protein Phosphatase 1 (PP1)-mediated dephosphorylation of the MCM complex. Deleting RIF1 partially compensates for the limited DDK activity in a cdc7-1 mutant strain by allowing increased, premature phosphorylation of Mcm4. PP1 interaction motifs within the Rif1 N-terminal domain are critical for its repressive effect on replication. We confirm that Rif1 interacts with PP1 and that PP1 prevents premature Mcm4 phosphorylation. Remarkably, our results suggest that replication repression by Rif1 is itself also DDK-regulated through phosphorylation near the PP1-interacting motifs. Based on our findings, we propose that Rif1 is a novel PP1 substrate targeting subunit that counteracts DDK-mediated phosphorylation during replication. Fission yeast and mammalian Rif1 proteins have also been implicated in regulating DNA replication. Since PP1 interaction sites are evolutionarily conserved within the Rif1 sequence, it is likely that replication control by Rif1 through PP1 is a conserved mechanism.

Defective replication initiation results in locus specific chromosome breakage and a ribosomal RNA deficiency in yeast.

A form of dwarfism known as Meier-Gorlin syndrome (MGS) is caused by recessive mutations in one of six different genes (ORC1, ORC4, ORC6, CDC6, CDT1, and MCM5). These genes encode components of the pre-replication complex, which assembles at origins of replication prior to S phase. Also, variants in two additional replication initiation genes have joined the list of causative mutations for MGS (Geminin and CDC45). The identity of the causative MGS genetic variants strongly suggests that some aspect of replication is amiss in MGS patients; however, little evidence has been obtained regarding what aspect of chromosome replication is faulty. Since the site of one of the missense mutations in the human ORC4 alleles is conserved between humans and yeast, we sought to determine in what way this single amino acid change affects the process of chromosome replication, by introducing the comparable mutation into yeast (orc4Y232C). We find that yeast cells with the orc4Y232C allele have a prolonged S-phase, due to compromised replication initiation at the ribosomal DNA (rDNA) locus located on chromosome XII. The inability to initiate replication at the rDNA locus results in chromosome breakage and a severely reduced rDNA copy number in the survivors, presumably helping to ensure complete replication of chromosome XII. Although reducing rDNA copy number may help ensure complete chromosome replication, orc4Y232C cells struggle to meet the high demand for ribosomal RNA synthesis. This finding provides additional evidence linking two essential cellular pathways-DNA replication and ribosome biogenesis.

The B. subtilis Accessory Helicase PcrA Facilitates DNA Replication through Transcription Units.

In bacteria the concurrence of DNA replication and transcription leads to potentially deleterious encounters between the two machineries, which can occur in either the head-on (lagging strand genes) or co-directional (leading strand genes) orientations. These conflicts lead to replication fork stalling and can destabilize the genome. Both eukaryotic and prokaryotic cells possess resolution factors that reduce the severity of these encounters. Though Escherichia coli accessory helicases have been implicated in the mitigation of head-on conflicts, direct evidence of these proteins mitigating co-directional conflicts is lacking. Furthermore, the endogenous chromosomal regions where these helicases act, and the mechanism of recruitment, have not been identified. We show that the essential Bacillus subtilis accessory helicase PcrA aids replication progression through protein coding genes of both head-on and co-directional orientations, as well as rRNA and tRNA genes. ChIP-Seq experiments show that co-directional conflicts at highly transcribed rRNA, tRNA, and head-on protein coding genes are major targets of PcrA activity on the chromosome. Partial depletion of PcrA renders cells extremely sensitive to head-on conflicts, linking the essential function of PcrA to conflict resolution. Furthermore, ablating PcrA's ATPase/helicase activity simultaneously increases its association with conflict regions, while incapacitating its ability to mitigate conflicts, and leads to cell death. In contrast, disruption of PcrA's C-terminal RNA polymerase interaction domain does not impact its ability to mitigate conflicts between replication and transcription, its association with conflict regions, or cell survival. Altogether, this work establishes PcrA as an essential factor involved in mitigating transcription-replication conflicts and identifies chromosomal regions where it routinely acts. As both conflicts and accessory helicases are found in all domains of life, these results are broadly relevant.

Origin-Dependent Inverted-Repeat Amplification: Tests of a Model for Inverted DNA Amplification.

DNA replication errors are a major driver of evolution--from single nucleotide polymorphisms to large-scale copy number variations (CNVs). Here we test a specific replication-based model to explain the generation of interstitial, inverted triplications. While no genetic information is lost, the novel inversion junctions and increased copy number of the included sequences create the potential for adaptive phenotypes. The model--Origin-Dependent Inverted-Repeat Amplification (ODIRA)-proposes that a replication error at pre-existing short, interrupted, inverted repeats in genomic sequences generates an extrachromosomal, inverted dimeric, autonomously replicating intermediate; subsequent genomic integration of the dimer yields this class of CNV without loss of distal chromosomal sequences. We used a combination of in vitro and in vivo approaches to test the feasibility of the proposed replication error and its downstream consequences on chromosome structure in the yeast Saccharomyces cerevisiae. We show that the proposed replication error-the ligation of leading and lagging nascent strands to create "closed" forks-can occur in vitro at short, interrupted inverted repeats. The removal of molecules with two closed forks results in a hairpin-capped linear duplex that we show replicates in vivo to create an inverted, dimeric plasmid that subsequently integrates into the genome by homologous recombination, creating an inverted triplication. While other models have been proposed to explain inverted triplications and their derivatives, our model can also explain the generation of human, de novo, inverted amplicons that have a 2:1 mixture of sequences from both homologues of a single parent--a feature readily explained by a plasmid intermediate that arises from one homologue and integrates into the other homologue prior to meiosis. Our tests of key features of ODIRA lend support to this mechanism and suggest further avenues of enquiry to unravel the origins of interstitial, inverted CNVs pivotal in human health and evolution.

GC-rich DNA elements enable replication origin activity in the methylotrophic yeast Pichia pastoris.

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins--a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

A natural polymorphism in rDNA replication origins links origin activation with calorie restriction and lifespan.

Aging and longevity are complex traits influenced by genetic and environmental factors. To identify quantitative trait loci (QTLs) that control replicative lifespan, we employed an outbred Saccharomyces cerevisiae model, generated by crossing a vineyard and a laboratory strain. The predominant QTL mapped to the rDNA, with the vineyard rDNA conferring a lifespan increase of 41%. The lifespan extension was independent of Sir2 and Fob1, but depended on a polymorphism in the rDNA origin of replication from the vineyard strain that reduced origin activation relative to the laboratory origin. Strains carrying vineyard rDNA origins have increased capacity for replication initiation at weak plasmid and genomic origins, suggesting that inability to complete genome replication presents a major impediment to replicative lifespan. Calorie restriction, a conserved mediator of lifespan extension that is also independent of Sir2 and Fob1, reduces rDNA origin firing in both laboratory and vineyard rDNA. Our results are consistent with the possibility that calorie restriction, similarly to the vineyard rDNA polymorphism, modulates replicative lifespan through control of rDNA origin activation, which in turn affects genome replication dynamics.

Maintaining replication origins in the face of genomic change.

Origins of replication present a paradox to evolutionary biologists. As a collection, they are absolutely essential genomic features, but individually are highly redundant and nonessential. It is therefore difficult to predict to what extent and in what regard origins are conserved over evolutionary time. Here, through a comparative genomic analysis of replication origins and chromosomal replication patterns in the budding yeasts Saccharomyces cerevisiae and Lachancea waltii, we assess to what extent replication origins survived genomic change produced from 150 million years of evolution. We find that L. waltii origins exhibit a core consensus sequence and nucleosome occupancy pattern highly similar to those of S. cerevisiae origins. We further observe that the overall progression of chromosomal replication is similar between L. waltii and S. cerevisiae. Nevertheless, few origins show evidence of being conserved in location between the two species. Among the conserved origins are those surrounding centromeres and adjacent to histone genes, suggesting that proximity to an origin may be important for their regulation. We conclude that, over evolutionary time, origins maintain sequence, structure, and regulation, but are continually being created and destroyed, with the result that their locations are generally not conserved.

Functional centromeres determine the activation time of pericentric origins of DNA replication in Saccharomyces cerevisiae.

The centromeric regions of all Saccharomyces cerevisiae chromosomes are found in early replicating domains, a property conserved among centromeres in fungi and some higher eukaryotes. Surprisingly, little is known about the biological significance or the mechanism of early centromere replication; however, the extensive conservation suggests that it is important for chromosome maintenance. Do centromeres ensure their early replication by promoting early activation of nearby origins, or have they migrated over evolutionary time to reside in early replicating regions? In Candida albicans, a neocentromere contains an early firing origin, supporting the first hypothesis but not addressing whether the new origin is intrinsically early firing or whether the centromere influences replication time. Because the activation time of individual origins is not an intrinsic property of S. cerevisiae origins, but is influenced by surrounding sequences, we sought to test the hypothesis that centromeres influence replication time by moving a centromere to a late replication domain. We used a modified Meselson-Stahl density transfer assay to measure the kinetics of replication for regions of chromosome XIV in which either the functional centromere or a point-mutated version had been moved near origins that reside in a late replication region. We show that a functional centromere acts in cis over a distance as great as 19 kb to advance the initiation time of origins. Our results constitute a direct link between establishment of the kinetochore and the replication initiation machinery, and suggest that the proposed higher-order structure of the pericentric chromatin influences replication initiation.

Genomic mapping of single-stranded DNA in hydroxyurea-challenged yeasts identifies origins of replication.

During DNA replication one or both strands transiently become single stranded: first at the sites where initiation of DNA synthesis occurs (known as origins of replication) and subsequently on the lagging strands of replication forks as discontinuous Okazaki fragments are generated. We report a genome-wide analysis of single-stranded DNA (ssDNA) formation in the presence of hydroxyurea during DNA replication in wild-type and checkpoint-deficient rad53 Saccharomyces cerevisiae cells. In wild-type cells, ssDNA was first observed at a subset of replication origins and later 'migrated' bi-directionally, suggesting that ssDNA formation is associated with continuously moving replication forks. In rad53 cells, ssDNA was observed at virtually every known origin, but remained there over time, suggesting that replication forks stall. Telomeric regions seemed to be particularly sensitive to the loss of Rad53 checkpoint function. Replication origins in Schizosaccharomyces pombe were also mapped using our method.

Hotspot of de novo telomere addition stabilizes linear amplicons in yeast grown in sulfate-limiting conditions.

Evolution is driven by the accumulation of competing mutations that influence survival. A broad form of genetic variation is the amplification or deletion of DNA (≥50 bp) referred to as copy number variation (CNV). In humans, CNV may be inconsequential, contribute to minor phenotypic differences, or cause conditions such as birth defects, neurodevelopmental disorders, and cancers. To identify mechanisms that drive CNV, we monitored the experimental evolution of Saccharomyces cerevisiae populations grown under sulfate-limiting conditions. Cells with increased copy number of the gene SUL1, which encodes a primary sulfate transporter, exhibit a fitness advantage. Previously, we reported interstitial inverted triplications of SUL1 as the dominant rearrangement in a haploid population. Here, in a diploid population, we find instead that small linear fragments containing SUL1 form and are sustained over several generations. Many of the linear fragments are stabilized by de novo telomere addition within a telomere-like sequence near SUL1 (within the SNF5 gene). Using an assay that monitors telomerase action following an induced chromosome break, we show that this region acts as a hotspot of de novo telomere addition and that required sequences map to a region of <250 base pairs. Consistent with previous work showing that association of the telomere-binding protein Cdc13 with internal sequences stimulates telomerase recruitment, mutation of a four-nucleotide motif predicted to associate with Cdc13 abolishes de novo telomere addition. Our study suggests that internal telomere-like sequences that stimulate de novo telomere addition can contribute to adaptation by promoting genomic plasticity.

Ribosomal DNA replication time coordinates completion of genome replication and anaphase in yeast.

Timely completion of genome replication is a prerequisite for mitosis, genome integrity, and cell survival. A challenge to this timely completion comes from the need to replicate the hundreds of untranscribed copies of rDNA that organisms maintain in addition to the copies required for ribosome biogenesis. Replication of these rDNA arrays is relegated to late S phase despite their large size, repetitive nature, and essentiality. Here, we show that, in Saccharomyces cerevisiae, reducing the number of rDNA repeats leads to early rDNA replication, which results in delaying replication elsewhere in the genome. Moreover, cells with early-replicating rDNA arrays and delayed genome-wide replication aberrantly release the mitotic phosphatase Cdc14 from the nucleolus and enter anaphase prematurely. We propose that rDNA copy number determines the replication time of the rDNA locus and that the release of Cdc14 upon completion of rDNA replication is a signal for cell cycle progression.

Genetic, genomic, and molecular tools for studying the protoploid yeast, L. waltii.

Sequencing of the yeast Kluyveromyces waltii (recently renamed Lachancea waltii) provided evidence of a whole genome duplication event in the lineage leading to the well-studied Saccharomyces cerevisiae. While comparative genomic analyses of these yeasts have proven to be extremely instructive in modeling the loss or maintenance of gene duplicates, experimental tests of the ramifications following such genome alterations remain difficult. To transform L. waltii from an organism of the computational comparative genomic literature into an organism of the functional comparative genomic literature, we have developed genetic, molecular and genomic tools for working with L. waltii. In particular, we have characterized basic properties of L. waltii (growth, ploidy, molecular karyotype, mating type and the sexual cycle), developed transformation, cell cycle arrest and synchronization protocols, and have created centromeric and non-centromeric vectors as well as a genome browser for L. waltii. We hope that these tools will be used by the community to follow up on the ideas generated by sequence data and lead to a greater understanding of eukaryotic biology and genome evolution.

Molecular analysis of the replication program in unicellular model organisms.

Eukaryotes have long been reported to show temporal programs of replication, different portions of the genome being replicated at different times in S phase, with the added possibility of developmentally regulated changes in this pattern depending on species and cell type. Unicellular model organisms, primarily the budding yeast Saccharomyces cerevisiae, have been central to our current understanding of the mechanisms underlying the regulation of replication origins and the temporal program of replication in particular. But what exactly is a temporal program of replication, and how might it arise? In this article, we explore this question, drawing again on the wealth of experimental information in unicellular model organisms.

Replication in hydroxyurea: it's a matter of time.

Hydroxyurea (HU) is a DNA replication inhibitor that negatively affects both the elongation and initiation phases of replication and triggers the "intra-S phase checkpoint." Previous work with budding yeast has shown that, during a short exposure to HU, MEC1/RAD53 prevent initiation at some late S phase origins. In this study, we have performed microarray experiments to follow the fate of all origins over an extended exposure to HU. We show that the genome-wide progression of DNA synthesis, including origin activation, follows the same pattern in the presence of HU as in its absence, although the time frames are very different. We find no evidence for a specific effect that excludes initiation from late origins. Rather, HU causes S phase to proceed in slow motion; all temporal classes of origins are affected, but the order in which they become active is maintained. We propose a revised model for the checkpoint response to HU that accounts for the continued but slowed pace of the temporal program of origin activation.

Structural changes in Mcm5 protein bypass Cdc7-Dbf4 function and reduce replication origin efficiency in Saccharomyces cerevisiae.

Eukaryotic chromosomal replication is a complicated process with many origins firing at different efficiencies and times during S phase. Prereplication complexes are assembled on all origins in G(1) phase, and yet only a subset of complexes is activated during S phase by DDK (for Dbf4-dependent kinase) (Cdc7-Dbf4). The yeast mcm5-bob1 (P83L) mutation bypasses DDK but results in reduced intrinsic firing efficiency at 11 endogenous origins and at origins located on minichromosomes. Origin efficiency may result from Mcm5 protein assuming an altered conformation, as predicted from the atomic structure of an archaeal MCM (for minichromosome maintenance) homologue. Similarly, an intragenic mutation in a residue predicted to interact with P83L suppresses the mcm5-bob1 bypass phenotype. We propose DDK phosphorylation of the MCM complex normally results in a single, highly active conformation of Mcm5, whereas the mcm5-bob1 mutation produces a number of conformations, only one of which is permissive for origin activation. Random adoption of these alternate states by the mcm5-bob1 protein can explain both how origin firing occurs independently of DDK and why origin efficiency is reduced. Because similar mutations in mcm2 and mcm4 cannot bypass DDK, Mcm5 protein may be a unique Mcm protein that is the final target of DDK regulation.

Challenges and Approaches to Genotyping Repetitive DNA.

Individuals within a species can exhibit vast variation in copy number of repetitive DNA elements. This variation may contribute to complex traits such as lifespan and disease, yet it is only infrequently considered in genotype-phenotype associations. Although the possible importance of copy number variation is widely recognized, accurate copy number quantification remains challenging. Here, we assess the technical reproducibility of several major methods for copy number estimation as they apply to the large repetitive ribosomal DNA array (rDNA). rDNA encodes the ribosomal RNAs and exists as a tandem gene array in all eukaryotes. Repeat units of rDNA are kilobases in size, often with several hundred units comprising the array, making rDNA particularly intractable to common quantification techniques. We evaluate pulsed-field gel electrophoresis, droplet digital PCR, and Nextera-based whole genome sequencing as approaches to copy number estimation, comparing techniques across model organisms and spanning wide ranges of copy numbers. Nextera-based whole genome sequencing, though commonly used in recent literature, produced high error. We explore possible causes for this error and provide recommendations for best practices in rDNA copy number estimation. We present a resource of high-confidence rDNA copy number estimates for a set of S. cerevisiae and C. elegans strains for future use. We furthermore explore the possibility for FISH-based copy number estimation, an alternative that could potentially characterize copy number on a cellular level.

The temporal program of chromosome replication: genomewide replication in clb5{Delta} Saccharomyces cerevisiae.

Temporal regulation of origin activation is widely thought to explain the pattern of early- and late-replicating domains in the Saccharomyces cerevisiae genome. Recently, single-molecule analysis of replication suggested that stochastic processes acting on origins with different probabilities of activation could generate the observed kinetics of replication without requiring an underlying temporal order. To distinguish between these possibilities, we examined a clb5Delta strain, where origin firing is largely limited to the first half of S phase, to ask whether all origins nonspecifically show decreased firing (as expected for disordered firing) or if only some origins ("late" origins) are affected. Approximately half the origins in the mutant genome show delayed replication while the remainder replicate largely on time. The delayed regions can encompass hundreds of kilobases and generally correspond to regions that replicate late in wild-type cells. Kinetic analysis of replication in wild-type cells reveals broad windows of origin firing for both early and late origins. Our results are consistent with a temporal model in which origins can show some heterogeneity in both time and probability of origin firing, but clustering of temporally like origins nevertheless yields a genome that is organized into blocks showing different replication times.