A reversible mutation in a genomic hotspot saves bacterial swarms from extinction.

Microbial adaptation to changing environmental conditions is frequently mediated by hypermutable sequences. Here we demonstrate that such a hypermutable hotspot within a gene encoding a flagellar unit of Paenibacillus glucanolyticus generated spontaneous non-swarming mutants with increased stress resistance. These mutants, which survived conditions that eliminated wild-type cultures, could be carried by their swarming siblings when the colony spread, consequently increasing their numbers at the spreading edge. Of interest, the hypermutable nature of the aforementioned sequence enabled the non-swarming mutants to serve as "seeds" for a new generation of wild-type cells through reversion of the mutation. Using a mathematical model, we examined the survival dynamics of P. glucanolyticus colonies under fluctuating environments. Our experimental and theoretical results suggest that the non-swarming, stress-resistant mutants can save the colony from extinction. Notably, we identified this hypermutable sequence in flagellar genes of additional Paenibacillus species, suggesting that this phenomenon could be wide-spread and ecologically important.

An Approximate Bayesian Computation Approach for Modeling Genome Rearrangements.

The inference of genome rearrangement events has been extensively studied, as they play a major role in molecular evolution. However, probabilistic evolutionary models that explicitly imitate the evolutionary dynamics of such events, as well as methods to infer model parameters, are yet to be fully utilized. Here, we developed a probabilistic approach to infer genome rearrangement rate parameters using an Approximate Bayesian Computation (ABC) framework. We developed two genome rearrangement models, a basic model, which accounts for genomic changes in gene order, and a more sophisticated one which also accounts for changes in chromosome number. We characterized the ABC inference accuracy using simulations and applied our methodology to both prokaryotic and eukaryotic empirical datasets. Knowledge of genome-rearrangement rates can help elucidate their role in evolution as well as help simulate genomes with evolutionary dynamics that reflect empirical genomes.

Genomic structural plasticity of rodent-associated Bartonella in nature.

Rodent-associated Bartonella species have shown a remarkable genetic diversity and pathogenic potential. To further explore the extent of the natural intraspecific genomic variation and its potential role as an evolutionary driver, we focused on a single genetically diverse Bartonella species, Bartonella krasnovii, which circulates among gerbils and their associated fleas. Twenty genomes from 16 different B. krasnovii genotypes were fully characterized through a genome sequencing assay (using short and long read sequencing), pulse field gel electrophoresis (PFGE), and PCR validation. Genomic analyses were performed in comparison to the B. krasnovii strain OE 1-1. While, single nucleotide polymorphism represented only a 0.3% of the genome variation, structural diversity was identified in these genomes, with an average of 51 ± 24 structural variation (SV) events per genome. Interestingly, a large proportion of the SVs (>40%) was associated with prophages. Further analyses revealed that most of the SVs, and prophage insertions were found at the chromosome replication termination site (ter), suggesting this site as a plastic zone of the B. krasnovii chromosome. Accordingly, six genomes were found to be unbalanced, and essential genes near the ter showed a shift between the leading and lagging strands, revealing the SV effect on these genomes. In summary, our findings demonstrate the extensive genomic diversity harbored by wild B. krasnovii strains and suggests that its diversification is initially promoted by structural changes, probably driven by phages. These events may constantly feed the system with novel genotypes that ultimately lead to inter- and intraspecies competition and adaptation.

Protein innovation through template switching in the Saccharomyces cerevisiae lineage.

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

Limited DNA Repair Gene Repertoire in Ascomycete Yeast Revealed by Comparative Genomics.

Ascomycota is the largest phylogenetic group of fungi that includes species important to human health and wellbeing. DNA repair is important for fungal survival and genome evolution. Here, we describe a detailed comparative genomic analysis of DNA repair genes in Ascomycota. We determined the DNA repair gene repertoire in Taphrinomycotina, Saccharomycotina, Leotiomycetes, Sordariomycetes, Dothideomycetes, and Eurotiomycetes. The subphyla of yeasts, Saccharomycotina and Taphrinomycotina, have a smaller DNA repair gene repertoire comparing to Pezizomycotina. Some genes were absent from most, if not all, yeast species. To study the conservation of these genes in Pezizomycotina, we used the Gain Loss Mapping Engine algorithm that provides the expectations of gain or loss of genes given the tree topology. Genes that were absent from most of the species of Taphrinomycotina or Saccharomycotina showed lower conservation in Pezizomycotina. This suggests that the absence of some DNA repair in yeasts is not random; genes with a tendency to be lost in other classes are missing. We ranked the conservation of DNA repair genes in Ascomycota. We found that Rad51 and its paralogs were less conserved than other recombinational proteins, suggesting that there is a redundancy between Rad51 and its paralogs, at least in some species. Finally, based on the repertoire of UV repair genes, we found conditions that differentially kill the wine pathogen Brettanomyces bruxellensis and not Saccharomyces cerevisiae. In summary, our analysis provides testable hypotheses to the role of DNA repair proteins in the genome evolution of Ascomycota.

Alternative Functional rad21 Paralogs in Fusarium oxysporum.

Cohesin, the sister chromatid cohesion complex, is an essential complex that ensures faithful sister chromatid segregation in eukaryotes. It also participates in DNA repair, transcription and maintenance of chromosome structure. Mitotic cohesin is composed of Smc1, Smc3, Scc3, and Rad21/Mcd1. The meiotic cohesin complex contains Rec8, a Rad21 paralog and not Rad21 itself. Very little is known about sister chromatid cohesion in fungal plant pathogens. Fusarium oxysporum is an important fungal plant pathogen without known sexual life cycle. Here, we describe that F. oxysporum encodes for three Rad21 paralogs; Rad21, Rec8, and the first alternative Rad21 paralog in the phylum of ascomycete. This last paralog is found only in several fungal plant pathogens from the Fusarium family and thus termed rad21nc (non-conserved). Conserved rad21 (rad21c), rad21nc, and rec8 genes are expressed in F. oxysporum although the expression of rad21c is much higher than the other paralogs. F. oxysporum strains deleted for the rad21nc or rec8 genes were analyzed for their role in fungal life cycle. Δrad21nc and Δrec8 single mutants were proficient in sporulation, conidia germination, hyphal growth and pathogenicity under optimal growth conditions. Interestingly, Δrad21nc and Δrec8 single mutants germinate less effectively than wild type (WT) strains under DNA replication and mitosis stresses. We provide here the first genetic analysis of alternative rad21nc and rec8 paralogs in filamentous fungi. Our results suggest that rad21nc and rec8 may have a unique role in cell cycle related functions of F. oxysporum.

Mosaic mitochondrial-plastid insertions into the nuclear genome show evidence of both non-homologous end joining and homologous recombination.

BACKGROUND: Mitochondrial and plastid DNA fragments are continuously transferred into eukaryotic nuclear genomes, giving rise to nuclear copies of mitochondrial DNA (numts) and nuclear copies of plastid DNA (nupts). Numts and nupts are classified as simple if they are composed of a single organelle fragment or as complex if they are composed of multiple fragments. Mosaic insertions are complex insertions composed of fragments of both mitochondrial and plastid DNA. Simple numts and nupts in eukaryotes have been extensively studied, their mechanism of insertion involves non-homologous end joining (NHEJ). Mosaic insertions have been less well-studied and their mechanisms of integration are unknown. RESULTS: Here we estimated the number of nuclear mosaic insertions (numins) in nine plant genomes. We show that numins compose up to 10% of the total nuclear insertions of organelle DNA in these plant genomes. The NHEJ hallmarks typical for numts and nupts were also identified in mosaic insertions. However, the number of identified insertions that integrated via NHEJ mechanism is underestimated, as NHEJ signatures are conserved only in recent insertions and mutationally eroded in older ones. A few complex insertions show signatures of long homology that cannot be attributed to NHEJ, a novel observation that implicates gene conversion or single strand annealing mechanisms in organelle nuclear insertions. CONCLUSIONS: The common NHEJ signature that was identified here reveals that, in plant cells, mitochondria and plastid fragments in numins must meet during or prior to integration into the nuclear genome.

Failure to Recover Major Events of Gene Flux in Real Biological Data Due to Method Misapplication.

In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation, and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional models designed to investigate gene family size evolution (Count) to support claims that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here, we show that the methodology used in analyses favoring continuous gene transfers to archaea was misapplied in other studies and does not recover known events of single simultaneous origin for many genes followed by differential loss in real data: plastid genomes. Using the same software and the same settings, we reanalyzed presence/absence pattern data for proteins encoded in plastid genomes and for eukaryotic protein families acquired from plastids. Contrary to expectations under a plastid origin model, we found that the methodology employed inferred that gene acquisitions occurred uniformly across the plant tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, the methodology that recovered gradual and continuous lateral gene transfer among lineages for archaea obtains the same result for plastids, even though it is known that massive gains followed by gradual differential loss is the true evolutionary process that generated plastid gene distribution data. Our findings caution against the use of models designed to study gene family size evolution for investigating gene transfer processes, especially when transfers involving more than one gene per event are possible.

The Prevalence and Evolutionary Conservation of Inverted Repeats in Proteobacteria.

Perfect short inverted repeats (IRs) are known to be enriched in a variety of bacterial and eukaryotic genomes. Currently, it is unclear whether perfect IRs are conserved over evolutionary time scales. In this study, we aimed to characterize the prevalence and evolutionary conservation of IRs across 20 proteobacterial strains. We first identified IRs in Escherichia coli K-12 substr MG1655 and showed that they are overabundant. We next aimed to test whether this overabundance is reflected in the conservation of IRs over evolutionary time scales. To this end, for each perfect IR identified in E. coli MG1655, we collected orthologous sequences from related proteobacterial genomes. We next quantified the evolutionary conservation of these IRs, that is, the presence of the exact same IR across orthologous regions. We observed high conservation of perfect IRs: out of the 234 examined orthologous regions, 145 were more conserved than expected, which is statistically significant even after correcting for multiple testing. Our results together with previous experimental findings support a model in which imperfect IRs are corrected to perfect IRs in a preferential manner via a template switching mechanism.

Whole Genome Sequence Analysis of Mutations Accumulated in rad27Δ Yeast Strains with Defects in the Processing of Okazaki Fragments Indicates Template-Switching Events.

Okazaki fragments that are formed during lagging strand DNA synthesis include an initiating primer consisting of both RNA and DNA. The RNA fragment must be removed before the fragments are joined. In Saccharomyces cerevisiae, a key player in this process is the structure-specific flap endonuclease, Rad27p (human homolog FEN1). To obtain a genomic view of the mutational consequence of loss of RAD27, a S. cerevisiae rad27Δ strain was subcultured for 25 generations and sequenced using Illumina paired-end sequencing. Out of the 455 changes observed in 10 colonies isolated the two most common types of events were insertions or deletions (INDELs) in simple sequence repeats (SSRs) and INDELs mediated by short direct repeats. Surprisingly, we also detected a previously neglected class of 21 template-switching events. These events were presumably generated by quasi-palindrome to palindrome correction, as well as palindrome elongation. The formation of these events is best explained by folding back of the stalled nascent strand and resumption of DNA synthesis using the same nascent strand as a template. Evidence of quasi-palindrome to palindrome correction that could be generated by template switching appears also in yeast genome evolution. Out of the 455 events, 55 events appeared in multiple isolates; further analysis indicates that these loci are mutational hotspots. Since Rad27 acts on the lagging strand when the leading strand should not contain any gaps, we propose a mechanism favoring intramolecular strand switching over an intermolecular mechanism. We note that our results open new ways of understanding template switching that occurs during genome instability and evolution.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health.

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.

Endosymbiotic origin and differential loss of eukaryotic genes.

Chloroplasts arose from cyanobacteria, mitochondria arose from proteobacteria. Both organelles have conserved their prokaryotic biochemistry, but their genomes are reduced, and most organelle proteins are encoded in the nucleus. Endosymbiotic theory posits that bacterial genes in eukaryotic genomes entered the eukaryotic lineage via organelle ancestors. It predicts episodic influx of prokaryotic genes into the eukaryotic lineage, with acquisition corresponding to endosymbiotic events. Eukaryotic genome sequences, however, increasingly implicate lateral gene transfer, both from prokaryotes to eukaryotes and among eukaryotes, as a source of gene content variation in eukaryotic genomes, which predicts continuous, lineage-specific acquisition of prokaryotic genes in divergent eukaryotic groups. Here we discriminate between these two alternatives by clustering and phylogenetic analysis of eukaryotic gene families having prokaryotic homologues. Our results indicate (1) that gene transfer from bacteria to eukaryotes is episodic, as revealed by gene distributions, and coincides with major evolutionary transitions at the origin of chloroplasts and mitochondria; (2) that gene inheritance in eukaryotes is vertical, as revealed by extensive topological comparison, sparse gene distributions stemming from differential loss; and (3) that continuous, lineage-specific lateral gene transfer, although it sometimes occurs, does not contribute to long-term gene content evolution in eukaryotic genomes.

Endosymbiotic gene transfer from prokaryotic pangenomes: Inherited chimerism in eukaryotes.

Endosymbiotic theory in eukaryotic-cell evolution rests upon a foundation of three cornerstone partners--the plastid (a cyanobacterium), the mitochondrion (a proteobacterium), and its host (an archaeon)--and carries a corollary that, over time, the majority of genes once present in the organelle genomes were relinquished to the chromosomes of the host (endosymbiotic gene transfer). However, notwithstanding eukaryote-specific gene inventions, single-gene phylogenies have never traced eukaryotic genes to three single prokaryotic sources, an issue that hinges crucially upon factors influencing phylogenetic inference. In the age of genomes, single-gene trees, once used to test the predictions of endosymbiotic theory, now spawn new theories that stand to eventually replace endosymbiotic theory with descriptive, gene tree-based variants featuring supernumerary symbionts: prokaryotic partners distinct from the cornerstone trio and whose existence is inferred solely from single-gene trees. We reason that the endosymbiotic ancestors of mitochondria and chloroplasts brought into the eukaryotic--and plant and algal--lineage a genome-sized sample of genes from the proteobacterial and cyanobacterial pangenomes of their respective day and that, even if molecular phylogeny were artifact-free, sampling prokaryotic pangenomes through endosymbiotic gene transfer would lead to inherited chimerism. Recombination in prokaryotes (transduction, conjugation, transformation) differs from recombination in eukaryotes (sex). Prokaryotic recombination leads to pangenomes, and eukaryotic recombination leads to vertical inheritance. Viewed from the perspective of endosymbiotic theory, the critical transition at the eukaryote origin that allowed escape from Muller's ratchet--the origin of eukaryotic recombination, or sex--might have required surprisingly little evolutionary innovation.

Nonrandom distribution of interhomolog recombination events induced by breakage of a dicentric chromosome in Saccharomyces cerevisiae.

Dicentric chromosomes undergo breakage in mitosis, resulting in chromosome deletions, duplications, and translocations. In this study, we map chromosome break sites of dicentrics in Saccharomyces cerevisiae by a mitotic recombination assay. The assay uses a diploid strain in which one homolog has a conditional centromere in addition to a wild-type centromere, and the other homolog has only the wild-type centromere; the conditional centromere is inactive when cells are grown in galactose and is activated when the cells are switched to glucose. In addition, the two homologs are distinguishable by multiple single-nucleotide polymorphisms (SNPs). Under conditions in which the conditional centromere is activated, the functionally dicentric chromosome undergoes double-stranded DNA breaks (DSBs) that can be repaired by mitotic recombination with the homolog. Such recombination events often lead to loss of heterozygosity (LOH) of SNPs that are centromere distal to the crossover. Using a PCR-based assay, we determined the position of LOH in multiple independent recombination events to a resolution of ∼4 kb. This analysis shows that dicentric chromosomes have recombination breakpoints that are broadly distributed between the two centromeres, although there is a clustering of breakpoints within 10 kb of the conditional centromere.

High-resolution genome-wide analysis of irradiated (UV and γ-rays) diploid yeast cells reveals a high frequency of genomic loss of heterozygosity (LOH) events.

In diploid eukaryotes, repair of double-stranded DNA breaks by homologous recombination often leads to loss of heterozygosity (LOH). Most previous studies of mitotic recombination in Saccharomyces cerevisiae have focused on a single chromosome or a single region of one chromosome at which LOH events can be selected. In this study, we used two techniques (single-nucleotide polymorphism microarrays and high-throughput DNA sequencing) to examine genome-wide LOH in a diploid yeast strain at a resolution averaging 1 kb. We examined both selected LOH events on chromosome V and unselected events throughout the genome in untreated cells and in cells treated with either γ-radiation or ultraviolet (UV) radiation. Our analysis shows the following: (1) spontaneous and damage-induced mitotic gene conversion tracts are more than three times larger than meiotic conversion tracts, and conversion tracts associated with crossovers are usually longer and more complex than those unassociated with crossovers; (2) most of the crossovers and conversions reflect the repair of two sister chromatids broken at the same position; and (3) both UV and γ-radiation efficiently induce LOH at doses of radiation that cause no significant loss of viability. Using high-throughput DNA sequencing, we also detected new mutations induced by γ-rays and UV. To our knowledge, our study represents the first high-resolution genome-wide analysis of DNA damage-induced LOH events performed in any eukaryote.

Directed networks reveal genomic barriers and DNA repair bypasses to lateral gene transfer among prokaryotes.

Lateral gene transfer (LGT) plays a major role in prokaryote evolution with only a few genes that are resistant to it; yet the nature and magnitude of barriers to lateral transfer are still debated. Here, we implement directed networks to investigate donor-recipient events of recent lateral gene transfer among 657 sequenced prokaryote genomes. For 2,129,548 genes investigated, we detected 446,854 recent lateral gene transfer events through nucleotide pattern analysis. Among these, donor-recipient relationships could be specified through phylogenetic reconstruction for 7% of the pairs, yielding 32,028 polarized recent gene acquisition events, which constitute the edges of our directed networks. We find that the frequency of recent LGT is linearly correlated both with genome sequence similarity and with proteome similarity of donor-recipient pairs. Genome sequence similarity accounts for 25% of the variation in gene-transfer frequency, with proteome similarity adding only 1% to the variability explained. The range of donor-recipient GC content similarity within the network is extremely narrow, with 86% of the LGTs occurring between donor-recipient pairs having ≤5% difference in GC content. Hence, genome sequence similarity and GC content similarity are strong barriers to LGT in prokaryotes. But they are not insurmountable, as we detected 1530 recent transfers between distantly related genomes. The directed network revealed that recipient genomes of distant transfers encode proteins of nonhomologous end-joining (NHEJ; a DNA repair mechanism) far more frequently than the recipient lacking that mechanism. This implicates NHEJ in genes spread across distantly related prokaryotes through bypassing the donor-recipient sequence similarity barrier.

Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes.

The natural transfer of DNA from mitochondria to the nucleus generates nuclear copies of mitochondrial DNA (numts) and is an ongoing evolutionary process, as genome sequences attest. In humans, five different numts cause genetic disease and a dozen human loci are polymorphic for the presence of numts, underscoring the rapid rate at which mitochondrial sequences reach the nucleus over evolutionary time. In the laboratory and in nature, numts enter the nuclear DNA via non-homolgous end joining (NHEJ) at double-strand breaks (DSBs). The frequency of numt insertions among 85 sequenced eukaryotic genomes reveal that numt content is strongly correlated with genome size, suggesting that the numt insertion rate might be limited by DSB frequency. Polymorphic numts in humans link maternally inherited mitochondrial genotypes to nuclear DNA haplotypes during the past, offering new opportunities to associate nuclear markers with mitochondrial markers back in time.

Linking big: the continuing promise of evolutionary synthesis.

Synthetic science promises an unparalleled ability to find new meaning in old data, extant results, or previously unconnected methods and concepts, but pursuing synthesis can be a difficult and risky endeavor. Our experience as biologists, informaticians, and educators at the National Evolutionary Synthesis Center has affirmed that synthesis can yield major insights, but also revealed that technological hurdles, prevailing academic culture, and general confusion about the nature of synthesis can hamper its progress. By presenting our view of what synthesis is, why it will continue to drive progress in evolutionary biology, and how to remove barriers to its progress, we provide a map to a future in which all scientists can engage productively in synthetic research.

Mitochondrial insertions into primate nuclear genomes suggest the use of numts as a tool for phylogeny.

Homoplasy-free characters are a valuable and highly desired tool for molecular systematics. Nuclear sequences of mitochondrial origin (numts) are fragments of mitochondrial DNA that have been transferred into the nuclear genome. numts are passively captured into genomes and have no transposition activity, which suggests they may have utility as phylogenetic markers. Here, five fully sequenced primate genomes (human, chimpanzee, orangutan, rhesus macaque, and marmoset) are used to reconstruct the evolutionary dynamics of recent numt accumulation in a phylogenetic context. The status of 367 numt loci is used as categorical data, and a maximum parsimony approach is used to trace numt insertions on different branches of the taxonomically undisputed primate phylogenetic tree. The presence of a given numt in related taxa implies orthologous integration, whereas the absence of a numt indicates the plesiomorphic condition prior to integration. An average rate of 5.65 numts per 1 My is estimated on the tree, but insertion rates vary significantly on different branches. Two instances in which the presence-absence pattern of numts does not agree with the phylogenetic tree were identified. These events may be the result of either lineage sorting or reversal. Using the numts reported here to reconstruct primate phylogeny produces the canonical primate tree topology with high bootstrap support. Moreover, numts identified in gorilla Supercontigs were used to test the human-chimp-gorilla trichotomy, yielding a high level of support for the sister relationship of human and chimpanzee. These analyses suggest that numts are valuable phylogenetic markers that can be used for molecular systematics. It remains to be tested whether numts are useful at deeper phylogenetic levels.

Numt-mediated double-strand break repair mitigates deletions during primate genome evolution.

Non-homologous end joining (NHEJ) is the major mechanism of double-strand break repair (DSBR) in mammalian cells. NHEJ has traditionally been inferred from experimental systems involving induced double strand breaks (DSBs). Whether or not the spectrum of repair events observed in experimental NHEJ reflects the repair of natural breaks by NHEJ during chromosomal evolution is an unresolved issue. In primate phylogeny, nuclear DNA sequences of mitochondrial origin, numts, are inserted into naturally occurring chromosomal breaks via NHEJ. Thus, numt integration sites harbor evidence for the mechanisms that act on the genome over evolutionary timescales. We have identified 35 and 55 lineage-specific numts in the human and chimpanzee genomes, respectively, using the rhesus monkey genome as an outgroup. One hundred and fifty two numt-chromosome fusion points were classified based on their repair patterns. Repair involving microhomology and repair leading to nucleotide additions were detected. These repair patterns are within the experimentally determined spectrum of classical NHEJ, suggesting that information from experimental systems is representative of broader genetic loci and end configurations. However, in incompatible DSBR events, small deletions always occur, whereas in 54% of numt integration events examined, no deletions were detected. Numts show a statistically significant reduction in deletion frequency, even in comparison to DSBR involving filler DNA. Therefore, numts show a unique mechanism of integration via NHEJ. Since the deletion frequency during numt insertion is low, native overhangs of chromosome breaks are preserved, allowing us to determine that 24% of the analyzed breaks are cohesive with overhangs of up to 11 bases. These data represent, to the best of our knowledge, the most comprehensive description of the structure of naturally occurring DSBs. We suggest a model in which the sealing of DSBs by numts, and probably by other filler DNA, prevents nuclear processing of DSBs that could result in deleterious repair.