Segmental duplications arise from Pol32-dependent repair of broken forks through two alternative replication-based mechanisms

Chromosome catastrophes involve replication mechanisms generating complex genomic rearrangements

Complex genomic rearrangements (CGRs) consisting of two or more breakpoint junctions have been observed in genomic disorders. Recently, a chromosome catastrophe phenomenon termed chromothripsis, in which numerous genomic rearrangements are apparently acquired in one single catastrophic event, was described in multiple cancers. Here, we show that constitutionally acquired CGRs share similarities with cancer chromothripsis. In the 17 CGR cases investigated, we observed localization and multiple copy number changes including deletions, duplications, and/or triplications, as well as extensive translocations and inversions. Genomic rearrangements involved varied in size and complexities; in one case, array comparative genomic hybridization revealed 18 copy number changes. Breakpoint sequencing identified characteristic features, including small templated insertions at breakpoints and microhomology at breakpoint junctions, which have been attributed to replicative processes. The resemblance between CGR and chromothripsis suggests similar mechanistic underpinnings. Such chromosome catastrophic events appear to reflect basic DNA metabolism operative throughout an organism's life cycle.

Intermittency as a universal characteristic of the complete chromosome DNA sequences of eukaryotes: from protozoa to human genomes

Large-scale dynamical properties of complete chromosome DNA sequences of eukaryotes are considered. Using the proposed deterministic models with intermittency and symbolic dynamics we describe a wide spectrum of large-scale patterns inherent in these sequences, such as segmental duplications, tandem repeats, and other complex sequence structures. It is shown that the recently discovered gene number balance on the strands is not of a random nature, and certain subsystems of a complete chromosome DNA sequence exhibit the properties of deterministic chaos.

Global chromosomal structural instability in a subpopulation of starving Escherichia coli cells

Copy-number variations (CNVs) constitute very common differences between individual humans and possibly all genomes and may therefore be important fuel for evolution, yet how they form remains elusive. In starving Escherichia coli, gene amplification is induced by stress, controlled by the general stress response. Amplification has been detected only encompassing genes that confer a growth advantage when amplified. We studied the structure of stress-induced gene amplification in starving cells in the Lac assay in Escherichia coli by array comparative genomic hybridization (aCGH), with polymerase chain reaction (pcr) and DNA sequencing to establish the structures generated. About 10% of 300 amplified isolates carried other chromosomal structural change in addition to amplification. Most of these were inversions and duplications associated with the amplification event. This complexity supports a mechanism similar to that seen in human non-recurrent copy number variants. We interpret these complex events in terms of repeated template switching during DNA replication. Importantly, we found a significant occurrence (6 out of 300) of chromosomal structural changes that were apparently not involved in the amplification event. These secondary changes were absent from 240 samples derived from starved cells not carrying amplification, suggesting that amplification happens in a differentiated subpopulation of stressed cells licensed for global chromosomal structural change and genomic instability. These data imply that chromosomal structural changes occur in bursts or showers of instability that may have the potential to drive rapid evolution.

Much of the difference between individual humans is in the number of copies of genes and lengths of genome. The mechanisms by which copy number variation arises are not well understood. We sought information on copy number change mechanisms by extensive use of array comparative genomic hybridization of whole genomes in bacteria selected for amplification of part of the genome. We report that about 10% of amplified isolates carried other chromosomal structural changes associated with the amplification, a result comparable to that seen in human copy number variants. Importantly, we found a significant occurrence of structural changes that were not involved in the amplification event. These were not seen in a control sample of stressed cells not carrying amplification. This establishes that chromosomal structural change happens in a subpopulation of cells apparently licensed to undergo these changes. Because the changes occur under the stress of starvation and require two of the cells' stress-response systems, we propose that licensing for cell-wide structural change in this subpopulation is a component of response to stress. This idea has implications for the mechanisms of evolution and cancer progression, suggesting that changes occur in a shower of events rather than as isolated random events.

The evolution of gene expression regulatory networks in yeasts

Gene regulation is a major source of phenotypic diversity between and within species. This aspect of evolution has long been addressed from the sole point of view of the genome sequence. The incredible development of transcriptomics approaches now allows one to actually study the topology and the properties of regulatory networks on an evolutionary perspective. This new discipline is called comparative functional genomics or comparative transcriptomics. This article reviews some of the main advances made in this field, using yeast species, and especially the species sequenced in the frame of the Genolevures program, as a model.

Repeat expansion in the budding yeast ribosomal DNA can occur independently of the canonical homologous recombination machinery

Major eukaryotic genomic elements, including the ribosomal DNA (rDNA), are composed of repeated sequences with well-defined copy numbers that must be maintained by regulated recombination. Although mechanisms that instigate rDNA recombination have been identified, none are directional and they therefore cannot explain precise repeat number control. Here, we show that yeast lacking histone chaperone Asf1 undergo reproducible rDNA repeat expansions. These expansions do not require the replication fork blocking protein Fob1 and are therefore independent of known rDNA expansion mechanisms. We propose the existence of a regulated rDNA repeat gain pathway that becomes constitutively active in asf1Δ mutants. Cells lacking ASF1 accumulate rDNA repeats with high fidelity in a processive manner across multiple cell divisions. The mechanism of repeat gain is dependent on highly repetitive sequence but, surprisingly, is independent of the homologous recombination proteins Rad52, Rad51 and Rad59. The expansion mechanism is compromised by mutations that decrease the processivity of DNA replication, which leads to progressive loss of rDNA repeats. Our data suggest that a novel mode of break-induced replication occurs in repetitive DNA that is dependent on high homology but does not require the canonical homologous recombination machinery.

Divergent roles for the two PolI-like organelle DNA polymerases of Arabidopsis

DNA polymerases play a central role in the process of DNA replication. Yet, the proteins in charge of the replication of plant organelle DNA have not been unambiguously identified. There are however many indications that a family of proteins homologous to bacterial DNA polymerase I (PolI) is implicated in organelle DNA replication. Here, we have isolated mutant lines of the PolIA and PolIB genes of Arabidopsis (Arabidopsis thaliana) to test this hypothesis. We find that mutation of both genes is lethal, thus confirming an essential and redundant role for these two proteins. However, the mutation of a single gene is sufficient to cause a reduction in the levels of DNA in both mitochondria and plastids. We also demonstrate that polIb, but not polIa mutant lines, are hypersensitive to ciprofloxacin, a small molecule that specifically induces DNA double-strand breaks in plant organelles, suggesting a function for PolIB in DNA repair. In agreement with this result, a cross between polIb and a plastid Whirly mutant line yielded plants with high levels of DNA rearrangements and severe growth defects, indicating impairments in plastid DNA repair pathways. Taken together, this work provides further evidences for the involvement of the plant PolI-like genes in organelle DNA replication and suggests an additional role for PolIB in DNA repair.

Ubiquitin family modifications and template switching

Homologous recombination plays an important role in the maintenance of genome integrity. Arrested forks and DNA lesions trigger strand annealing events, called template switching, which can provide for accurate damage bypass, but can also lead to chromosome rearrangements. Advances have been made in understanding the underlying mechanisms for these events and in elucidating the factors involved. Ubiquitin- and SUMO-mediated modification pathways have emerged as key players in regulating damage-induced template switching. Here I review the biological significance of template switching at the nexus of DNA replication and recombination, and the role of ubiquitin-like modifications in mediating and controlling this process.

RAF gene fusion breakpoints in pediatric brain tumors are characterized by significant enrichment of sequence microhomology

Gene fusions involving members of the RAF family of protein kinases have recently been identified as characteristic aberrations of low-grade astrocytomas, the most common tumors of the central nervous system in children. While it has been shown that these fusions cause constitutive activation of the ERK/MAPK pathway, very little is known about their formation. Here, we present a detailed analysis of RAF gene fusion breakpoints from a well-characterized cohort of 43 low-grade astrocytomas. Our findings show that the rearrangements that generate these RAF gene fusions may be simple or complex and that both inserted nucleotides and microhomology are common at the DNA breakpoints. Furthermore, we identify novel enrichment of microhomologous sequences in the regions immediately flanking the breakpoints. We thus provide evidence that the tandem duplications responsible for these fusions are generated by microhomology-mediated break-induced replication (MMBIR). Although MMBIR has previously been implicated in the pathogenesis of other diseases and the evolution of eukaryotic genomes, we demonstrate here that the proposed details of MMBIR are consistent with a recurrent rearrangement in cancer. Our analysis of repetitive elements, Z-DNA and sequence motifs in the fusion partners identified significant enrichment of the human minisatellite conserved sequence/χ-like element at one side of the breakpoint. Therefore, in addition to furthering our understanding of low-grade astrocytomas, this study provides insights into the molecular mechanistic details of MMBIR and the sequence of events that occur in the formation of genomic rearrangements.

Recurrence, submicroscopic complexity, and potential clinical relevance of copy gains detected by array CGH that are shown to be unbalanced insertions by FISH

Insertions occur when a segment of one chromosome is translocated and inserted into a new region of the same chromosome or a non-homologous chromosome. We report 71 cases with unbalanced insertions identified using array CGH and FISH in 4909 cases referred to our laboratory for array CGH and found to have copy-number abnormalities. Although the majority of insertions were non-recurrent, several recurrent unbalanced insertions were detected, including three der(Y)ins(Y;18)(q?11.2;p11.32p11.32)pat inherited from parents carrying an unbalanced insertion. The clinical significance of these recurrent rearrangements is unclear, although the small size, limited gene content, and inheritance pattern of each suggests that the phenotypic consequences may be benign. Cryptic, submicroscopic duplications were observed at or near the insertion sites in two patients, further confounding the clinical interpretation of these insertions. Using FISH, linear amplification, and array CGH, we identified a 126-kb duplicated region from 19p13.3 inserted into MECP2 at Xq28 in a patient with symptoms of Rett syndrome. Our results demonstrate that although the interpretation of most non-recurrent insertions is unclear without high-resolution insertion site characterization, the potential for an otherwise benign duplication to result in a clinically relevant outcome through the disruption of a gene necessitates the use of FISH to determine whether copy-number gains detected by array CGH represent tandem duplications or unbalanced insertions. Further follow-up testing using techniques such as linear amplification or sequencing should be used to determine gene involvement at the insertion site after FISH has identified the presence of an insertion.

Mus81 and Yen1 promote reciprocal exchange during mitotic recombination to maintain genome integrity in budding yeast

Holliday junction (HJ) resolution is required for segregation of chromosomes and for formation of crossovers during homologous recombination. The identity of the resolvase(s) that functions in vivo has yet to be established, although several proteins able to cut HJs in vitro have been identified as candidates in yeasts and mammals. Using an assay to detect unselected products of mitotic recombination, we found a significant decrease in crossovers in the Saccharomyces cerevisiae mus81Δ mutant. Yen1 serves a backup function responsible for resolving intermediates in mus81Δ mutants, or when conversion tracts are short. In the absence of both Mus81 and Yen1, intermediates are not channeled exclusively to noncrossover recombinants, but instead are processed by Pol32-dependent break-induced replication (BIR). The channeling of recombination from reciprocal exchange to BIR results in greatly increased spontaneous loss of heterozygosity (LOH) and chromosome mis-segregation in the mus81Δ yen1Δ mutant, typical of the genomic instability found in tumor cells.

A human genome structural variation sequencing resource reveals insights into mutational mechanisms

Understanding the prevailing mutational mechanisms responsible for human genome structural variation requires uniformity in the discovery of allelic variants and precision in terms of breakpoint delineation. We develop a resource based on capillary end sequencing of 13.8 million fosmid clones from 17 human genomes and characterize the complete sequence of 1054 large structural variants corresponding to 589 deletions, 384 insertions, and 81 inversions. We analyze the 2081 breakpoint junctions and infer potential mechanism of origin. Three mechanisms account for the bulk of germline structural variation: microhomology-mediated processes involving short (2-20 bp) stretches of sequence (28%), nonallelic homologous recombination (22%), and L1 retrotransposition (19%). The high quality and long-range continuity of the sequence reveals more complex mutational mechanisms, including repeat-mediated inversions and gene conversion, that are most often missed by other methods, such as comparative genomic hybridization, single nucleotide polymorphism microarrays, and next-generation sequencing.

An Sp185/333 gene cluster from the purple sea urchin and putative microsatellite-mediated gene diversification

Background

The immune system of the purple sea urchin, Strongylocentrotus purpuratus, is complex and sophisticated. An important component of sea urchin immunity is the Sp185/333 gene family, which is significantly upregulated in immunologically challenged animals. The Sp185/333 genes are less than 2 kb with two exons and are members of a large diverse family composed of greater than 40 genes. The S. purpuratus genome assembly, however, contains only six Sp185/333 genes. This underrepresentation could be due to the difficulties that large gene families present in shotgun assembly, where multiple similar genes can be collapsed into a single consensus gene.

Results

To understand the genomic organization of the Sp185/333 gene family, a BAC insert containing Sp185/333 genes was assembled, with careful attention to avoiding artifacts resulting from collapse or artificial duplication/expansion of very similar genes. Twelve candidate BAC assemblies were generated with varying parameters and the optimal assembly was identified by PCR, restriction digests, and subclone sequencing. The validated assembly contained six Sp185/333 genes that were clustered in a 34 kb region at one end of the BAC with five of the six genes tightly clustered within 20 kb. The Sp185/333 genes in this cluster were no more similar to each other than to previously sequenced Sp185/333 genes isolated from three different animals. This was unexpected given their proximity and putative effects of gene homogenization in closely linked, similar genes. All six genes displayed significant similarity including both 5' and 3' flanking regions, which were bounded by microsatellites. Three of the Sp185/333 genes and their flanking regions were tandemly duplicated such that each repeated segment consisted of a gene plus 0.7 kb 5' and 2.4 kb 3' of the gene (4.5 kb total). Both edges of the segmental duplications were bounded by different microsatellites.

Conclusions

The high sequence similarity of the Sp185/333 genes and flanking regions, suggests that the microsatellites may promote genomic instability and are involved with gene duplication and/or gene conversion and the extraordinary sequence diversity of this family.

Genome stability control by checkpoint regulation of tRNA gene transcription

The RNA polymerase III pre-initiation complex (PIC) assembled on yeast tRNA genes naturally causes replication fork pausing that contributes to genome instability. Mechanistic coupling of the fork pausing activity of tRNA genes to replication has long been considered likely, but only recently demonstrated. In contrast to the expectation that this coupling might occur by a passive mechanism such as direct disruption of transcription factor-DNA complexes by a component of the replisome, it turns out that disassembly of the RNA polymerase III PIC is actively controlled by the replication stress checkpoint signal transduction pathway. This advance supports a new model in which checkpoint-dependent disassembly of the transcription machinery at tRNA genes is a vital component of an overall system of genome stability control that also targets replication and DNA repair proteins.

The role of replication bypass pathways in dicentric chromosome formation in budding yeast

Gross chromosomal rearrangements (GCRs) are large scale changes to chromosome structure and can lead to human disease. We previously showed in Saccharomyces cerevisiae that nearby inverted repeat sequences (∼20-200 bp of homology, separated by ∼1-5 kb) frequently fuse to form unstable dicentric and acentric chromosomes. Here we analyzed inverted repeat fusion in mutants of three sets of genes. First, we show that genes in the error-free postreplication repair (PRR) pathway prevent fusion of inverted repeats, while genes in the translesion branch have no detectable role. Second, we found that siz1 mutants, which are defective for Srs2 recruitment to replication forks, and srs2 mutants had opposite effects on instability. This may reflect separate roles for Srs2 in different phases of the cell cycle. Third, we provide evidence for a faulty template switch model by studying mutants of DNA polymerases; defects in DNA pol delta (lagging strand polymerase) and Mgs1 (a pol delta interacting protein) lead to a defect in fusion events as well as allelic recombination. Pol delta and Mgs1 may collaborate either in strand annealing and/or DNA replication involved in fusion and allelic recombination events. Fourth, by studying genes implicated in suppression of GCRs in other studies, we found that inverted repeat fusion has a profile of genetic regulation distinct from these other major forms of GCR formation.

Loss of DNA replication control is a potent inducer of gene amplification

Eukaryotic cells use numerous mechanisms to ensure that no segment of their DNA is inappropriately re-replicated, but the importance of this stringent control on genome stability has not been tested. Here we show that re-replication in Saccharomyces cerevisiae can strongly induce the initial step of gene amplification, increasing gene copy number from one to two or more. The resulting amplicons consist of large internal chromosomal segments that are bounded by Ty repetitive elements and are intrachromosomally arrayed at their endogenous locus in direct head-to-tail orientation. These re-replication-induced gene amplifications are mediated by nonallelic homologous recombination between the repetitive elements. We suggest that re-replication may be a contributor to gene copy number changes, which are important in fields such as cancer biology, evolution, and human genetics.

Formation of complex and unstable chromosomal translocations in yeast

Genome instability, associated with chromosome breakage syndromes and most human cancers, is still poorly understood. In the yeast Saccharomyces cerevisiae, numerous genes with roles in the preservation of genome integrity have been identified. DNA-damage-checkpoint-deficient yeast cells that lack Sgs1, a RecQ-like DNA helicase related to the human Bloom's-syndrome-associated helicase BLM, show an increased rate of genome instability, and we have previously shown that they accumulate recurring chromosomal translocations between three similar genes, CAN1, LYP1 and ALP1. Here, the chromosomal location, copy number and sequence similarity of the translocation targets ALP1 and LYP1 were altered to gain insight into the formation of complex translocations. Among 844 clones with chromosomal rearrangements, 93 with various types of simple and complex translocations involving CAN1, LYP1 and ALP1 were identified. Breakpoint sequencing and mapping showed that the formation of complex translocation types is strictly dependent on the location of the initiating DNA break and revealed that complex translocations arise via a combination of interchromosomal translocation and template-switching, as well as from unstable dicentric intermediates. Template-switching occurred between sequences on the same chromosome, but was inhibited if the genes were transferred to different chromosomes. Unstable dicentric translocations continuously gave rise to clones with multiple translocations in various combinations, reminiscent of intratumor heterogeneity in human cancers. Base substitutions and evidence of DNA slippage near rearrangement breakpoints revealed that translocation formation can be accompanied by point mutations, and their presence in different translocation types within the same clone provides evidence that some of the different translocation types are derived from each other rather than being formed de novo. These findings provide insight into eukaryotic genome instability, especially the formation of translocations and the sources of intraclonal heterogeneity, both of which are often associated with human cancers.

Increased mutagenesis and unique mutation signature associated with mitotic gene conversion

To examine the fidelity of DNA synthesis during double-strand break (DSB) repair in Saccharomyces cerevisiae we studied gene conversion in which both strands of DNA are newly synthesized. The mutation rate increases up to 1400 times over spontaneous events, with a significantly different mutation signature. Especially prominent are microhomology-mediated template switches. Recombination-induced mutations are largely independent of mismatch repair, by DNA polymerases Polzeta, Poleta, and Pol32, but result from errors made by Poldelta and Polepsilon. These observations suggest that increased DSB frequencies in oncogene-activated mammalian cells may also increase the probability of acquiring mutations required for transition to a cancerous state.

Copy-number variations on the X chromosome in Japanese patients with mental retardation detected by array-based comparative genomic hybridization analysis

X-linked mental retardation (XLMR) is a common, clinically complex and genetically heterogeneous disease arising from many mutations along the X chromosome. Although research during the past decade has identified >90 XLMR genes, many more remain uncharacterized. In this study, copy-number variations (CNVs) were screened in individuals with MR from 144 families by array-based comparative genomic hybridization (aCGH) using a bacterial artificial chromosome-based X-tiling array. Candidate pathogenic CNVs (pCNVs) were detected in 10 families (6.9%). Five of the families had pCNVs involving known XLMR genes, duplication of Xq28 containing MECP2 in three families, duplication of Xp11.22-p11.23 containing FTSJ1 and PQBP1 in one family, and deletion of Xp11.22 bearing SHROOM4 in one family. New candidate pCNVs were detected in five families as follows: identical complex pCNVs involved in dup(X)(p22.2) and dup(X)(p21.3) containing part of REPS2, NHS and IL1RAPL1 in two unrelated families, duplication of Xp22.2 including part of FRMPD4, duplication of Xq21.1 including HDX and deletion of Xq24 noncoding region in one family, respectively. Both parents and only mother samples were available in six and three families, respectively, and pCNVs were inherited from each of their mothers in those families other than a family of the proband with deletion of SHROOM4. This study should help to identify the novel XLMR genes and mechanisms leading to MR and reveal the clinical conditions and genomic background of XLMR.

Intragenic rearrangements in NRXN1 in three families with autism spectrum disorder, developmental delay, and speech delay

NRXN1 is highly expressed in brain and has been shown recently to be associated with ASD, schizophrenia, cognitive and behavioral abnormalities, and alcohol and nicotine dependence. We present three families, in whom we identified intragenic rearrangements within NRXN1 using a clinical targeted oligonucleotide array CGH. An approximately 380 kb deletion was identified in a woman with Asperger syndrome, anxiety, and depression and in all four of her children affected with autism, anxiety, developmental delay, and speech delay but not in an unaffected child. An approximately 180 kb tandem duplication was found in a patient with autistic disorder and cognitive delays, and in his mother and younger brother who have speech delay. An approximately 330 kb tandem duplication was identified in a patient with autistic features. As predicted by conceptual translation, all three genomic rearrangements led to the premature truncation of NRXN1. Our data support previous observations that NRXN1 may be pathogenic in a wide variety of psychiatric diseases, including autism spectrum disorder, global developmental delay, anxiety, and depression.

Mechanisms for nonrecurrent genomic rearrangements associated with CMT1A or HNPP: rare CNVs as a cause for missing heritability

Genomic rearrangements involving the peripheral myelin protein gene (PMP22) in human chromosome 17p12 are associated with neuropathy: duplications cause Charcot-Marie-Tooth disease type 1A (CMT1A), whereas deletions lead to hereditary neuropathy with liability to pressure palsies (HNPP). Our previous studies showed that >99% of these rearrangements are recurrent and mediated by nonallelic homologous recombination (NAHR). Rare copy number variations (CNVs) generated by nonrecurrent rearrangements also exist in 17p12, but their underlying mechanisms are not well understood. We investigated 21 subjects with rare CNVs associated with CMT1A or HNPP by oligonucleotide-based comparative genomic hybridization microarrays and breakpoint sequence analyses, and we identified 17 unique CNVs, including two genomic deletions, ten genomic duplications, two complex rearrangements, and three small exonic deletions. Each of these CNVs includes either the entire PMP22 gene, or exon(s) only, or ultraconserved potential regulatory sequences upstream of PMP22, further supporting the contention that PMP22 is the critical gene mediating the neuropathy phenotypes associated with 17p12 rearrangements. Breakpoint sequence analysis reveals that, different from the predominant NAHR mechanism in recurrent rearrangement, various molecular mechanisms, including nonhomologous end joining, Alu-Alu-mediated recombination, and replication-based mechanisms (e.g., FoSTeS and/or MMBIR), can generate nonrecurrent 17p12 rearrangements associated with neuropathy. We document a multitude of ways in which gene function can be altered by CNVs. Given the characteristics, including small size, structural complexity, and location outside of coding regions, of selected rare CNVs, their identification remains a challenge for genome analysis. Rare CNVs may potentially represent an important portion of "missing heritability" for human diseases.

Break-induced replication requires all essential DNA replication factors except those specific for pre-RC assembly

Break-induced replication (BIR) is an efficient homologous recombination (HR) pathway employed to repair a DNA double-strand break (DSB) when homology is restricted to one end. All three major replicative DNA polymerases are required for BIR, including the otherwise nonessential Pol32 subunit. Here we show that BIR requires the replicative DNA helicase (Cdc45, the GINS, and Mcm2-7 proteins) as well as Cdt1. In contrast, both subunits of origin recognition complex (ORC) and Cdc6, which are required to create a prereplication complex (pre-RC), are dispensable. The Cdc7 kinase, required for both initiation of DNA replication and post-replication repair (PRR), is also required for BIR. Ubiquitination and sumoylation of the DNA processivity clamp PCNA play modest roles; in contrast, PCNA alleles that suppress pol32Delta's cold sensitivity fail to suppress its role in BIR, and are by themselves dominant inhibitors of BIR. These results suggest that origin-independent BIR involves cross-talk between normal DNA replication factors and PRR.

Origin and fate of pseudogenes in Hemiascomycetes: a comparative analysis

Background

Pseudogenes are ubiquitous genetic elements that derive from functional genes after mutational inactivation. Characterization of pseudogenes is important to understand genome dynamics and evolution, and its significance increases when several genomes of related organisms can be compared. Among yeasts, only the genome of the S. cerevisiae reference strain has been analyzed so far for pseudogenes.

Results

We present here the first comparative analysis of pseudogenes within the fully sequenced and annotated genomes of eight yeast species, spanning the entire phylogenetic range of Hemiascomycetes. A total of 871 pseudogenes were found, out of which mutational degradation patterns and consequences on the genetic repertoire of each species could be identified. We found that most pseudogenes in yeasts originate from mutational degradation of gene copies formed after species-specific duplications but duplications of pseudogenes themselves are also encountered. In all yeasts, except in Y. lipolytica, pseudogenes tend to cluster in subtelomeric regions where they can outnumber the number of functional genes from 3 to 16 times. Pseudogenes are generally not conserved between the yeast species studied (except in two cases), consistent with their large evolutionary distances, but tend to be conserved among S. cerevisiae strains. Reiterated pseudogenization of some genes is often observed in different lineages and may affect functions essential in S. cerevisiae, which are, therefore, lost in other species. Although a variety of functions are affected by pseudogenization, there is a bias towards functions involved in the adaptation of the yeasts to their environment, and towards genes of unknown functions.

Conclusions

Our work illustrates for the first time the formation of pseudogenes in different branches of hemiascomycetous yeasts, showing their limited conservation and how they testify for the adaptation of the yeasts functional repertoires.

Identification of uncommon recurrent Potocki-Lupski syndrome-associated duplications and the distribution of rearrangement types and mechanisms in PTLS

Nonallelic homologous recombination (NAHR) can mediate recurrent rearrangements in the human genome and cause genomic disorders. Smith-Magenis syndrome (SMS) and Potocki-Lupski syndrome (PTLS) are genomic disorders associated with a 3.7 Mb deletion and its reciprocal duplication in 17p11.2, respectively. In addition to these common recurrent rearrangements, an uncommon recurrent 5 Mb SMS-associated deletion has been identified. However, its reciprocal duplication predicted by the NAHR mechanism had not been identified. Here we report the molecular assays on 74 subjects with PTLS-associated duplications, 35 of whom are newly investigated. By both oligonucleotide-based comparative genomic hybridization and recombination hot spot analyses, we identified two cases of the predicted 5 Mb uncommon recurrent PTLS-associated duplication. Interestingly, the crossovers occur in proximity to a recently delineated allelic homologous recombination (AHR) hot spot-associated sequence motif, further documenting the common hot spot features shared between NAHR and AHR. An additional eight subjects with nonrecurrent PTLS duplications were identified. The smallest region of overlap (SRO) for all of the 74 PTLS duplications examined is narrowed to a 125 kb interval containing only RAI1, a gene recently further implicated in autism. Sequence complexities consistent with DNA replication-based mechanisms were identified in four of eight (50%) newly identified nonrecurrent PTLS duplications. Our findings of the uncommon recurrent PTLS-associated duplication at a relative prevalence reflecting the de novo mutation rate and the distribution of 17p11.2 duplication types in PTLS reveal insights into both the contributions of new mutations and the different underlying mechanisms that generate genomic rearrangements causing genomic disorders.

Differences in the DNA replication of unicellular eukaryotes and metazoans: known unknowns

Although the basic mechanisms of DNA synthesis are conserved across species, there are differences between simple and complex organisms. In contrast to lower eukaryotes, replication origins in complex eukaryotes lack DNA sequence specificity, can be activated in response to stressful conditions and require poorly conserved factors for replication firing. The response to replication fork damage is monitored by conserved proteins, such as the TIPIN-TIM-CLASPIN complex. The absence of this complex induces severe effects on yeast replication, whereas in higher eukaryotes it is only crucial when the availability of replication origins is limiting. Finally, the dependence of DNA replication on homologous recombination proteins such as RAD51 and the MRE11-RAD50-NBS1 complex is also different; they are dispensable for yeast S-phase but essential for accurate DNA replication in metazoans under unchallenged conditions. The reasons for these differences are not yet understood. Here, we focus on some of these known unknowns of DNA replication.

Alterations in DNA replication and histone levels promote histone gene amplification in Saccharomyces cerevisiae

Gene amplification, a process that increases the copy number of a gene or a genomic region to two or more, is utilized by many organisms in response to environmental stress or decreased levels of a gene product. Our previous studies in Saccharomyces cerevisiae identified the amplification of a histone H2A-H2B gene pair, HTA2-HTB2, in response to the deletion of the other H2A-H2B gene pair, HTA1-HTB1. This amplification arises from a recombination event between two flanking Ty1 elements to form a new, stable circular chromosome and occurs at a frequency higher than has been observed for other Ty1-Ty1 recombination events. To understand the regulation of this amplification event, we screened the S. cerevisiae nonessential deletion set for mutations that alter the amplification frequency. Among the deletions that increase HTA2-HTB2 amplification frequency, we identified those that either decrease DNA replication fork progression (rrm3Delta, dpb3Delta, dpb4Delta, and clb5Delta) or that reduce histone H3-H4 levels (hht2-hhf2Delta). These two classes are related because reduced histone H3-H4 levels increase replication fork pauses, and impaired replication forks cause a reduction in histone levels. Consistent with our mutant screen, we found that the introduction of DNA replication stress by hydroxyurea induces the HTA2-HTB2 amplification event. Taken together, our results suggest that either reduced histone levels or slowed replication forks stimulate the HTA2-HTB2 amplification event, contributing to the restoration of normal chromatin structure.

Whole-genome sequencing of a laboratory-evolved yeast strain

Background

Experimental evolution of microbial populations provides a unique opportunity to study evolutionary adaptation in response to controlled selective pressures. However, until recently it has been difficult to identify the precise genetic changes underlying adaptation at a genome-wide scale. New DNA sequencing technologies now allow the genome of parental and evolved strains of microorganisms to be rapidly determined.

Results

We sequenced >93.5% of the genome of a laboratory-evolved strain of the yeast Saccharomyces cerevisiae and its ancestor at >28× depth. Both single nucleotide polymorphisms and copy number amplifications were found, with specific gains over array-based methodologies previously used to analyze these genomes. Applying a segmentation algorithm to quantify structural changes, we determined the approximate genomic boundaries of a 5× gene amplification. These boundaries guided the recovery of breakpoint sequences, which provide insights into the nature of a complex genomic rearrangement.

Conclusions

This study suggests that whole-genome sequencing can provide a rapid approach to uncover the genetic basis of evolutionary adaptations, with further applications in the study of laboratory selections and mutagenesis screens. In addition, we show how single-end, short read sequencing data can provide detailed information about structural rearrangements, and generate predictions about the genomic features and processes that underlie genome plasticity.

Leaping forks at inverted repeats

Genome rearrangements are often associated with genome instability observed in cancer and other pathological disorders. Different types of repeat elements are common in genomes and are prone to instability. S-phase checkpoints, recombination, and telomere maintenance pathways have been implicated in suppressing chromosome rearrangements, but little is known about the molecular mechanisms and the chromosome intermediates generating such genome-wide instability. In the December 15, 2009, issue of Genes & Development, two studies by Paek and colleagues (2861-2875) and Mizuno and colleagues (pp. 2876-2886), demonstrate that nearby inverted repeats in budding and fission yeasts recombine spontaneously and frequently to form dicentric and acentric chromosomes. The recombination mechanism underlying this phenomenon does not appear to require double-strand break formation, and is likely caused by a replication mechanism involving template switching.

Fusion of nearby inverted repeats by a replication-based mechanism leads to formation of dicentric and acentric chromosomes that cause genome instability in budding yeast

Large-scale changes (gross chromosomal rearrangements [GCRs]) are common in genomes, and are often associated with pathological disorders. We report here that a specific pair of nearby inverted repeats in budding yeast fuse to form a dicentric chromosome intermediate, which then rearranges to form a translocation and other GCRs. We next show that fusion of nearby inverted repeats is general; we found that many nearby inverted repeats that are present in the yeast genome also fuse, as does a pair of synthetically constructed inverted repeats. Fusion occurs between inverted repeats that are separated by several kilobases of DNA and share >20 base pairs of homology. Finally, we show that fusion of inverted repeats, surprisingly, does not require genes involved in double-strand break (DSB) repair or genes involved in other repeat recombination events. We therefore propose that fusion may occur by a DSB-independent, DNA replication-based mechanism (which we term "faulty template switching"). Fusion of nearby inverted repeats to form dicentrics may be a major cause of instability in yeast and in other organisms.

Nearby inverted repeats fuse to generate acentric and dicentric palindromic chromosomes by a replication template exchange mechanism

Gene amplification plays important roles in the progression of cancer and contributes to acquired drug resistance during treatment. Amplification can initiate via dicentric palindromic chromosome production and subsequent breakage-fusion-bridge cycles. Here we show that, in fission yeast, acentric and dicentric palindromic chromosomes form by homologous recombination protein-dependent fusion of nearby inverted repeats, and that these fusions occur frequently when replication forks arrest within the inverted repeats. Genetic and molecular analyses suggest that these acentric and dicentric palindromic chromosomes arise not by previously described mechanisms, but by a replication template exchange mechanism that does not involve a DNA double-strand break. We thus propose an alternative mechanism for the generation of palindromic chromosomes dependent on replication fork arrest at closely spaced inverted repeats.

Role of double-stranded DNA translocase activity of human HLTF in replication of damaged DNA

Unrepaired DNA lesions can block the progression of the replication fork, leading to genomic instability and cancer in higher-order eukaryotes. In Saccharomyces cerevisiae, replication through DNA lesions can be mediated by translesion synthesis DNA polymerases, leading to error-free or error-prone damage bypass, or by Rad5-mediated template switching to the sister chromatid that is inherently error free. While translesion synthesis pathways are highly conserved from yeast to humans, very little is known of a Rad5-like pathway in human cells. Here we show that a human homologue of Rad5, HLTF, can facilitate fork regression and has a role in replication of damaged DNA. We found that HLTF is able to reverse model replication forks, a process which depends on its double-stranded DNA translocase activity. Furthermore, from analysis of isolated dually labeled chromosomal fibers, we demonstrate that in vivo, HLTF promotes the restart of replication forks blocked at DNA lesions. These findings suggest that HLTF can promote error-free replication of damaged DNA and support a role for HLTF in preventing mutagenesis and carcinogenesis, providing thereby for its potential tumor suppressor role.

The replication fork's five degrees of freedom, their failure and genome rearrangements

Genome rearrangements are important in pathology and evolution. The thesis of this review is that the genome is in peril when replication forks stall, and stalled forks are normally rescued by error-free mechanisms. Failure of error-free mechanisms results in large-scale chromosome changes called gross chromosomal rearrangements, GCRs, by the aficionados. In this review we discuss five error-free mechanisms a replication fork may use to overcome blockage, mechanisms that are still poorly understood. We then speculate on how genome rearrangements may occur when such mechanisms fail. Replication fork recovery failure may be an important feature of the oncogenic process. (Feedback to the authors on topics discussed herein is welcome.).

The origins and impact of primate segmental duplications

Duplicated sequences are substrates for the emergence of new genes and are an important source of genetic instability associated with rare and common diseases. Analyses of primate genomes have shown an increase in the proportion of interspersed segmental duplications (SDs) within the genomes of humans and great apes. This contrasts with other mammalian genomes that seem to have their recently duplicated sequences organized in a tandem configuration. In this review, we focus on the mechanistic origin and impact of this difference with respect to evolution, genetic diversity and primate phenotype. Although many genomes will be sequenced in the future, resolution of this aspect of genomic architecture still requires high quality sequences and detailed analyses.

Fragile genomic sites are associated with origins of replication

Genome rearrangements are mediators of evolution and disease. Such rearrangements are frequently bounded by transfer RNAs (tRNAs), transposable elements, and other repeated elements, suggesting a functional role for these elements in creating or repairing breakpoints. Though not well explored, there is evidence that origins of replication also colocalize with breakpoints. To investigate a potential correlation between breakpoints and origins, we analyzed evolutionary breakpoints defined between Saccharomyces cerevisiae and Kluyveromyces waltii and S. cerevisiae and a hypothetical ancestor of both yeasts, as well as breakpoints reported in the experimental literature. We find that origins correlate strongly with both evolutionary breakpoints and those described in the literature. Specifically, we find that origins firing earlier in S phase are more strongly correlated with breakpoints than are later-firing origins. Despite origins being located in genomic regions also bearing tRNAs and Ty elements, the correlation we observe between origins and breakpoints appears to be independent of these genomic features. This study lays the groundwork for understanding the mechanisms by which origins of replication may impact genome architecture and disease.

Gene duplication, transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.

Whirly proteins maintain plastid genome stability in Arabidopsis

Maintenance of genome stability is essential for the accurate propagation of genetic information and cell growth and survival. Organisms have therefore developed efficient strategies to prevent DNA lesions and rearrangements. Much of the information concerning these strategies has been obtained through the study of bacterial and nuclear genomes. Comparatively, little is known about how organelle genomes maintain a stable structure. Here, we report that the plastid-localized Whirly ssDNA-binding proteins are required for plastid genome stability in Arabidopsis. We show that a double KO of the genes AtWhy1 and AtWhy3 leads to the appearance of plants with variegated green/white/yellow leaves, symptomatic of nonfunctional chloroplasts. This variegation is maternally inherited, indicating defects in the plastid genome. Indeed, in all variegated lines examined, reorganized regions of plastid DNA are amplified as circular and/or head-tail concatemers. All amplified regions are delimited by short direct repeats of 10-18 bp, strongly suggesting that these regions result from illegitimate recombination between repeated sequences. This type of recombination occurs frequently in plants lacking both Whirlies, to a lesser extent in single KO plants and rarely in WT individuals. Maize mutants for the ZmWhy1 Whirly protein also show an increase in the frequency of illegitimate recombination. We propose a model where Whirlies contribute to plastid genome stability by protecting against illegitimate repeat-mediated recombination.

Mechanisms of change in gene copy number

Deletions and duplications of chromosomal segments (copy number variants, CNVs) are a major source of variation between individual humans and are an underlying factor in human evolution and in many diseases, including mental illness, developmental disorders and cancer. CNVs form at a faster rate than other types of mutation, and seem to do so by similar mechanisms in bacteria, yeast and humans. Here we review current models of the mechanisms that cause copy number variation. Non-homologous end-joining mechanisms are well known, but recent models focus on perturbation of DNA replication and replication of non-contiguous DNA segments. For example, cellular stress might induce repair of broken replication forks to switch from high-fidelity homologous recombination to non-homologous repair, thus promoting copy number change.

Mechanisms of change in gene copy number

Deletions and duplications of chromosomal segments (copy number variants, CNVs) are a major source of variation between individual humans and are an underlying factor in human evolution and in many diseases, including mental illness, developmental disorders and cancer. CNVs form at a faster rate than other types of mutation, and seem to do so by similar mechanisms in bacteria, yeast and humans. Here we review current models of the mechanisms that cause copy number variation. Non-homologous end-joining mechanisms are well known, but recent models focus on perturbation of DNA replication and replication of non-contiguous DNA segments. For example, cellular stress might induce repair of broken replication forks to switch from high-fidelity homologous recombination to non-homologous repair, thus promoting copy number change.

Complex human chromosomal and genomic rearrangements

Copy number variation (CNV) is a major source of genetic variation among humans. In addition to existing as benign polymorphisms, CNVs can also convey clinical phenotypes, including genomic disorders, sporadic diseases and complex human traits. CNV results from genomic rearrangements that can represent simple deletion or duplication of a genomic segment, or be more complex. Complex chromosomal rearrangements (CCRs) have been known for some time but their mechanisms have remained elusive. Recent technology advances and high-resolution human genome analyses have revealed that complex genomic rearrangements can account for a large fraction of non-recurrent rearrangements at a given locus. Various mechanisms, most of which are DNA-replication-based, for example fork stalling and template switching (FoSTeS) and microhomology-mediated break-induced replication (MMBIR), have been proposed for generating such complex genomic rearrangements and are probably responsible for CCR.

The DNA replication FoSTeS/MMBIR mechanism can generate genomic, genic and exonic complex rearrangements in humans

We recently proposed a DNA replication-based mechanism of fork stalling and template switching (FoSTeS) to explain the complex genomic rearrangements associated with a dysmyelinating central nervous system disorder in humans. The FoSTeS mechanism has been further generalized and molecular mechanistic details have been provided in the microhomology-mediated break-induced replication (MMBIR) model that may underlie many structural variations in genomes from all domains of life. Here we provide evidence that human genomic rearrangements ranging in size from several megabases to a few hundred base pairs can be generated by FoSTeS/MMBIR. Furthermore, we show that FoSTeS/MMBIR-mediated rearrangements can occur mitotically and can result in duplication or triplication of individual genes or even rearrangements of single exons. The FoSTeS/MMBIR mechanism can explain both the gene duplication-divergence hypothesis and exon shuffling, suggesting an important role in both genome and single-gene evolution.

Comparative genomics of protoploid Saccharomycetaceae

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

Repair-mediated duplication by capture of proximal chromosomal DNA has shaped vertebrate genome evolution

DNA double-strand breaks (DSBs) are a common form of cellular damage that can lead to cell death if not repaired promptly. Experimental systems have shown that DSB repair in eukaryotic cells is often imperfect and may result in the insertion of extra chromosomal DNA or the duplication of existing DNA at the breakpoint. These events are thought to be a source of genomic instability and human diseases, but it is unclear whether they have contributed significantly to genome evolution. Here we developed an innovative computational pipeline that takes advantage of the repetitive structure of genomes to detect repair-mediated duplication events (RDs) that occurred in the germline and created insertions of at least 50 bp of genomic DNA. Using this pipeline we identified over 1,000 probable RDs in the human genome. Of these, 824 were intra-chromosomal, closely linked duplications of up to 619 bp bearing the hallmarks of the synthesis-dependent strand-annealing repair pathway. This mechanism has duplicated hundreds of sequences predicted to be functional in the human genome, including exons, UTRs, intron splice sites and transcription factor binding sites. Dating of the duplication events using comparative genomics and experimental validation revealed that the mechanism has operated continuously but with decreasing intensity throughout primate evolution. The mechanism has produced species-specific duplications in all primate species surveyed and is contributing to genomic variation among humans. Finally, we show that RDs have also occurred, albeit at a lower frequency, in non-primate mammals and other vertebrates, indicating that this mechanism has been an important force shaping vertebrate genome evolution.

The repair of DNA double-strand breaks (DSBs) is essential for the maintenance of genome integrity. The mechanisms by which DSBs are repaired have been the subject of intense experimental investigations. It has emerged that several imperfect repair pathways exist in eukaryotes that have the potential to result in chromosomal alterations, including genomic duplications. However, it remains unclear to what extent these imperfect repair events have contributed to shaping genomes throughout evolution. Here we introduce an innovative computational approach that takes advantage of the repetitive nature of eukaryotic genomes to identify repair-mediated duplications (RD) that occurred during evolution. We discovered over one thousand RDs in the human genome, with two-thirds resulting from the capture of a chromosomal DNA segment located in close proximity to the presumed site of the DSB, giving rise to local genomic duplications. Comparative genomic analyses reveal that the mechanism has operated continuously, but with decreasing intensity during primate evolution, generating species-specific duplications in all primates surveyed and generating genomic variation among humans. Finally, we show that RDs have also occurred in non-primate mammals and other vertebrates, indicating that this is a previously under-appreciated force shaping vertebrate genomes.

Additions, losses, and rearrangements on the evolutionary route from a reconstructed ancestor to the modern Saccharomyces cerevisiae genome

Comparative genomics can be used to infer the history of genomic rearrangements that occurred during the evolution of a species. We used the principle of parsimony, applied to aligned synteny blocks from 11 yeast species, to infer the gene content and gene order that existed in the genome of an extinct ancestral yeast about 100 Mya, immediately before it underwent whole-genome duplication (WGD). The reconstructed ancestral genome contains 4,703 ordered loci on eight chromosomes. The reconstruction is complete except for the subtelomeric regions. We then inferred the series of rearrangement steps that led from this ancestor to the current Saccharomyces cerevisiae genome; relative to the ancestral genome we observe 73 inversions, 66 reciprocal translocations, and five translocations involving telomeres. Some fragile chromosomal sites were reused as evolutionary breakpoints multiple times. We identified 124 genes that have been gained by S. cerevisiae in the time since the WGD, including one that is derived from a hAT family transposon, and 88 ancestral loci at which S. cerevisiae did not retain either of the gene copies that were formed by WGD. Sites of gene gain and evolutionary breakpoints both tend to be associated with tRNA genes and, to a lesser extent, with origins of replication. Many of the gained genes in S. cerevisiae have functions associated with ethanol production, growth in hypoxic environments, or the uptake of alternative nutrient sources.

Telomere maintenance and survival in saccharomyces cerevisiae in the absence of telomerase and RAD52

Telomeres are essential features of linear genomes that are crucial for chromosome stability. Telomeric DNA is usually replenished by telomerase. Deletion of genes encoding telomerase components leads to telomere attrition with each cycle of DNA replication, eventually causing cell senescence or death. In the Saccharomyces cerevisiae strain W303, telomerase-null populations bypass senescence and, unless EXO1 is also deleted, this survival is RAD52 dependent. Unexpectedly, we found that the S. cerevisiae strain S288C could survive the removal of RAD52 and telomerase at a low frequency without additional gene deletions. These RAD52-independent survivors were propagated stably and exhibited a telomere organization typical of recombination between telomeric DNA tracts, and in diploids behaved as a multigenic trait. The polymerase-delta subunit Pol32 was dispensable for the maintenance of RAD52-independent survivors. The incidence of this rare escape was not affected by deletion of other genes necessary for RAD52-dependent survival, but correlated with initial telomere length. If W303 strains lacking telomerase and RAD52 first underwent telomere elongation, rare colonies could then bypass senescence. We suggest that longer telomeres provide a more proficient substrate for a novel telomere maintenance mechanism that does not rely on telomerase, RAD52, or POL32.

Complex rearrangements in patients with duplications of MECP2 can occur by fork stalling and template switching

Duplication at the Xq28 band including the MECP2 gene is one of the most common genomic rearrangements identified in neurodevelopmentally delayed males. Such duplications are non-recurrent and can be generated by a non-homologous end joining (NHEJ) mechanism. We investigated the potential mechanisms for MECP2 duplication and examined whether genomic architectural features may play a role in their origin using a custom designed 4-Mb tiling-path oligonucleotide array CGH assay. Each of the 30 patients analyzed showed a unique duplication varying in size from approximately 250 kb to approximately 2.6 Mb. Interestingly, in 77% of these non-recurrent duplications, the distal breakpoints grouped within a 215 kb genomic interval, located 47 kb telomeric to the MECP2 gene. The genomic architecture of this region contains both direct and inverted low-copy repeat (LCR) sequences; this same region undergoes polymorphic structural variation in the general population. Array CGH revealed complex rearrangements in eight patients; in six patients the duplication contained an embedded triplicated segment, and in the other two, stretches of non-duplicated sequences occurred within the duplicated region. Breakpoint junction sequencing was achieved in four duplications and identified an inversion in one patient, demonstrating further complexity. We propose that the presence of LCRs in the vicinity of the MECP2 gene may generate an unstable DNA structure that can induce DNA strand lesions, such as a collapsed fork, and facilitate a Fork Stalling and Template Switching event producing the complex rearrangements involving MECP2.

Replication stress induces genome-wide copy number changes in human cells that resemble polymorphic and pathogenic variants

Copy number variants (CNVs) are an important component of genomic variation in humans and other mammals. Similar de novo deletions and duplications, or copy number changes (CNCs), are now known to be a major cause of genetic and developmental disorders and to arise somatically in many cancers. A major mechanism leading to both CNVs and disease-associated CNCs is meiotic unequal crossing over, or nonallelic homologous recombination (NAHR), mediated by flanking repeated sequences or segmental duplications. Others appear to involve nonhomologous end joining (NHEJ) or aberrant replication suggesting a mitotic cell origin. Here we show that aphidicolin-induced replication stress in normal human cells leads to a high frequency of CNCs of tens to thousands of kilobases across the human genome that closely resemble CNVs and disease-associated CNCs. Most deletion and duplication breakpoint junctions were characterized by short (<6 bp) microhomologies, consistent with the hypothesis that these rearrangements were formed by NHEJ or a replication-coupled process, such as template switching. This is a previously unrecognized consequence of replication stress and suggests that replication fork stalling and subsequent error-prone repair are important mechanisms in the formation of CNVs and pathogenic CNCs in humans.