Molecular keys to speciation: DNA polymorphism and the control of genetic exchange in enterobacteria

Genomic divergence across the tree of life

Nucleotide sequence data are being harnessed to identify species, even in cases in which organisms themselves are neither in hand nor witnessed. But how genome-wide sequence divergence maps to species status is far from clear. While gene sequence divergence is commonly used to delineate bacterial species, its correspondence to established species boundaries has yet to be explored across eukaryotic taxa. Because the processes underlying gene flow differ fundamentally between prokaryotes and eukaryotes, these domains are likely to differ in the relationship between reproductive isolation and genome-wide sequence divergence. In prokaryotes, homologous recombination, the basis of gene flow, depends directly on the degree of genomic sequence divergence, whereas in sexually reproducing eukaryotes, reproductive incompatibility can stem from changes in very few genes. Guided by measures of genome-wide sequence divergence in bacteria, we gauge how genomic criteria correspond to species boundaries in eukaryotes. In recognized species of eukaryotes, levels of gene sequence divergence within species are typically very small, averaging <1% across protein-coding regions in most animals, plants, and fungi. There are even instances in which divergence between sister species is the same or less than that among conspecifics. In contrast, bacterial species, defined as populations exchanging homologous genes, show levels of divergence both within and between species that are considerably higher. Although no single threshold delineates species, eukaryotic populations with >1% genome-wide sequence divergence are likely separate species, whereas prokaryotic populations with 1% divergence are still able to recombine and thus can be considered the same species.

Investigating Additive and Replacing Horizontal Gene Transfers Using Phylogenies and Whole Genomes

Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing nonhomologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between "additive" and "replacing" HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT.

Capsules and their traits shape phage susceptibility and plasmid conjugation efficiency

Bacterial evolution is affected by mobile genetic elements like phages and conjugative plasmids, offering new adaptive traits while incurring fitness costs. Their infection is affected by the bacterial capsule. Yet, its importance has been difficult to quantify because of the high diversity of confounding mechanisms in bacterial genomes such as anti-viral systems and surface receptor modifications. Swapping capsule loci between Klebsiella pneumoniae strains allowed us to quantify their impact on plasmid and phage infection independently of genetic background. Capsule swaps systematically invert phage susceptibility, revealing serotypes as key determinants of phage infection. Capsule types also influence conjugation efficiency in both donor and recipient cells, a mechanism shaped by capsule volume and conjugative pilus structure. Comparative genomics confirmed that more permissive serotypes in the lab correspond to the strains acquiring more conjugative plasmids in nature. The least capsule-sensitive pili (F-like) are the most frequent in the species' plasmids, and are the only ones associated with both antibiotic resistance and virulence factors, driving the convergence between virulence and antibiotics resistance in the population. These results show how traits of cellular envelopes define slow and fast lanes of infection by mobile genetic elements, with implications for population dynamics and horizontal gene transfer.

Co-evolution of gene transfer agents and their alphaproteobacterial hosts

Gene transfer agents (GTAs) are enigmatic elements that resemble small viruses and are known to be produced during nutritional stress by some bacteria and archaea. The production of GTAs is regulated by quorum sensing, under which a small fraction of the population acts as GTA producers, while the rest becomes GTA recipients. In contrast to canonical viruses, GTAs cannot propagate themselves because they package pieces of the producing cell's genome. In alphaproteobacteria, GTAs are mostly vertically inherited and reside in their hosts' genomes for hundreds of millions of years. While GTAs' ability to transfer genetic material within a population and their long-term preservation suggest an increased fitness of GTA-producing microbes, the associated benefits and type of selection that maintains GTAs are poorly understood. By comparing rates of evolutionary change in GTA genes to the rates in gene families abundantly present across 293 alphaproteobacterial genomes, we detected 59 gene families that likely co-evolve with GTA genes. These gene families are predominantly involved in stress response, DNA repair, and biofilm formation. We hypothesize that biofilm formation enables the physical proximity of GTA-producing cells, limiting GTA-derived benefits only to a group of closely related cells. We further conjecture that the population structure of biofilm-forming sub-populations ensures that the trait of GTA production is maintained despite the inevitable rise of "cheating" genotypes. Because release of GTA particles kills the producing cell, maintenance of GTAs is an exciting example of social evolution in a microbial population.IMPORTANCEGene transfer agents (GTAs) are viruses domesticated by some archaea and bacteria as vehicles for carrying pieces of the host genome. Produced under certain environmental conditions, GTA particles can deliver DNA to neighboring, closely related cells. The function of GTAs remains uncertain. While making GTAs is suicidal for a cell, GTA-encoding genes are widespread in genomes of alphaproteobacteria. Such GTA persistence implies functional benefits but raises questions about how selection maintains this lethal trait. By showing that GTA genes co-evolve with genes involved in stress response, DNA repair, and biofilm formation, we provide support for the hypothesis that GTAs facilitate DNA exchange during the stress conditions and present a model for how GTAs persist in biofilm-forming bacterial populations despite being lethal.

Dynamics of bacterial recombination in the human gut microbiome

Horizontal gene transfer (HGT) is a ubiquitous force in microbial evolution. Previous work has shown that the human gut is a hotspot for gene transfer between species, but the more subtle exchange of variation within species-also known as recombination-remains poorly characterized in this ecosystem. Here, we show that the genetic structure of the human gut microbiome provides an opportunity to measure recent recombination events from sequenced fecal samples, enabling quantitative comparisons across diverse commensal species that inhabit a common environment. By analyzing recent recombination events in the core genomes of 29 human gut bacteria, we observed widespread heterogeneities in the rates and lengths of transferred fragments, which are difficult to explain by existing models of ecological isolation or homology-dependent recombination rates. We also show that natural selection helps facilitate the spread of genetic variants across strain backgrounds, both within individual hosts and across the broader population. These results shed light on the dynamics of in situ recombination, which can strongly constrain the adaptability of gut microbial communities.

Genome-wide transformation reveals extensive exchange across closely related Bacillus species

Bacterial transformation is an important mode of horizontal gene transfer that helps spread genetic material across species boundaries. Yet, the factors that pose barriers to genome-wide cross-species gene transfer are poorly characterized. Here, we develop a replacement accumulation assay to study the effects of genomic distance on transfer dynamics. Using Bacillus subtilis as recipient and various species of the genus Bacillus as donors, we find that the rate of orthologous replacement decreases exponentially with the divergence of their core genomes. We reveal that at least 96% of the B. subtilis core genes are accessible to replacement by alleles from Bacillus spizizenii. For the more distantly related Bacillus atrophaeus, gene replacement events cluster at genomic locations with high sequence identity and preferentially replace ribosomal genes. Orthologous replacement also creates mosaic patterns between donor and recipient genomes, rearranges the genome architecture, and governs gain and loss of accessory genes. We conclude that cross-species gene transfer is dominated by orthologous replacement of core genes which occurs nearly unrestricted between closely related species. At a lower rate, the exchange of accessory genes gives rise to more complex genome dynamics.

Sequence-discrete species for prokaryotes and other microbes: A historical perspective and pending issues

Whether prokaryotes, and other microorganisms, form distinct clusters that can be recognized as species remains an issue of paramount theoretical as well as practical consequence in identifying, regulating, and communicating about these organisms. In the past decade, comparisons of thousands of genomes of isolates and hundreds of metagenomes have shown that prokaryotic diversity may be predominantly organized in such sequence-discrete clusters, albeit organisms of intermediate relatedness between the identified clusters are also frequently found. Accumulating evidence suggests, however, that the latter "intermediate" organisms show enough ecological and/or functional distinctiveness to be considered different species. Notably, the area of discontinuity between clusters often-but not always-appears to be around 85%-95% genome-average nucleotide identity, consistently among different taxa. More recent studies have revealed remarkably similar diversity patterns for viruses and microbial eukaryotes as well. This high consistency across taxa implies a specific mechanistic process that underlies the maintenance of the clusters. The underlying mechanism may be a substantial reduction in the efficiency of homologous recombination, which mediates (successful) horizontal gene transfer, around 95% nucleotide identity. Deviations from the 95% threshold (e.g., species showing lower intraspecies diversity) may be caused by ecological differentiation that imposes barriers to otherwise frequent gene transfer. While this hypothesis that clusters are driven by ecological differentiation coupled to recombination frequency (i.e., higher recombination within vs. between groups) is appealing, the supporting evidence remains anecdotal. The data needed to rigorously test the hypothesis toward advancing the species concept are also outlined.

Phages for treatment of Escherichia coli infections

Diseases due to infections by pathogenic Escherichia coli strains are on the rise and with the growing antimicrobial resistance among bacterial pathogens, including this group. Thus, alternative therapeutic options are actively investigated. Among these alternatives is phage therapy. In the case of E. coli, the combination of the well understood biology of this species and its bacteriophages represents a good guiding example for the establishment of phage therapy principles against this and other pathogenic bacteria. In this chapter, the procedures toward the development of phage therapy against pathogenic E. coli with the use of T-even group of phages are discussed. These steps involve the isolation, purification, characterisation and large-scale production of these phages, with formulation of phage cocktails for in vitro and in vivo studies. The main emphasis is made on phage therapy of enteropathogenic E. coli O157:H, which is one of the prominent human pathogens but persists as a commensal bacterium in many food animals. The implementation of phage therapy against E. coli O157:H within the One Health framework in carrier animals and for treatment of meat, vegetables, fruits and other agricultural produce thus would allow controlling and interrupting the transmission routes of this pathogen to the human food chain and preventing human disease. Examples of successful control and elimination of E. coli O157:H are given, while the problems encountered in phage treatment of this pathogen are also discussed.

Horizontal gene transfer among host-associated microbes

Horizontal gene transfer is an important evolutionary force, facilitating bacterial diversity. It is thought to be pervasive in host-associated microbiomes, where bacterial densities are high and mobile elements are frequent. These genetic exchanges are also key for the rapid dissemination of antibiotic resistance. Here, we review recent studies that have greatly extended our knowledge of the mechanisms underlying horizontal gene transfer, the ecological complexities of a network of interactions involving bacteria and their mobile elements, and the effect of host physiology on the rates of genetic exchanges. Furthermore, we discuss other, fundamental challenges in detecting and quantifying genetic exchanges in vivo, and how studies have contributed to start overcoming these challenges. We highlight the importance of integrating novel computational approaches and theoretical models with experimental methods where multiple strains and transfer elements are studied, both in vivo and in controlled conditions that mimic the intricacies of host-associated environments.

Selfish, promiscuous and sometimes useful: how mobile genetic elements drive horizontal gene transfer in microbial populations

Horizontal gene transfer (HGT) drives microbial adaptation but is often under the control of mobile genetic elements (MGEs) whose interests are not necessarily aligned with those of their hosts. In general, transfer is costly to the donor cell while potentially beneficial to the recipients. The diversity and plasticity of cell–MGEs interactions, and those among MGEs, result in complex evolutionary processes where the source, or even the existence of selection for maintaining a function in the genome, is often unclear. For example, MGE-driven HGT depends on cell envelope structures and defense systems, but many of these are transferred by MGEs themselves. MGEs can spur periods of intense gene transfer by increasing their own rates of horizontal transmission upon communicating, eavesdropping, or sensing the environment and the host physiology. This may result in high-frequency transfer of host genes unrelated to the MGE. Here, we review how MGEs drive HGT and how their transfer mechanisms, selective pressures and genomic traits affect gene flow, and therefore adaptation, in microbial populations. The encoding of many adaptive niche-defining microbial traits in MGEs means that intragenomic conflicts and alliances between cells and their MGEs are key to microbial functional diversification.

This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

L’impact des mutations neutres sur l’évolvabilité et l’évolution des génomes

Les mutations bénéfiques à forts effets sont rares et les mutations délétères sont éliminées par la sélection naturelle. La majorité des mutations qui s’accumulent dans les génomes ont donc des effets sélectifs très faibles, voire nuls ; elles sont alors appelées mutations neutres. Au cours des deux dernières décennies, il a été montré que les mutations, même en l’absence d’effet sur la valeur sélective des organismes, affectent leur évolvabilité, en donnant accès à de nouveaux phénotypes par le biais de mutations apparaissant ultérieurement, et qui n’auraient pas été disponibles autrement. En plus de cet effet, de nombreuses mutations neutres – indépendamment de leurs effets sélectifs – peuvent affecter la mutabilité de séquences d’ADN voisines, et moduler l’efficacité de la recombinaison homologue. De telles mutations ne modifient pas le spectre des phénotypes accessibles, mais plutôt la vitesse à laquelle de nouveaux phénotypes seront produits, un processus qui a des conséquences à long terme mais aussi potentiellement à court terme, en lien avec l’émergence de cancers.

Core genes can have higher recombination rates than accessory genes within global microbial populations

Recombination is essential to microbial evolution, and is involved in the spread of antibiotic resistance, antigenic variation, and adaptation to the host niche. However, assessing the impact of homologous recombination on accessory genes which are only present in a subset of strains of a given species remains challenging due to their complex phylogenetic relationships. Quantifying homologous recombination for accessory genes (which are important for niche-specific adaptations) in comparison to core genes (which are present in all strains and have essential functions) is critical to understanding how selection acts on variation to shape species diversity and genome structures of bacteria. Here, we apply a computationally efficient, non-phylogenetic approach to measure homologous recombination rates in the core and accessory genome using >100,000 whole genome sequences from Streptococcus pneumoniae and several additional species. By analyzing diverse sets of sequence clusters, we show that core genes often have higher recombination rates than accessory genes, and for some bacterial species the associated effect sizes for these differences are pronounced. In a subset of species, we find that gene frequency and homologous recombination rate are positively correlated. For S. pneumoniae and several additional species, we find that while the recombination rate is higher for the core genome, the mutational divergence is lower, indicating that divergence-based homologous recombination barriers could contribute to differences in recombination rates between the core and accessory genome. Homologous recombination may therefore play a key role in increasing the efficiency of selection in the most conserved parts of the genome.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Horizontal Gene Transfer in Archaea-From Mechanisms to Genome Evolution

Archaea remains the least-studied and least-characterized domain of life despite its significance not just to the ecology of our planet but also to the evolution of eukaryotes. It is therefore unsurprising that research into horizontal gene transfer (HGT) in archaea has lagged behind that of bacteria. Indeed, several archaeal lineages may owe their very existence to large-scale HGT events, and thus understanding both the molecular mechanisms and the evolutionary impact of HGT in archaea is highly important. Furthermore, some mechanisms of gene exchange, such as plasmids that transmit themselves via membrane vesicles and the formation of cytoplasmic bridges that allows transfer of both chromosomal and plasmid DNA, may be archaea specific. This review summarizes what we know about HGT in archaea, and the barriers that restrict it, highlighting exciting recent discoveries and pointing out opportunities for future research. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Horizontal gene transfer and adaptive evolution in bacteria

Horizontal gene transfer (HGT) is arguably the most conspicuous feature of bacterial evolution. Evidence for HGT is found in most bacterial genomes. Although HGT can considerably alter bacterial genomes, not all transfer events may be biologically significant and may instead represent the outcome of an incessant evolutionary process that only occasionally has a beneficial purpose. When adaptive transfers occur, HGT and positive selection may result in specific, detectable signatures in genomes, such as gene-specific sweeps or increased transfer rates for genes that are ecologically relevant. In this Review, we first discuss the various mechanisms whereby HGT occurs, how the genetic signatures shape patterns of genomic variation and the distinct bioinformatic algorithms developed to detect these patterns. We then discuss the evolutionary theory behind HGT and positive selection in bacteria, and discuss the approaches developed over the past decade to detect transferred DNA that may be involved in adaptation to new environments.

Speciation of Genes and Genomes: Conservation of DNA Polymorphism by Barriers to Recombination Raised by Mismatch Repair System

Some basic aspects of human and animal biology and evolution involve the establishment of biological uniqueness of species and individuals within their huge variety. The discrimination among closely related species occurs in their offspring at the level of chromosomal DNA sequence homology, which is required for fertility as the hallmark of species. Biological identification of individuals, i.e., of their biological “self”, occurs at the level of protein sequences presented by the MHC/HLA complex as part of the immune system that discriminates non-self from self. Here, a mechanistic molecular model is presented that can explain how DNA sequence divergence and the activity of key mismatch repair proteins, MutS and MutL, lead to 1) genetic separation of closely related species (sympatric speciation) (Fitch and Ayala, Proceedings of the National Academy of Sciences, 1994, 91, 6717–6720), 2) the stability of genomes riddled by diverged repeated sequences, and 3) conservation of highly polymorphic DNA sequence blocks that constitute the immunological self. All three phenomena involve suppression of recombination between diverged homologies, resulting in prevention of gene sharing between closely related genomes (evolution of new species) as well as sequence sharing between closely related genes within a genome (e.g., evolution of immunoglobulin, MHC, and other gene families bearing conserved polymorphisms).

Adaptive evolution of hybrid bacteria by horizontal gene transfer

Horizontal gene transfer (HGT) is an important factor in bacterial evolution that can act across species boundaries. Yet, we know little about rate and genomic targets of cross-lineage gene transfer and about its effects on the recipient organism's physiology and fitness. Here, we address these questions in a parallel evolution experiment with two Bacillus subtilis lineages of 7% sequence divergence. We observe rapid evolution of hybrid organisms: gene transfer swaps ∼12% of the core genome in just 200 generations, and 60% of core genes are replaced in at least one population. By genomics, transcriptomics, fitness assays, and statistical modeling, we show that transfer generates adaptive evolution and functional alterations in hybrids. Specifically, our experiments reveal a strong, repeatable fitness increase of evolved populations in the stationary growth phase. By genomic analysis of the transfer statistics across replicate populations, we infer that selection on HGT has a broad genetic basis: 40% of the observed transfers are adaptive. At the level of functional gene networks, we find signatures of negative, positive, and epistatic selection, consistent with hybrid incompatibilities and adaptive evolution of network functions. Our results suggest that gene transfer navigates a complex cross-lineage fitness landscape, bridging epistatic barriers along multiple high-fitness paths.

Genome-wide dissection reveals diverse pathogenic roles of bacterial Tc toxins

Tc toxins were originally identified in entomopathogenic bacteria, which are important as biological pest control agents. Tc toxins are heteromeric exotoxins composed of three subunit types, TcA, TcB, and TcC. The C-terminal portion of the TcC protein encodes the actual toxic domain, which is translocated into host cells by an injectosome nanomachine comprising the other subunits. Currently the pathogenic roles and distribution of Tc toxins among different bacterial genera remain unclear. Here we have performed a comprehensive genome-wide analysis, and established a database that includes 1,608 identified Tc loci containing 2,528 TcC proteins in 1,421 Gram-negative and positive bacterial genomes. Our findings indicate that TcCs conform to the architecture of typical polymorphic toxins, with C-terminal hypervariable regions (HVR) encoding more than 100 different classes of putative toxic domains, most of which have not been previously recognized. Based on further analysis of Tc loci in the genomes of all Salmonella and Yersinia strains in EnteroBase, a “two-level” evolutionary dynamics scenario is proposed for TcC homologues. This scenario implies that the conserved TcC RHS core domain plays a critical role in the taxonomical specific distribution of TcC HVRs. This study provides an extensive resource for the future development of Tc toxins as valuable agrochemical tools. It furthermore implies that Tc proteins, which are encoded by a wide range of pathogens, represent an important versatile toxin superfamily with diverse pathogenic mechanisms.

Entomopathogenic bacteria deploy a range of toxins to combat their insect hosts. The Tc toxins were first identified in Photorhabdus as having potent oral toxicity to insects, with a mode of action distinct from the well-studied Bacillus thuringiensis Cry toxins. As such the Tc toxins have been considered as potential candidates for novel crop protection strategies. This could mitigate against the potential risks of pest insects developing resistance to the traditionally used Cry toxin-based systems. To date, the generality of diverse Tc toxins and their related pathogenic roles has remained mainly obscure. Our analysis has showed Tc toxins are widely distributed among Gram-negative and positive bacterial genomes. A database was constructed including thousands of Tc loci with hundreds of different putative TcC toxic domains, any one of which might represent candidates for the development of future pest control systems. Moreover, the findings of this study are of wider significance because Tc toxin homologues have been shown to be encoded by a range of human pathogens. These include Salmonella and Yersinia, suggesting their potential roles in human infectious diseases. Together, this study describes the characteristics and distribution of Tc toxins among diverse bacterial genera, and provides a new insight into their roles in different pathogenesis mechanisms. This study also describes findings of potential importance to their development as tools for biotechnological applications.

Whole genome phylogenies reflect the distributions of recombination rates for many bacterial species

Although recombination is accepted to be common in bacteria, for many species robust phylogenies with well-resolved branches can be reconstructed from whole genome alignments of strains, and these are generally interpreted to reflect clonal relationships. Using new methods based on the statistics of single-nucleotide polymorphism (SNP) splits, we show that this interpretation is incorrect. For many species, each locus has recombined many times along its line of descent, and instead of many loci supporting a common phylogeny, the phylogeny changes many thousands of times along the genome alignment. Analysis of the patterns of allele sharing among strains shows that bacterial populations cannot be approximated as either clonal or freely recombining but are structured such that recombination rates between lineages vary over several orders of magnitude, with a unique pattern of rates for each lineage. Thus, rather than reflecting clonal ancestry, whole genome phylogenies reflect distributions of recombination rates.

CoreSimul: a forward-in-time simulator of genome evolution for prokaryotes modeling homologous recombination

Background

Prokaryotes are asexual, but these organisms frequently engage in homologous recombination, a process that differs from meiotic recombination in sexual organisms. Most tools developed to simulate genome evolution either assume sexual reproduction or the complete absence of DNA flux in the population. As a result, very few simulators are adapted to model prokaryotic genome evolution while accounting for recombination. Moreover, many simulators are based on the coalescent, which assumes a neutral model of genomic evolution, and those are best suited for organisms evolving under weak selective pressures, such as animals and plants. In contrast, prokaryotes are thought to be evolving under much stronger selective pressures, suggesting that forward-in-time simulators are better suited for these organisms.

Results

Here, I present CoreSimul, a forward-in-time simulator of core genome evolution for prokaryotes modeling homologous recombination. Simulations are guided by a phylogenetic tree and incorporate different substitution models, including models of codon selection.

Conclusions

CoreSimul is a flexible forward-in-time simulator that constitutes a significant addition to the limited list of available simulators applicable to prokaryote genome evolution.

The Impact of Neutral Mutations on Genome Evolvability

Beneficial mutations are rare and deleterious mutations are purged by natural selection. As a result, the vast majority of mutations that accumulate in genomes belong to the class of neutral mutations. Over the last two decades, neutral mutations, despite their null effect on fitness, have been shown to affect evolvability by providing access to new phenotypes through subsequent mutations that would not have been available otherwise. Here we propose that in addition, many mutations - independent of their selective effects - can affect the mutability of neighboring DNA sequences and modulate the efficacy of homologous recombination. Such mutations do not change the spectrum of accessible phenotypes, but rather the rate at which new phenotypes will be produced. Therefore, neutral mutations that accumulate in genomes have an important long-term impact on the evolutionary fate of genomes.

Consistent Metagenome-Derived Metrics Verify and Delineate Bacterial Species Boundaries

There is controversy about whether bacterial diversity is clustered into distinct species groups or exists as a continuum. To address this issue, we analyzed bacterial genome databases and reports from several previous large-scale environment studies and identified clear discrete groups of species-level bacterial diversity in all cases. Genetic analysis further revealed that quasi-sexual reproduction via horizontal gene transfer is likely a key evolutionary force that maintains bacterial species integrity. We next benchmarked over 100 metrics to distinguish these bacterial species from each other and identified several genes encoding ribosomal proteins with high species discrimination power. Overall, the results from this study provide best practices for bacterial species delineation based on genome content and insight into the nature of bacterial species population genetics.

Longstanding questions relate to the existence of naturally distinct bacterial species and genetic approaches to distinguish them. Bacterial genomes in public databases form distinct groups, but these databases are subject to isolation and deposition biases. To avoid these biases, we compared 5,203 bacterial genomes from 1,457 environmental metagenomic samples to test for distinct clouds of diversity and evaluated metrics that could be used to define the species boundary. Bacterial genomes from the human gut, soil, and the ocean all exhibited gaps in whole-genome average nucleotide identities (ANI) near the previously suggested species threshold of 95% ANI. While genome-wide ratios of nonsynonymous and synonymous nucleotide differences (dN/dS) decrease until ANI values approach ∼98%, two methods for estimating homologous recombination approached zero at ∼95% ANI, supporting breakdown of recombination due to sequence divergence as a species-forming force. We evaluated 107 genome-based metrics for their ability to distinguish species when full genomes are not recovered. Full-length 16S rRNA genes were least useful, in part because they were underrecovered from metagenomes. However, many ribosomal proteins displayed both high metagenomic recoverability and species discrimination power. Taken together, our results verify the existence of sequence-discrete microbial species in metagenome-derived genomes and highlight the usefulness of ribosomal genes for gene-level species discrimination.

IMPORTANCE There is controversy about whether bacterial diversity is clustered into distinct species groups or exists as a continuum. To address this issue, we analyzed bacterial genome databases and reports from several previous large-scale environment studies and identified clear discrete groups of species-level bacterial diversity in all cases. Genetic analysis further revealed that quasi-sexual reproduction via horizontal gene transfer is likely a key evolutionary force that maintains bacterial species integrity. We next benchmarked over 100 metrics to distinguish these bacterial species from each other and identified several genes encoding ribosomal proteins with high species discrimination power. Overall, the results from this study provide best practices for bacterial species delineation based on genome content and insight into the nature of bacterial species population genetics.

Systematics: The Cohesive Nature of Bacterial Species Taxa

A survey of bacterial genomes suggests that the diversity within recognized species is constrained by a force of cohesion. However, recognized bacterial species do not adhere to another species-like property-that of being the newest lineages that can coexist indefinitely.

Gene gain and loss push prokaryotes beyond the homologous recombination barrier and accelerate genome sequence divergence

Bacterial and archaeal evolution involve extensive gene gain and loss. Thus, phylogenetic trees of prokaryotes can be constructed both by traditional sequence-based methods (gene trees) and by comparison of gene compositions (genome trees). Comparing the branch lengths in gene and genome trees with identical topologies for 34 clusters of closely related bacterial and archaeal genomes, we show here that terminal branches of gene trees are systematically compressed compared to those of genome trees. Thus, sequence evolution is delayed compared to genome evolution by gene gain and loss. The extent of this delay differs widely among bacteria and archaea. Mathematical modeling shows that the divergence delay can result from sequence homogenization by homologous recombination. The model explains how homologous recombination maintains the cohesiveness of the core genome of a species while allowing extensive gene gain and loss within the accessory genome. Once evolving genomes become isolated by barriers impeding homologous recombination, gene and genome evolution processes settle into parallel trajectories, and genomes diverge, resulting in speciation.

A significant proportion of the molecular evolution of bacteria and archaea occurs through gene gain and loss. Here Iranzo et al. develop a mathematical model that explains observed differential patterns of sequence evolution vs. gene content evolution as a consequence of homologous recombination.

Mismatch recognition and subsequent processing have distinct effects on mitotic recombination intermediates and outcomes in yeast

The post-replicative mismatch repair (MMR) system has anti-recombination activity that limits interactions between diverged sequences by recognizing mismatches in strand-exchange intermediates. In contrast to their equivalent roles during replication-error repair, mismatch recognition is more important for anti-recombination than subsequent mismatch processing. To obtain insight into this difference, ectopic substrates with 2% sequence divergence were used to examine mitotic recombination outcome (crossover or noncrossover; CO and NCO, respectively) and to infer molecular intermediates formed during double-strand break repair in Saccharomyces cerevisiae. Experiments were performed in an MMR-proficient strain, a strain with compromised mismatch-recognition activity (msh6Δ) and a strain that retained mismatch-recognition activity but was unable to process mismatches (mlh1Δ). While the loss of either mismatch binding or processing elevated the NCO frequency to a similar extent, CO events increased only when mismatch binding was compromised. The molecular features of NCOs, however, were altered in fundamentally different ways depending on whether mismatch binding or processing was eliminated. These data suggest a model in which mismatch recognition reverses strand-exchange intermediates prior to the initiation of end extension, while subsequent mismatch processing that is linked to end extension specifically destroys NCO intermediates that contain conflicting strand-discrimination signals for mismatch removal.

Adapting the engine to the fuel: mutator populations can reduce the mutational load by reorganizing their genome structure

Background

Mutators are common in bacterial populations, both in natural isolates and in the lab. The fate of these lineages, which mutation rate is increased up to 100 ×, has long been studied using population genetics models, showing that they can spread in a population following an environmental change. However in stable conditions, they suffer from the increased mutational load, hence being overcome by non-mutators. However, these results don’t take into account the fact that an elevated mutation rate can impact the genetic structure, hence changing the sensitivity of the population to mutations. Here we used Aevol, an in silico experimental evolution platform in which genomic structures are free to evolve, in order to study the fate of mutator populations evolving for a long time in constant conditions.

Results

Starting from wild-types that were pre-evolved for 300,000 generations, we let 100 mutator populations (point mutation rate ×100) evolve for 100,000 further generations in constant conditions. As expected all populations initially undergo a fitness loss. However, after that the mutator populations started to recover. Most populations ultimately recovered their ancestors fitness, and a significant fraction became even fitter than the non-mutator control clones that evolved in parallel. By analyzing the genomes of the mutators, we show that the fitness recovery is due to two mechanisms: i. an increase in robustness through compaction of the coding part of the mutator genomes, ii. an increase of the selection coefficient that decreases the mean-fitness of the population. Strikingly the latter is due to the accumulation of non-coding sequences in the mutators genomes.

Conclusion

Our results show that the mutational burden that is classically thought to be associated with mutator phenotype is escapable. On the long run mutators adapted their genomes and reshaped the distribution of mutation effects. Therewith the lineage is able to recover fitness even though the population still suffers the elevated mutation rate. Overall these results change our view of mutator dynamics: by being able to reduce the deleterious effect of the elevated mutation rate, mutator populations may be able to last for a very long time; A situation commonly observed in nature.

The speciation and hybridization history of the genus Salmonella

Bacteria and archaea make up most of natural diversity, but the mechanisms that underlie the origin and maintenance of prokaryotic species are poorly understood. We investigated the speciation history of the genus Salmonella, an ecologically diverse bacterial lineage, within which S. enterica subsp. enterica is responsible for important human food-borne infections. We performed a survey of diversity across a large reference collection using multilocus sequence typing, followed by genome sequencing of distinct lineages. We identified 11 distinct phylogroups, 3 of which were previously undescribed. Strains assigned to S. enterica subsp. salamae are polyphyletic, with two distinct lineages that we designate Salamae A and B. Strains of the subspecies houtenae are subdivided into two groups, Houtenae A and B, and are both related to Selander’s group VII. A phylogroup we designate VIII was previously unknown. A simple binary fission model of speciation cannot explain observed patterns of sequence diversity. In the recent past, there have been large-scale hybridization events involving an unsampled ancestral lineage and three distantly related lineages of the genus that have given rise to Houtenae A, Houtenae B and VII. We found no evidence for ongoing hybridization in the other eight lineages, but detected subtler signals of ancient recombination events. We are unable to fully resolve the speciation history of the genus, which might have involved additional speciation-by-hybridization or multi-way speciation events. Our results imply that traditional models of speciation by binary fission and divergence are not sufficient to account for Salmonella evolution.

Bacillus subtilis MutS Modulates RecA-Mediated DNA Strand Exchange Between Divergent DNA Sequences

The efficiency of horizontal gene transfer, which contributes to acquisition and spread of antibiotic resistance and pathogenicity traits, depends on nucleotide sequence and different mismatch-repair (MMR) proteins participate in this process. To study how MutL and MutS MMR proteins regulate recombination across species boundaries, we have studied natural chromosomal transformation with DNA up to ∼23% sequence divergence. We show that Bacillus subtilis natural chromosomal transformation decreased logarithmically with increased sequence divergence up to 15% in wild type (wt) cells or in cells lacking MutS2 or mismatch repair proteins (MutL, MutS or both). Beyond 15% sequence divergence, the chromosomal transformation efficiency is ∼100-fold higher in ΔmutS and ΔmutSL than in ΔmutS2 or wt cells. In the first phase of the biphasic curve (up to 15% sequence divergence), RecA-catalyzed DNA strand exchange contributes to the delineation of species, and in the second phase, homology-facilitated illegitimate recombination might aid in the restoration of inactivated genes. To understand how MutS modulates the integration process, we monitored DNA strand exchange reactions using a circular single-stranded DNA and a linear double-stranded DNA substrate with an internal 77-bp region with ∼16% or ∼54% sequence divergence in an otherwise homologous substrate. The former substrate delayed, whereas the latter halted RecA-mediated strand exchange. Interestingly, MutS addition overcame the heterologous barrier. We propose that MutS assists DNA strand exchange by facilitating RecA disassembly, and indirectly re-engagement with the homologous 5′-end of the linear duplex. Our data supports the idea that MutS modulates bidirectional RecA-mediated integration of divergent sequences and this is important for speciation.

Cross-Contamination Explains "Inter and Intraspecific Horizontal Genetic Transfers" between Asexual Bdelloid Rotifers

A few metazoan lineages are thought to have persisted for millions of years without sexual reproduction. If so, they would offer important clues to the evolutionary paradox of sex itself [1, 2]. Most "ancient asexuals" are subject to ongoing doubt because extant populations continue to invest in males [3-9]. However, males are famously unknown in bdelloid rotifers, a class of microscopic invertebrates comprising hundreds of species [10-12]. Bdelloid genomes have acquired an unusually high proportion of genes from non-metazoans via horizontal transfer [13-17]. This well-substantiated finding has invited speculation [13] that homologous horizontal transfer between bdelloid individuals also may occur, perhaps even "replacing" sex [14]. In 2016, Current Biology published an article claiming to supply evidence for this idea. Debortoli et al. [18] sampled rotifers from natural populations and sequenced one mitochondrial and four nuclear loci. Species assignments were incongruent among loci for several samples, which was interpreted as evidence of "interspecific horizontal genetic transfers." Here, we use sequencing chromatograms supplied by the authors to demonstrate that samples treated as individuals actually contained two or more highly divergent mitochondrial and ribosomal sequences, revealing cross-contamination with DNA from multiple animals of different species. Other chromatograms indicate contamination with DNA from conspecific animals, explaining genetic and genomic evidence for "intraspecific horizontal exchanges" reported in the same study. Given the clear evidence of contamination, the data and findings of Debortoli et al. [18] provide no reliable support for their conclusions that DNA is transferred horizontally between or within bdelloid species.

Role of the DNA Mismatch Repair Gene MutS4 in Driving the Evolution of Mycobacterium yongonense Type I via Homologous Recombination

We recently showed that Mycobacterium yongonense could be divided into two genotypes: Type I, in which the rpoB gene has been transferred from Mycobacterium parascrofulaceum, and Type II, in which the rpoB gene has not been transferred. Comparative genome analysis of three M. yongonense Type I, two M. yongonense Type II and M. parascrofulaceum type strains were performed in this study to gain insight into gene transfer from M. parascrofulaceum into M. yongonense Type I strains. We found two genome regions transferred from M. parascrofulaceum: one contained 3 consecutive genes, including the rpoBC operon, and the other contained 57 consecutive genes that had been transferred into M. yongonense Type I genomes via homologous recombination. Further comparison between the M. yongonense Type I and II genomes revealed that Type I, but not Type II has a distinct DNA mismatch repair gene (MutS4 subfamily) that was possibly transferred via non-homologous recombination from other actinomycetes. We hypothesized that it could facilitate homologous recombination from the M. parascrofulaceum to the M. yongonense Type I genomes. We therefore generated recombinant Mycobacterium smegmatis containing a MutS4 operon of M. yongonense. We found that the M. tuberculosis rpoB fragment with a rifampin resistance-conferring mutation was more frequently inserted into recombinant M. smegmatis than the wild type, suggesting that MutS4 is a driving force in the gene transfer from M. parascrofulaceum to M. yongonense Type I strains via homologous recombination. In conclusion, our data indicated that MutS4 in M. yongonense Type I genomes may drive gene transfer from M. parascrofulaceum via homologous recombination, resulting in division of M. yongonense into two genotypes, Type I and II.

Steady at the wheel: conservative sex and the benefits of bacterial transformation

Many bacteria are highly sexual, but the reasons for their promiscuity remain obscure. Did bacterial sex evolve to maximize diversity and facilitate adaptation in a changing world, or does it instead help to retain the bacterial functions that work right now? In other words, is bacterial sex innovative or conservative? Our aim in this review is to integrate experimental, bioinformatic and theoretical studies to critically evaluate these alternatives, with a main focus on natural genetic transformation, the bacterial equivalent of eukaryotic sexual reproduction. First, we provide a general overview of several hypotheses that have been put forward to explain the evolution of transformation. Next, we synthesize a large body of evidence highlighting the numerous passive and active barriers to transformation that have evolved to protect bacteria from foreign DNA, thereby increasing the likelihood that transformation takes place among clonemates. Our critical review of the existing literature provides support for the view that bacterial transformation is maintained as a means of genomic conservation that provides direct benefits to both individual bacterial cells and to transformable bacterial populations. We examine the generality of this view across bacteria and contrast this explanation with the different evolutionary roles proposed to maintain sex in eukaryotes. This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Genomic divergence and cohesion in a species of pelagic freshwater bacteria

Background

In many prokaryotic genera a clustered phylogeny is observed, akin to the occurrence of species in sexually reproducing organisms. For some taxa, homologous recombination has been invoked as the underlying mechanism providing genomic cohesion among conspecific individuals. Whether this mechanism is applicable to prokaryotes in freshwaters with low habitat connectivity – i.e. elevated geographic barriers to gene flow – is unclear. To investigate further we studied genomic trends within the globally abundant PnecC cluster (genus Polynucleobacter, Betaproteobacteria) and analyzed homologous recombination within the affiliated species P. asymbioticus.

Results

Comparisons among 20 PnecC genomes revealed a clearly discontinuous distribution of nucleotide sequence similarities. Among the nine conspecific individuals (P. asymbioticus) all average nucleotide identity (ANI) values were greater than 97%, whereas all other comparisons exhibited ANI values lower than 85%. The reconstruction of recombination and mutation events for the P. asymbioticus core genomes yielded an r/m ratio of 7.4, which is clearly above estimated thresholds for recombination to act as a cohesive force. Hotspots of recombination were found to be located in the flanking regions of genomic islands. Even between geographically separated habitats a high flux of recombination was evident. While a biogeographic population structure was suggested from MLST data targeting rather conserved loci, such a structure was barely visible when whole genome data was considered. However, both MLST and whole genome data showed evidence of differentiation between two lineages of P. asymbioticus. The ratios of non-synonymous to synonymous substitution rates as well as growth rates in transplantation experiments suggested that this divergence was not selectively neutral.

Conclusions

The high extent of homologous recombination among P. asymbioticus bacteria can act as a cohesive force that effectively counteracts genetic divergence. At least on a regional scale, homologous recombination can act across geographically separated ecosystems and therefore plays an important role in the evolution and consistency of bacterial freshwater species. A species model akin to the biological species concept may be applicable for P. asymbioticus. Nonetheless, two genetically distinct lineages have emerged and further research may clarify if their divergence has been initiated by reinforced geographical barriers or has been evolving in sympatry.

Electronic supplementary material

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

Recombination-Driven Genome Evolution and Stability of Bacterial Species

While bacteria divide clonally, horizontal gene transfer followed by homologous recombination is now recognized as an important contributor to their evolution. However, the details of how the competition between clonality and recombination shapes genome diversity remains poorly understood. Using a computational model, we find two principal regimes in bacterial evolution and identify two composite parameters that dictate the evolutionary fate of bacterial species. In the divergent regime, characterized by either a low recombination frequency or strict barriers to recombination, cohesion due to recombination is not sufficient to overcome the mutational drift. As a consequence, the divergence between pairs of genomes in the population steadily increases in the course of their evolution. The species lacks genetic coherence with sexually isolated clonal subpopulations continuously formed and dissolved. In contrast, in the metastable regime, characterized by a high recombination frequency combined with low barriers to recombination, genomes continuously recombine with the rest of the population. The population remains genetically cohesive and temporally stable. Notably, the transition between these two regimes can be affected by relatively small changes in evolutionary parameters. Using the Multi Locus Sequence Typing (MLST) data, we classify a number of bacterial species to be either the divergent or the metastable type. Generalizations of our framework to include selection, ecologically structured populations, and horizontal gene transfer of nonhomologous regions are discussed as well.

Speciation trajectories in recombining bacterial species

It is generally agreed that bacterial diversity can be classified into genetically and ecologically cohesive units, but what produces such variation is a topic of intensive research. Recombination may maintain coherent species of frequently recombining bacteria, but the emergence of distinct clusters within a recombining species, and the impact of habitat structure in this process are not well described, limiting our understanding of how new species are created. Here we present a model of bacterial evolution in overlapping habitat space. We show that the amount of habitat overlap determines the outcome for a pair of clusters, which may range from fast clonal divergence with little interaction between the clusters to a stationary population structure, where different clusters maintain an equilibrium distance between each other for an indefinite time. We fit our model to two data sets. In Streptococcus pneumoniae, we find a genomically and ecologically distinct subset, held at a relatively constant genetic distance from the majority of the population through frequent recombination with it, while in Campylobacter jejuni, we find a minority population we predict will continue to diverge at a higher rate. This approach may predict and define speciation trajectories in multiple bacterial species.

Evolution of commensal bacteria in the intestinal tract of mice

Hundreds of different bacterial species inhabit our intestines and contribute to our health status, with significant loss of species diversity typically observed in disease conditions. Within each microbial species a great deal of diversity is hidden and such intra-specific variation is also key to the proper homeostasis between the host and its microbial inhabitants. Indeed, it is at this level that new mechanisms of antibiotic resistance emerge and pathogenic characteristics evolve. Yet, our knowledge on intra-species variation in the gut is still limited and an understanding of the evolutionary mechanisms acting on it is extremely reduced. Here we review recent work that has begun to reveal that adaptation of commensal bacteria to the mammalian intestine may be fast and highly repeatable, and that the time scales of evolutionary and ecological change can be very similar in these ecosystems.

Recombination Processes and Nonlinear Markov Chains

Bacteria are known to exchange genetic information by horizontal gene transfer. Since the frequency of homologous recombination depends on the similarity between the recombining segments, several studies examined whether this could lead to the emergence of subspecies. Most of them simulated fixed-size Wright-Fisher populations, in which the genetic drift should be taken into account. Here, we use nonlinear Markov processes to describe a bacterial population evolving under mutation and recombination. We consider a population structure as a probability measure on the space of genomes. This approach implies the infinite population size limit, and thus, the genetic drift is not assumed. We prove that under these conditions, the emergence of subspecies is impossible.

Regulation of genetic flux between bacteria by restriction-modification systems

Restriction-modification (R-M) systems are often regarded as bacteria's innate immune systems, protecting cells from infection by mobile genetic elements (MGEs). Their diversification has been recently associated with the emergence of particularly virulent lineages. However, we have previously found more R-M systems in genomes carrying more MGEs. Furthermore, it has been suggested that R-M systems might favor genetic transfer by producing recombinogenic double-stranded DNA ends. To test whether R-M systems favor or disfavor genetic exchanges, we analyzed their frequency with respect to the inferred events of homologous recombination and horizontal gene transfer within 79 bacterial species. Genetic exchanges were more frequent in bacteria with larger genomes and in those encoding more R-M systems. We created a recognition target motif predictor for Type II R-M systems that identifies genomes encoding systems with similar restriction sites. We found more genetic exchanges between these genomes, independently of their evolutionary distance. Our results reconcile previous studies by showing that R-M systems are more abundant in promiscuous species, wherein they establish preferential paths of genetic exchange within and between lineages with cognate R-M systems. Because the repertoire and/or specificity of R-M systems in bacterial lineages vary quickly, the preferential fluxes of genetic transfer within species are expected to constantly change, producing time-dependent networks of gene transfer.

How clonal are bacteria over time?

Bacteria and archaea reproduce clonally (vertical descent), but exchange genes by recombination (horizontal transfer). Recombination allows adaptive mutations or genes to spread rapidly within (or even between) species, and reduces the burden of deleterious mutations. Clonality-defined here as the balance between vertical and horizontal inheritance-is therefore a key microbial trait, determining how quickly a population can adapt and the size of its gene pool. Here, I discuss whether clonality varies over time and if it can be considered a stable trait of a given population. I show that, in some cases, clonality is clearly not static. For example, non-clonal (highly recombining) populations can give rise to clonal expansions, often of pathogens. However, an analysis of time-course metagenomic data from a lake suggests that a bacterial population's past clonality (as measured by its genetic diversity) is a good predictor of its future clonality. Clonality therefore appears to be relatively-but not completely-stable over evolutionary time.

Recombination produces coherent bacterial species clusters in both core and accessory genomes

BACKGROUND

Population samples show bacterial genomes can be divided into a core of ubiquitous genes and accessory genes that are present in a fraction of isolates. The ecological significance of this variation in gene content remains unclear. However, microbiologists agree that a bacterial species should be 'genomically coherent', even though there is no consensus on how this should be determined.

RESULTS

We use a parsimonious model combining diversification in both the core and accessory genome, including mutation, homologous recombination (HR) and horizontal gene transfer (HGT) introducing new loci, to produce a population of interacting clusters of strains with varying genome content. New loci introduced by HGT may then be transferred on by HR. The model fits well to a systematic population sample of 616 pneumococcal genomes, capturing the major features of the population structure with parameter values that agree well with empirical estimates.

CONCLUSIONS

The model does not include explicit selection on individual genes, suggesting that crude comparisons of gene content may be a poor predictor of ecological function. We identify a clearly divergent subpopulation of pneumococci that are inconsistent with the model and may be considered genomically incoherent with the rest of the population. These strains have a distinct disease tropism and may be rationally defined as a separate species. We also find deviations from the model that may be explained by recent population bottlenecks or spatial structure.

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is ‘open’, with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.

Microbial Speciation

What are species? How do they arise? These questions are not easy to answer and have been particularly controversial in microbiology. Yet, for those microbiologists studying environmental questions or dealing with clinical issues, the ability to name and recognize species, widely considered the fundamental units of ecology, can be practically useful. On a more fundamental level, the speciation problem, the focus here, is more mechanistic and conceptual. What is the origin of microbial species, and what evolutionary and ecological mechanisms keep them separate once they begin to diverge? To what extent are these mechanisms universal across diverse types of microbes, and more broadly across the entire the tree of life? Here, we propose that microbial speciation must be viewed in light of gene flow, which defines units of genetic similarity, and of natural selection, which defines units of phenotype and ecological function. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and population genomics studies. These studies suggest a continuous "speciation spectrum," which microbial populations traverse in different ways depending on their balance of gene flow and natural selection.

Recombinant transfer in the basic genome of Escherichia coli

An approximation to the ∼4-Mbp basic genome shared by 32 strains of Escherichia coli representing six evolutionary groups has been derived and analyzed computationally. A multiple alignment of the 32 complete genome sequences was filtered to remove mobile elements and identify the most reliable ∼90% of the aligned length of each of the resulting 496 basic-genome pairs. Patterns of single base-pair mutations (SNPs) in aligned pairs distinguish clonally inherited regions from regions where either genome has acquired DNA fragments from diverged genomes by homologous recombination since their last common ancestor. Such recombinant transfer is pervasive across the basic genome, mostly between genomes in the same evolutionary group, and generates many unique mosaic patterns. The six least-diverged genome pairs have one or two recombinant transfers of length ∼40-115 kbp (and few if any other transfers), each containing one or more gene clusters known to confer strong selective advantage in some environments. Moderately diverged genome pairs (0.4-1% SNPs) show mosaic patterns of interspersed clonal and recombinant regions of varying lengths throughout the basic genome, whereas more highly diverged pairs within an evolutionary group or pairs between evolutionary groups having >1.3% SNPs have few clonal matches longer than a few kilobase pairs. Many recombinant transfers appear to incorporate fragments of the entering DNA produced by restriction systems of the recipient cell. A simple computational model can closely fit the data. Most recombinant transfers seem likely to be due to generalized transduction by coevolving populations of phages, which could efficiently distribute variability throughout bacterial genomes.

Horizontal gene transfer, dispersal and haloarchaeal speciation

The Halobacteria are a well-studied archaeal class and numerous investigations are showing how their diversity is distributed amongst genomes and geographic locations. Evidence indicates that recombination between species continuously facilitates the arrival of new genes, and within species, it is frequent enough to spread acquired genes amongst all individuals in the population. To create permanent independent diversity and generate new species, barriers to recombination are probably required. The data support an interpretation that rates of evolution (e.g., horizontal gene transfer and mutation) are faster at creating geographically localized variation than dispersal and invasion are at homogenizing genetic differences between locations. Therefore, we suggest that recurrent episodes of dispersal followed by variable periods of endemism break the homogenizing forces of intrapopulation recombination and that this process might be the principal stimulus leading to divergence and speciation in Halobacteria.

Systematic identification and quantification of phase variation in commensal and pathogenic Escherichia coli

Bacteria have been shown to generate constant genetic variation in a process termed phase variation. We present a tool based on whole genome sequencing that allows detection and quantification of coexisting genotypes mediated by genomic inversions in bacterial cultures. We tested our method on widely used strains of Escherichia coli, and detected stable and reproducible phase variation in several invertible loci. These are shown here to be responsible for maintaining constant variation in populations grown from a single colony. Applying this tool on other bacterial strains can shed light on how pathogens adjust to hostile environments by diversifying their genomes.