A High-Resolution View of Genome-Wide Pneumococcal Transformation

Transformation is an important mechanism of microbial evolution through which bacteria have been observed to rapidly adapt in response to clinical interventions; examples include facilitating vaccine evasion and the development of penicillin resistance in the major respiratory pathogen Streptococcus pneumoniae. To characterise the process in detail, the genomes of 124 S. pneumoniae isolates produced through in vitro transformation were sequenced and recombination events detected. Those recombinations importing the selected marker were independent of unselected events elsewhere in the genome, the positions of which were not significantly affected by local sequence similarity between donor and recipient or mismatch repair processes. However, both types of recombinations were sometimes mosaic, with multiple non-contiguous segments originating from the same molecule of donor DNA. The lengths of the unselected events were exponentially distributed with a mean of 2.3 kb, implying that recombinations are stochastically resolved with a fixed per base probability of 4.4×10⁻⁴ bp⁻¹. This distribution of recombination sizes, coupled with an observed under representation of large insertions within transferred sequence, suggests transformation has the potential to reduce the size of bacterial genomes, and is unlikely to act as an efficient mechanism for the uptake of accessory genomic loci.     

Genetic approaches for studying transgene inheritance and genetic recombination in three successive generations of transformed tobacco

Transgene integration into plant genomes is a complex process accompanied by molecular rearrangements. Classic methods that are normally used to study transgenic population genetics are generally inadequate for assessing such integration. Two major characteristics of transgenic populations are that a transgenic genome may harbor many copies of the transgene and that molecular rearrangements can create an unstable transgenic locus. In this work, we examined the segregation of T1, T2 and T3 transgenic tobacco progenies. Since transfer DNA (T-DNA) contains the NptII selectable marker gene that confers resistance to kanamycin, we used this characteristic in developing a method to estimate the number of functional inserts integrated into the genome. This approach was based on calculation of the theoretical segregation ratios in successive generations. Mendelian ratios of 3:1, 15:1 and 63:1 were confirmed for five transformation events whereas six transformation events yielded non-segregating progenies, a finding that raised questions about causal factors. A second approach based on a maximum likelihood method was performed to estimate recombination frequencies between linked inserts. Recombination estimates varied among transformation events and over generations. Some transgenic loci were unstable and evolved continuously to segregate independently in the T3 generation. Recombination and amplification of the transgene and filler DNA yielded additional transformed genotypes.     

Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.     

Homologous gene targeting in Caenorhabditis elegans by biolistic transformation

Targeted homologous recombination is a powerful approach for genome manipulation that is widely used for gene alteration and knockouts in mouse and yeast. In Caenorhabditis elegans, several methods of target-selected mutagenesis have been implemented but none of them provides the opportunity of introducing exact predefined changes into the genome. Although anecdotal cases of homologous gene targeting in C.elegans have been reported, no practical technique of gene targeting has been developed so far. In this work we demonstrate that transformation of C.elegans by microparticle bombardment (biolistic transformation) can result in homologous recombination between introduced DNA and the chromosomal locus. We describe a scaled up version of biolistic transformation that can be used as a method for homologous gene targeting in the worm.     

Suppression of Repeat-Mediated Gross Mitochondrial Genome Rearrangements by RecA in the Moss Physcomitrella patens

RecA and its ubiquitous homologs are crucial components in homologous recombination. Besides their eukaryotic nuclear counterparts, plants characteristically possess several bacterial-type RecA proteins localized to chloroplasts and/or mitochondria, but their roles are poorly understood. Here, we analyzed the role of the only mitochondrial RecA in the moss Physcomitrella patens. Disruption of the P. patens mitochondrial recA gene RECA1 caused serious defects in plant growth and development and abnormal mitochondrial morphology. Analyses of mitochondrial DNA in disruptants revealed that frequent DNA rearrangements occurred at multiple loci. Structural analysis suggests that the rearrangements, which in some cases were associated with partial deletions and amplifications of mitochondrial DNA, were due to aberrant recombination between short (<100 bp) direct and inverted repeats in which the sequences were not always identical. Such repeats are abundant in the mitochondrial genome, and interestingly many are located in group II introns. These results suggest that RECA1 does not promote but rather suppresses recombination among short repeats scattered throughout the mitochondrial genome, thereby maintaining mitochondrial genome stability. We propose that RecA-mediated homologous recombination plays a crucial role in suppression of short repeat-mediated genome rearrangements in plant mitochondria.     

Generation of Genic Diversity among Streptococcus pneumoniae Strains via Horizontal Gene Transfer during a Chronic Polyclonal Pediatric Infection

Although there is tremendous interest in understanding the evolutionary roles of horizontal gene transfer (HGT) processes that occur during chronic polyclonal infections, to date there have been few studies that directly address this topic. We have characterized multiple HGT events that most likely occurred during polyclonal infection among nasopharyngeal strains of Streptococcus pneumoniae recovered from a child suffering from chronic upper respiratory and middle-ear infections. Whole genome sequencing and comparative genomics were performed on six isolates collected during symptomatic episodes over a period of seven months. From these comparisons we determined that five of the isolates were genetically highly similar and likely represented a dominant lineage. We analyzed all genic and allelic differences among all six isolates and found that all differences tended to occur within contiguous genomic blocks, suggestive of strain evolution by homologous recombination. From these analyses we identified three strains (two of which were recovered on two different occasions) that appear to have been derived sequentially, one from the next, each by multiple recombination events. We also identified a fourth strain that contains many of the genomic segments that differentiate the three highly related strains from one another, and have hypothesized that this fourth strain may have served as a donor multiple times in the evolution of the dominant strain line. The variations among the parent, daughter, and grand-daughter recombinant strains collectively cover greater than seven percent of the genome and are grouped into 23 chromosomal clusters. While capturing in vivo HGT, these data support the distributed genome hypothesis and suggest that a single competence event in pneumococci can result in the replacement of DNA at multiple non-adjacent loci.     

Presynaptic filament dynamics in homologous recombination and DNA repair

Homologous recombination (HR) is an essential genome stability mechanism used for high-fidelity repair of DNA double-strand breaks and for the recovery of stalled or collapsed DNA replication forks. The crucial homology search and DNA strand exchange steps of HR are catalyzed by presynaptic filaments—helical filaments of a recombinase enzyme bound to single-stranded DNA (ssDNA). Presynaptic filaments are fundamentally dynamic structures, the assembly, catalytic turnover, and disassembly of which must be closely coordinated with other elements of the DNA recombination, repair, and replication machinery in order for genome maintenance functions to be effective. Here, we reviewed the major dynamic elements controlling the assembly, activity, and disassembly of presynaptic filaments; some intrinsic such as recombinase ATP-binding and hydrolytic activities, others extrinsic such as ssDNA-binding proteins, mediator proteins, and DNA motor proteins. We examined dynamic behavior on multiple levels, including atomic- and filament-level structural changes associated with ATP binding and hydrolysis as evidenced in crystal structures, as well as subunit binding and dissociation events driven by intrinsic and extrinsic factors. We examined the biochemical properties of recombination proteins from four model systems (T4 phage, Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens), demonstrating how their properties are tailored for the context-specific requirements in these diverse species. We proposed that the presynaptic filament has evolved to rely on multiple external factors for increased multilevel regulation of HR processes in genomes with greater structural and sequence complexity.     

Dialects of the DNA Uptake Sequence in Neisseriaceae

In all sexual organisms, adaptations exist that secure the safe reassortment of homologous alleles and prevent the intrusion of potentially hazardous alien DNA. Some bacteria engage in a simple form of sex known as transformation. In the human pathogen Neisseria meningitidis and in related bacterial species, transformation by exogenous DNA is regulated by the presence of a specific DNA Uptake Sequence (DUS), which is present in thousands of copies in the respective genomes. DUS affects transformation by limiting DNA uptake and recombination in favour of homologous DNA. The specific mechanisms of DUS–dependent genetic transformation have remained elusive. Bioinformatic analyses of family Neisseriaceae genomes reveal eight distinct variants of DUS. These variants are here termed DUS dialects, and their effect on interspecies commutation is demonstrated. Each of the DUS dialects is remarkably conserved within each species and is distributed consistent with a robust Neisseriaceae phylogeny based on core genome sequences. The impact of individual single nucleotide transversions in DUS on meningococcal transformation and on DNA binding and uptake is analysed. The results show that a DUS core 5′-CTG-3′ is required for transformation and that transversions in this core reduce DNA uptake more than two orders of magnitude although the level of DNA binding remains less affected. Distinct DUS dialects are efficient barriers to interspecies recombination in N. meningitidis, N. elongata, Kingella denitrificans, and Eikenella corrodens, despite the presence of the core sequence. The degree of similarity between the DUS dialect of the recipient species and the donor DNA directly correlates with the level of transformation and DNA binding and uptake. Finally, DUS–dependent transformation is documented in the genera Eikenella and Kingella for the first time. The results presented here advance our understanding of the function and evolution of DUS and genetic transformation in bacteria, and define the phylogenetic relationships within the Neisseriaceae family.     

Extensive homologous recombination between introduced and native regulatory plastid DNA elements in transplastomic plants

Homologous recombination within plastids directs plastid genome transformation for foreign gene expression and study of plastid gene function. Though transgenes are generally efficiently targeted to their desired insertion site, unintended homologous recombination events have been observed during plastid transformation. To understand the nature and abundance of these recombination events, we analyzed transplastomic tobacco lines derived from three different plastid transformation vectors utilizing two different loci for foreign gene insertion. Two unintended recombinant plastid DNA species were formed from each regulatory plastid DNA element included in the transformation vector. Some of these recombinant DNA species accumulated to as much as 10–60% of the amount of the desired integrated transgenic sequence in T0 plants. Some of the recombinant DNA species undergo further, “secondary” recombination events, resulting in an even greater number of recombinant plastid DNA species. The abundance of novel recombinant DNA species was higher in T0 plants than in T1 progeny, indicating that the ancillary recombination events described here may have the greatest impact during selection and regeneration of transformants. A line of transplastomic tobacco was identified containing an antibiotic resistance gene unlinked from the intended transgene insertion as a result of an unintended recombination event, indicating that the homologous recombination events described here may hinder efficient recovery of plastid transformants containing the desired transgene. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Gene deletion from Plasmodium falciparum using FLP and Cre recombinases: Implications for applied site-specific recombination

The ability to manipulate the genome and induce site-specific recombination using either Flippase (FLP) or Cre recombinase has been useful in many systems including Plasmodium berghei for specific deletion events or to obtain conditional gene expression. To test whether these recombinases are active in Plasmodium falciparum we constructed gene knockouts that contain sequences recognised as templates for site-specific recombination. We tested the ability of FLP and Cre recombinases, expressed conditionally in P. falciparum, to mediate deletion of the human dihydrofolate reductase (hdhfr) drug resistance gene. We show that Cre recombinase is capable of efficient removal of hdhfr by site-specific recombination. In contrast, FLP recombinase is very inefficient, even at the optimum temperature of 30 °C for this enzyme. These results demonstrate that Cre recombinase can be utilised in P. falciparum for deletion of specific sequences such as drug resistance genes. This can be exploited for recycling of drug resistance cassettes and for the design of specific recombination events in P. falciparum.     

Transfection of plant mitochondria and in organello gene integration

Investigation and manipulation of mitochondrial genetics in animal and plant cells remains restricted by the lack of an efficient in vivo transformation methodology. Mitochondrial transfection in whole cells and maintenance of the transfected DNA are main issues on this track. We showed earlier that isolated mitochondria from different organisms can import DNA. Exploiting this mechanism, we assessed the possibility to maintain exogenous DNA in plant organelles. Whereas homologous recombination is scarce in the higher plant nuclear compartment, recombination between large repeats generates the multipartite structure of the plant mitochondrial genome. These processes are under strict surveillance to avoid extensive genomic rearrangements. Nevertheless, following transfection of isolated organelles with constructs composed of a partial gfp gene flanked by fragments of mitochondrial DNA, we demonstrated in organello homologous recombination of the imported DNA with the resident DNA and integration of the reporter gene. Recombination yielded insertion of a continuous exogenous DNA fragment including the gfp sequence and at least 0.5 kb of flanking sequence on each side. According to our observations, transfection constructs carrying multiple sequences homologous to the mitochondrial DNA should be suitable and targeting of most regions in the organelle genome should be feasible, making the approach of general interest.     

The aadA Gene of Plasmid R100 Confers Resistance to Spectinomycin and Streptomycin in Myxococcus xanthus

Plasmids with the aadA gene from plasmid R100, which confers resistance to the aminoglycosides spectinomycin and streptomycin in Escherchia coli, can be introduced into wild-type Myxococcus xanthus, strain DK1622, by electroporation. Recombinant M. xanthus strains with integrated plasmids carrying the aadA gene acquire resistance to high levels of these antibiotics. Selection for aadA in M. xanthus can be carried out independently of, or simultaneously with, selection for resistance to kanamycin. The kinds and frequencies of recombination events observed between integrative plasmids with aadA and the M. xanthus chromosome are similar to those observed after the transformation of yeast. Cleavage of integrative plasmid DNA at a site adjacent to a region of homology between the plasmid and the M. xanthus genome favors the targeted disruption of M. xanthus genes by allele replacement.     

Marker rescue from the Nicotiana tabacum plastid genome using a plastid/Escherichia coli shuttle vector

We recently reported an 868-bp plastid DNA minicircle, NICE1, that formed during transformation in a transplastomic Nicotiana tabacum line. Shuttle plasmids containing NICEI sequences were maintained extrachromosomally in plastids and shown to undergo recombination with NICE1 sequences on the plastid genome. To prove the general utility of the shuttle plasmids, we tested whether plastid genes outside the NICE1 region could be rescued in Escherichia coli. The NICE1-based rescue plasmid, pNICER1, carries NICE1 sequences for maintenance in plastids, the CoIE1 ori for maintenance in E. coli and a spectinomcyin resistance gene (aadA) for selection in both systems. In addition, pNICERl carries a defective kanamycin resistance gene, kan*, to target the rescue of a functional kanamycin resistance gene, kan, from the recipient plastid genome. pNICERl was introduced into plastids where recombination could occur between the homologous kan/kan* sequences, and subsequently rescued in E. coli to recover the products of recombination. Based on the expression of kanamycin resistance in E. coli and the analysis of three restriction fragment polymorphisms, recombinant kan genes were recovered at a high frequency. Efficient rescue of kan from the plastid genome in E. coli indicates that NICE 1-based plasmids are suitable for rescuing mutations from any part of the plastid genome, expanding the repertoire of genetic tools available for plastid biology.     

Diversity of Prdm9 Zinc Finger Array in Wild Mice Unravels New Facets of the Evolutionary Turnover of this Coding Minisatellite

In humans and mice, meiotic recombination events cluster into narrow hotspots whose genomic positions are defined by the PRDM9 protein via its DNA binding domain constituted of an array of zinc fingers (ZnFs). High polymorphism and rapid divergence of the Prdm9 gene ZnF domain appear to involve positive selection at DNA-recognition amino-acid positions, but the nature of the underlying evolutionary pressures remains a puzzle. Here we explore the variability of the Prdm9 ZnF array in wild mice, and uncovered a high allelic diversity of both ZnF copy number and identity with the caracterization of 113 alleles. We analyze features of the diversity of ZnF identity which is mostly due to non-synonymous changes at codons −1, 3 and 6 of each ZnF, corresponding to amino-acids involved in DNA binding. Using methods adapted to the minisatellite structure of the ZnF array, we infer a phylogenetic tree of these alleles. We find the sister species Mus spicilegus and M. macedonicus as well as the three house mouse (Mus musculus) subspecies to be polyphyletic. However some sublineages have expanded independently in Mus musculus musculus and M. m. domesticus, the latter further showing phylogeographic substructure. Compared to random genomic regions and non-coding minisatellites, none of these patterns appears exceptional. In silico prediction of DNA binding sites for each allele, overlap of their alignments to the genome and relative coverage of the different families of interspersed repeated elements suggest a large diversity between PRDM9 variants with a potential for highly divergent distributions of recombination events in the genome with little correlation to evolutionary distance. By compiling PRDM9 ZnF protein sequences in Primates, Muridae and Equids, we find different diversity patterns among the three amino-acids most critical for the DNA-recognition function, suggesting different diversification timescales.     

The Pif1 helicase is actively inhibited during meiotic recombination which restrains gene conversion tract length

Meiotic recombination ensures proper chromosome segregation to form viable gametes and results in gene conversions events between homologs. Conversion tracts are shorter in meiosis than in mitotically dividing cells. This results at least in part from the binding of a complex, containing the Mer3 helicase and the MutLβ heterodimer, to meiotic recombination intermediates. The molecular actors inhibited by this complex are elusive. The Pif1 DNA helicase is known to stimulate DNA polymerase delta (Pol δ) -mediated DNA synthesis from D-loops, allowing long synthesis required for break-induced replication. We show that Pif1 is also recruited genome wide to meiotic DNA double-strand break (DSB) sites. We further show that Pif1, through its interaction with PCNA, is required for the long gene conversions observed in the absence of MutLβ recruitment to recombination sites. In vivo, Mer3 interacts with the PCNA clamp loader RFC, and in vitro, Mer3-MutLβ ensemble inhibits Pif1-stimulated D-loop extension by Pol δ and RFC-PCNA. Mechanistically, our results suggest that Mer3-MutLβ may compete with Pif1 for binding to RFC-PCNA. Taken together, our data show that Pif1’s activity that promotes meiotic DNA repair synthesis is restrained by the Mer3-MutLβ ensemble which in turn prevents long gene conversion tracts and possibly associated mutagenesis.     

Genetic Analysis of Variation in Human Meiotic Recombination

The number of recombination events per meiosis varies extensively among individuals. This recombination phenotype differs between female and male, and also among individuals of each gender. In this study, we used high-density SNP genotypes of over 2,300 individuals and their offspring in two datasets to characterize recombination landscape and to map the genetic variants that contribute to variation in recombination phenotypes. We found six genetic loci that are associated with recombination phenotypes. Two of these (RNF212 and an inversion on chromosome 17q21.31) were previously reported in the Icelandic population, and this is the first replication in any other population. Of the four newly identified loci (KIAA1462, PDZK1, UGCG, NUB1), results from expression studies provide support for their roles in meiosis. Each of the variants that we identified explains only a small fraction of the individual variation in recombination. Notably, we found different sequence variants associated with female and male recombination phenotypes, suggesting that they are regulated by different genes. Characterization of genetic variants that influence natural variation in meiotic recombination will lead to a better understanding of normal meiotic events as well as of non-disjunction, the primary cause of pregnancy loss.     

Novel Conserved Genotypes Correspond to Antibiotic Resistance Phenotypes of E. coli Clinical Isolates

Current efforts to understand antibiotic resistance on the whole genome scale tend to focus on known genes even as high throughput sequencing strategies uncover novel mechanisms. To identify genomic variations associated with antibiotic resistance, we employed a modified genome-wide association study; we sequenced genomic DNA from pools of E. coli clinical isolates with similar antibiotic resistance phenotypes using SOLiD technology to uncover single nucleotide polymorphisms (SNPs) unanimously conserved in each pool. The multidrug-resistant pools were genotypically similar to SMS-3-5, a previously sequenced multidrug-resistant isolate from a polluted environment. The similarity was evenly spread across the entire genome and not limited to plasmid or pathogenicity island loci. Among the pools of clinical isolates, genomic variation was concentrated adjacent to previously reported inversion and duplication differences between the SMS-3-5 isolate and the drug-susceptible laboratory strain, DH10B. SNPs that result in non-synonymous changes in gyrA (encoding the well-known S83L allele associated with fluoroquinolone resistance), mutM, ligB, and recG were unanimously conserved in every fluoroquinolone-resistant pool. Alleles of the latter three genes are tightly linked among most sequenced E. coli genomes, and had not been implicated in antibiotic resistance previously. The changes in these genes map to amino acid positions in alpha helices that are involved in DNA binding. Plasmid-encoded complementation of null strains with either allelic variant of mutM or ligB resulted in variable responses to ultraviolet light or hydrogen peroxide treatment as markers of induced DNA damage, indicating their importance in DNA metabolism and revealing a potential mechanism for fluoroquinolone resistance. Our approach uncovered evidence that additional DNA binding enzymes may contribute to fluoroquinolone resistance and further implicate environmental bacteria as a reservoir for antibiotic resistance.     

Bacterial transformation: ComFA is a DNA‐dependent ATPase that forms complexes with ComFC and DprA

Pneumococcal natural transformation contributes to genomic plasticity, antibiotic resistance development and vaccine escape. Streptococcus pneumoniae, like many other naturally transformable species, has evolved sophisticated protein machinery for the binding and uptake of DNA. Two proteins encoded by the comF operon, ComFA and ComFC, are involved in transformation but their exact molecular roles remain unknown. In this study, we provide experimental evidence that ComFA binds to single stranded DNA (ssDNA) and has ssDNA‐dependent ATPase activity. We show that both ComFA and ComFC are essential for the transformation process in pneumococci. Moreover, we show that these proteins interact with each other and with other proteins involved in homologous recombination, such as DprA, thus placing the ComFA‐ComFC duo at the interface between DNA uptake and DNA recombination during transformation.     

Whole genome sequencing of penicillin-resistant Streptococcus pneumoniae reveals mutations in penicillin-binding proteins and in a putative iron permease

Background

Penicillin resistance in Streptococcus pneumoniae is mediated by a mosaic of genes encoding altered penicillin-binding proteins (PBPs). Nonetheless, S. pneumoniae has also developed non-PBP mechanisms implicated in penicillin resistance. In this study, whole genome sequencing of resistant organisms was used to discover mutations implicated in resistance to penicillin.

Results

We sequenced two S. pneumoniae isolates selected for resistance to penicillin in vitro. The analysis of the genome assemblies revealed that six genes were mutated in both mutants. These included three pbp genes, and three non-pbp genes, including a putative iron permease, spr1178. The nonsense mutation in spr1178 always occurred in the first step of the selection process. Although the mutants had increased resistance to penicillin, the introduction of altered versions of PBPs into a penicillin-susceptible strain by sequential transformation led to strains with a minimal increase in resistance, thus implicating other genes in resistance. The introduction by transformation of the non-PBP recurrent mutations did not increase penicillin resistance, but the introduction of the nonsense mutation in the putative iron permease spr1178 led to a reduced accumulation of reactive oxygen species following exposure to penicillin and to other bactericidal antibiotics as well.

Conclusions

This study indicates that the selection of resistance to penicillin in S. pneumoniae involves the acquisition of mutations conferring tolerance to the antibiotic-induced accumulation of oxidants, which translates into an increased survival that putatively enables the selection of major resistance determinants such as mutations in PBPs.     

The Complete Moss Mitochondrial Genome in the Angiosperm Amborella Is a Chimera Derived from Two Moss Whole-Genome Transfers

Sequencing of the 4-Mb mitochondrial genome of the angiosperm Amborella trichopoda has shown that it contains unprecedented amounts of foreign mitochondrial DNA, including four blocks of sequences that together correspond almost perfectly to one entire moss mitochondrial genome. This implies whole-genome transfer from a single moss donor but conflicts with phylogenetic results from an earlier, PCR-based study that suggested three different moss donors to Amborella. To resolve this conflict, we conducted an expanded set of phylogenetic analyses with respect to both moss lineages and mitochondrial loci. The moss DNA in Amborella was consistently placed in either of two positions, depending on the locus analyzed, as sister to the Ptychomniales or within the Hookeriales. This agrees with two of the three previously suggested donors, whereas the third is no longer supported. These results, combined with synteny analyses and other considerations, lead us to favor a model involving two successive moss-to-Amborella whole-genome transfers, followed by recombination that produced a single intact and chimeric moss mitochondrial genome integrated in the Amborella mitochondrial genome. Eight subsequent recombination events account for the state of fragmentation, rearrangement, duplication, and deletion of this chimeric moss mitochondrial genome as it currently exists in Amborella. Five of these events are associated with short-to-intermediate sized repeats. Two of the five probably occurred by reciprocal homologous recombination, whereas the other three probably occurred in a non-reciprocal manner via microhomology-mediated break-induced replication (MMBIR). These findings reinforce and extend recent evidence for an important role of MMBIR in plant mitochondrial DNA evolution.