Ecological speciation in bacteria: reverse ecology approaches reveal the adaptive part of bacterial cladogenesis

Combined reference-free and multi-reference based GWAS uncover cryptic variation underlying rapid adaptation in a fungal plant pathogen

Microbial pathogens often harbor substantial functional diversity driven by structural genetic variation. Rapid adaptation from such standing variation threatens global food security and human health. Genome-wide association studies (GWAS) provide a powerful approach to identify genetic variants underlying recent pathogen adaptation. However, the reliance on single reference genomes and single nucleotide polymorphisms (SNPs) obscures the true extent of adaptive genetic variation. Here, we show quantitatively how a combination of multiple reference genomes and reference-free approaches captures substantially more relevant genetic variation compared to single reference mapping. We performed reference-genome based association mapping across 19 reference-quality genomes covering the diversity of the species. We contrasted the results with a reference-free (i.e., k-mer) approach using raw whole-genome sequencing data in a panel of 145 strains collected across the global distribution range of the fungal wheat pathogen Zymoseptoria tritici. We mapped the genetic architecture of 49 life history traits including virulence, reproduction and growth in multiple stressful environments. The inclusion of additional reference genome SNP datasets provides a nearly linear increase in additional loci mapped through GWAS. Variants detected through the k-mer approach explained a higher proportion of phenotypic variation than a reference genome-based approach and revealed functionally confirmed loci that classic GWAS approaches failed to map. The power of GWAS in microbial pathogens can be significantly enhanced by comprehensively capturing structural genetic variation. Our approach is generalizable to a large number of species and will uncover novel mechanisms driving rapid adaptation of pathogens.

Leveraging phylogenetic signal to unravel microbiome function and assembly rules

Clarifying the general rules behind microbial community assembly will foster the development of microbiome-based technological solutions. Here, we study microbial community assembly through a computational analysis of phylogenetic core groups (PCGs): discrete portions of the bacterial phylogeny with high prevalence in the ecosystem under study. We first show that the existence of PCGs was a predominant feature of the varied set of microbial ecosystems studied. Then, we re-analyzed an in vitro experimental dataset using a PCG-based approach, drawing only from its community composition data and from publicly available genomic databases. Using mainly genome scale metabolic models and population dynamics modeling, we obtained ecological insights on metabolic niche structure and population dynamics comparable to those gained after canonical experimentation. Thus, leveraging phylogenetic signal to help unravel microbiome function and assembly rules offers a potential avenue to gain further insight on Earth's microbial ecosystems.

Phylogroup-specific variation shapes the clustering of antimicrobial resistance genes and defence systems across regions of genome plasticity in Pseudomonas aeruginosa

BACKGROUND

Pseudomonas aeruginosa is an opportunistic pathogen consisting of three phylogroups (hereafter named A, B, and C). Here, we assessed phylogroup-specific evolutionary dynamics across available and also new P. aeruginosa genomes.

METHODS

In this genomic analysis, we first generated new genome assemblies for 18 strains of the major P. aeruginosa clone type (mPact) panel, comprising a phylogenetically diverse collection of clinical and environmental isolates for this species. Thereafter, we combined these new genomes with 1991 publicly available P. aeruginosa genomes for a phylogenomic and comparative analysis. We specifically explored to what extent antimicrobial resistance (AMR) genes, defence systems, and virulence genes vary in their distribution across regions of genome plasticity (RGPs) and "masked" (RGP-free) genomes, and to what extent this variation differs among the phylogroups.

FINDINGS

We found that members of phylogroup B possess larger genomes, contribute a comparatively larger number of pangenome families, and show lower abundance of CRISPR-Cas systems. Furthermore, AMR and defence systems are pervasive in RGPs and integrative and conjugative/mobilizable elements (ICEs/IMEs) from phylogroups A and B, and the abundance of these cargo genes is often significantly correlated. Moreover, inter- and intra-phylogroup interactions occur at the accessory genome level, suggesting frequent recombination events. Finally, we provide here the mPact panel of diverse P. aeruginosa strains that may serve as a valuable reference for functional analyses.

INTERPRETATION

Altogether, our results highlight distinct pangenome characteristics of the P. aeruginosa phylogroups, which are possibly influenced by variation in the abundance of CRISPR-Cas systems and are shaped by the differential distribution of other defence systems and AMR genes.

FUNDING

German Science Foundation, Max-Planck Society, Leibniz ScienceCampus Evolutionary Medicine of the Lung, BMBF program Medical Infection Genomics, Kiel Life Science Postdoc Award.

Sympatric genetic divergence between early- and late-season weedy rice populations

Sympatric genetic divergence is the most appealing and controversial pattern in the theory of ecological speciation. Examples that support sympatric genetic divergence in plant species are extremely rare. Solid evidence of sympatric genetic divergence will provide deep insights for revealing the underlying mechanisms of ecological speciation. We analysed the total genomic DNA sequences of 120 weedy rice (WR; Oryza sativa f. spontanea) plants, representing three WR population pairs separately from three early- and late-season rice fields, in comparison with those of the co-occurring rice cultivars and other rice materials. We detected substantial genetic divergence within the pairs of the sympatric early- and late-season WR populations, although genetic divergence was unevenly distributed across the genomes. Restricted gene flow was determined between the sympatric WR populations, resulting in their distinct genetic structures. We also detected relatively low genetic diversity that was likely to be associated with stronger selection in early-season WR populations. Our findings provide strong evidence for sympatric genetic divergence between the WR populations in the same fields but in different seasons. We conclude that temporal isolation plays an important role in creating genetic divergence between sympatric populations/species in plants.

N limit as a switch node between positive and negative plant-soil feedback: A meta-analysis based on the covariant diagnosis of plant growth and soil factors

Mounting evidence has confirmed the existence of plant-soil feedback, a reflection of plant-soil interaction. However, analysis of ecological feedback pathways remains a challenge. In this study, single and mixed plant communities in different soil ecosystems were screened using strict control systems in global ecosystems to identify the positive or negative feedback effects in indicator plants. Furthermore, the plant components and biomass were identified in each pathway. The significantly changed components indicated pathway factors. As negative feedback increased, the InRR (Response Ratio) of soil organic matter, soil total N, microbial alpha diversity and the symbiotic fungi proportion were significantly up-regulated (P < 0.05). In contrast, the stoichiometric ratio (C: N), water content, and the pathogenic bacteria proportion were downregulated (P < 0.05). However, the positive feedback showed the opposite trend. Importantly, N limit as a transform node between positive and negative plant-soil feedback predicted by Akaike information criterion (AIC > 0.8). Therefore, it has become an important evaluation standard for the inter-species relationship and ecological environment changes under the background of global N deposition. Finally, the feedback values of each sampling site were recalculated over the next 20 years, 50 years, and 100 years based on the global temperature rise and changing rainfall patterns. We also found that global warming and extreme rainfall may change the distribution of interspecies relationships on a global scale, with global warming having the greatest recognisable effect and decreasing the negative feedback layout by 21.7% (P < 0.05). Therefore, this work promotes the cognition of relationship of soil environment, microbial abundance and function, plant diversity and plant- soil feedback model. Meanwhile, it is of great significance to protect species diversity and restore environmental degradation.

Inter-cluster competition and resource partitioning may govern the ecology of Frankia

Microbes live in a complex communal ecosystem. The structural complexity of microbial community reflects diversity, functionality, as well as habitat type. Delineation of ecologically important microbial populations along with exploration of their roles in environmental adaptation or host-microbe interaction has a crucial role in modern microbiology. In this scenario, reverse ecology (the use of genomics to study ecology) plays a pivotal role. Since the co-existence of two different genera in one small niche should maintain a strict direct interaction, it will be interesting to utilize the concept of reverse ecology in this scenario. Here, we exploited an 'R' package, the RevEcoR, to resolve the issue of co-existing microbes which are proven to be a crucial tool for identifying the nature of their relationship (competition or complementation) persisting among them. Our target organism here is Frankia, a nitrogen-fixing actinobacterium popular for its genetic and host-specific nature. According to their plant host, Frankia has already been sub-divided into four clusters C-I, C-II, C-III, and C-IV. Our results revealed a strong competing nature of CI Frankia. Among the clusters of Frankia studied, the competition index between C-I and C-III was the largest. The other interesting result was the co-occurrence of C-II and C-IV groups. It was revealed that these two groups follow the theory of resource partitioning in their lifestyle. Metabolic analysis along with their differential transporter machinery validated our hypothesis of resource partitioning among C-II and C-IV groups.

Taxonomy of Rhizobiaceae revisited: proposal of a new framework for genus delimitation

The alphaproteobacterial family Rhizobiaceae is highly diverse, with 168 species with validly published names classified into 17 genera with validly published names. Most named genera in this family are delineated based on genomic relatedness and phylogenetic relationships, but some historically named genera show inconsistent distribution and phylogenetic breadth. The most problematic is Rhizobium , which is notorious for being highly paraphyletic, as most newly described species in the family are assigned to this genus without consideration of their proximity to existing genera, or the need to create novel genera. Moreover, many Rhizobiaceae genera lack synapomorphic traits that would give them biological and ecological significance. We propose a common framework for genus delimitation within the family Rhizobiaceae , wherein genera are defined as monophyletic groups in a core-genome gene phylogeny, that are separated from related species using a pairwise core-proteome average amino acid identity (cpAAI) threshold of approximately 86 %. We further propose that additional genomic or phenotypic evidence can justify division of species into separate genera even if they share greater than 86 % cpAAI. Applying this framework, we propose to reclassify Rhizobium rhizosphaerae and Rhizobium oryzae into Xaviernesmea gen. nov. Data is also provided to support the formation of Peteryoungia aggregata comb. nov., Endobacterium yantingense comb. nov., Neorhizobium petrolearium comb. nov., Pararhizobium arenae comb. nov., Pseudorhizobium tarimense comb. nov. and Mycoplana azooxidifex comb. nov. Lastly, we present arguments that the unification of the genera Ensifer and Sinorhizobium in Opinion 84 of the Judicial Commission is no longer justified by current genomic and phenotypic data. Despite pairwise cpAAI values for all Ensifer species and all Sinorhizobium species being >86 %, additional genomic and phenotypic data suggest that they significantly differ in their biology and ecology. We therefore propose emended descriptions of Ensifer and Sinorhizobium , which we argue should be considered as separate genera.

Habitat Adaptation Drives Speciation of a Streptomyces Species with Distinct Habitats and Disparate Geographic Origins

Microbial diversification is driven by geographic and ecological factors, but how the relative importance of these factors varies among species, geographic scales, and habitats remains unclear. Streptomyces, a genus of antibiotic-producing, spore-forming, and widespread bacteria, offers a robust model for identifying the processes underlying population differentiation. We examined the population structure of 37 Streptomyces olivaceus strains isolated from various sources, showing that they diverged into two habitat-associated (free-living and insect-associated) and geographically disparate lineages. More frequent gene flow within than between the lineages confirmed genetic isolation in S. olivaceus. Geographic isolation could not explain the genetic isolation; instead, habitat type was a strong predictor of genetic distance when controlling for geographic distance. The identification of habitat-specific genetic variations, including genes involved in regulation, resource use, and secondary metabolism, suggested a significant role of habitat adaptation in the diversification process. Physiological assays revealed fitness trade-offs under different environmental conditions in the two lineages. Notably, insect-associated isolates could outcompete free-living isolates in a free-iron-deficient environment. Furthermore, substrate (e.g., sialic acid and glycogen) utilization but not thermal traits differentiated the two lineages. Overall, our results argue that adaptive processes drove ecological divergence among closely related streptomycetes, eventually leading to dispersal limitation and gene flow barriers between the lineages. S. olivaceus may best be considered a species complex consisting of two cryptic species. IMPORTANCE Both isolation by distance and isolation by environment occur in bacteria, and different diversification patterns may apply to different species. Streptomyces species, typified by producing useful natural products, are widespread in nature and possess high genetic diversity. However, the ecological processes and evolutionary mechanisms that shape their distribution are not well understood. Here, we show that the population structure of a ubiquitous Streptomyces species complex matches its habitat distribution and can be defined by gene flow discontinuities. Using comparative genomics and physiological assays, we reveal that gains and losses of specific genomic traits play a significant role in the transition between free-living and host-associated lifestyles, driving speciation of the species. These results provide new insights into the evolutionary trajectory of Streptomyces and the notion of species.

Different evolutionary trends form the twilight zone of the bacterial pan-genome

The pan-genome is defined as the combined set of all genes in the gene pool of a species. Pan-genome analyses have been very useful in helping to understand different evolutionary dynamics of bacterial species: an open pan-genome often indicates a free-living lifestyle with metabolic versatility, while closed pan-genomes are linked to host-restricted, ecologically specialized bacteria. A detailed understanding of the species pan-genome has also been instrumental in tracking the phylodynamics of emerging drug resistance mechanisms and drug-resistant pathogens. However, current approaches to analyse a species’ pan-genome do not take the species population structure into account, nor do they account for the uneven sampling of different lineages, as is commonplace due to over-sampling of clinically relevant representatives. Here we present the application of a population structure-aware approach for classifying genes in a pan-genome based on within-species distribution. We demonstrate our approach on a collection of 7500 Escherichia coli genomes, one of the most-studied bacterial species and used as a model for an open pan-genome. We reveal clearly distinct groups of genes, clustered by different underlying evolutionary dynamics, and provide a more biologically informed and accurate description of the species’ pan-genome.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Eco-evolutionary feedbacks mediated by bacterial membrane vesicles

Bacterial membrane vesicles (BMVs) are spherical extracellular organelles whose cargo is enclosed by a biological membrane. The cargo can be delivered to distant parts of a given habitat in a protected and concentrated manner. This review presents current knowledge about BMVs in the context of bacterial eco-evolutionary dynamics among different environments and hosts. BMVs may play an important role in establishing and stabilizing bacterial communities in such environments; for example, bacterial populations may benefit from BMVs to delay the negative effect of certain evolutionary trade-offs that can result in deleterious phenotypes. BMVs can also perform ecosystem engineering by serving as detergents, mediators in biochemical cycles, components of different biofilms, substrates for cross-feeding, defense systems against different dangers and enzyme-delivery mechanisms that can change substrate availability. BMVs further contribute to bacteria as mediators in different interactions, with either other bacterial species or their hosts. In short, BMVs extend and deliver phenotypic traits that can have ecological and evolutionary value to both their producers and the ecosystem as a whole.

Phylogenomics reveals the basis of adaptation of Pseudorhizobium species to extreme environments and supports a taxonomic revision of the genus

The family Rhizobiaceae includes many genera of soil bacteria, often isolated for their association with plants. Herein, we investigate the genomic diversity of a group of Rhizobium species and unclassified strains isolated from atypical environments, including seawater, rock matrix or polluted soil. Based on whole-genome similarity and core genome phylogeny, we show that this group corresponds to the genus Pseudorhizobium. We thus reclassify Rhizobium halotolerans, R. marinum, R. flavum and R. endolithicum as P. halotolerans sp. nov., P. marinum comb. nov., P. flavum comb. nov. and P. endolithicum comb. nov., respectively, and show that P. pelagicum is a synonym of P. marinum. We also delineate a new chemolithoautotroph species, P. banfieldiae sp. nov., whose type strain is NT-26^T (=DSM 106348^T=CFBP 8663^T). This genome-based classification was supported by a chemotaxonomic comparison, with increasing taxonomic resolution provided by fatty acid, protein and metabolic profiles. In addition, we used a phylogenetic approach to infer scenarios of duplication, horizontal transfer and loss for all genes in the Pseudorhizobium pangenome. We thus identify the key functions associated with the diversification of each species and higher clades, shedding light on the mechanisms of adaptation to their respective ecological niches. Respiratory proteins acquired at the origin of Pseudorhizobium were combined with clade-specific genes to enable different strategies for detoxification and nutrition in harsh, nutrient-poor environments.

Co-occurrence of rhizobacteria with nitrogen fixation and/or 1-aminocyclopropane-1-carboxylate deamination abilities in the maize rhizosphere

The plant microbiota may differ depending on soil type, but these microbiota probably share the same functions necessary for holobiont fitness. Thus, we tested the hypothesis that phytostimulatory microbial functional groups are likely to co-occur in the rhizosphere, using groups corresponding to nitrogen fixation (nifH) and 1-aminocyclopropane-1-carboxylate deamination (acdS), i.e. two key modes of action in plant-beneficial rhizobacteria. The analysis of three maize fields in two consecutive years showed that quantitative PCR numbers of nifH and of acdS alleles differed according to field site, but a positive correlation was found overall when comparing nifH and acdS numbers. Metabarcoding analyses in the second year indicated that the diversity level of acdS but not nifH rhizobacteria in the rhizosphere differed across fields. Furthermore, between-class analysis showed that the three sites differed from one another based on nifH or acdS sequence data (or rrs data), and the bacterial genera contributing most to field differentiation were not the same for the three bacterial groups. However, co-inertia analysis indicated that the genetic structures of both functional groups and of the whole bacterial community were similar across the three fields. Therefore, results point to co-selection of rhizobacteria harboring nitrogen fixation and/or 1-aminocyclopropane-1-carboxylate deamination abilities.

Horizontal Gene Transfer as a Source of Conflict and Cooperation in Prokaryotes

Horizontal gene transfer (HGT) is one of the most important processes in prokaryote evolution. The sharing of DNA can spread neutral or beneficial genes, as well as genetic parasites across populations and communities, creating a large proportion of the variability acted on by natural selection. Here, we highlight the role of HGT in enhancing the opportunities for conflict and cooperation within and between prokaryote genomes. We discuss how horizontally acquired genes can cooperate or conflict both with each other and with a recipient genome, resulting in signature patterns of gene co-occurrence, avoidance, and dependence. We then describe how interactions involving horizontally transferred genes may influence cooperation and conflict at higher levels (populations, communities, and symbioses). Finally, we consider the benefits and drawbacks of HGT for prokaryotes and its fundamental role in understanding conflict and cooperation from the gene-gene to the microbiome level.

Ancestral Genome Estimation Reveals the History of Ecological Diversification in Agrobacterium

Horizontal gene transfer (HGT) is considered as a major source of innovation in bacteria, and as such is expected to drive adaptation to new ecological niches. However, among the many genes acquired through HGT along the diversification history of genomes, only a fraction may have actively contributed to sustained ecological adaptation. We used a phylogenetic approach accounting for the transfer of genes (or groups of genes) to estimate the history of genomes in Agrobacterium biovar 1, a diverse group of soil and plant-dwelling bacterial species. We identified clade-specific blocks of cotransferred genes encoding coherent biochemical pathways that may have contributed to the evolutionary success of key Agrobacterium clades. This pattern of gene coevolution rejects a neutral model of transfer, in which neighboring genes would be transferred independently of their function and rather suggests purifying selection on collectively coded acquired pathways. The acquisition of these synapomorphic blocks of cofunctioning genes probably drove the ecological diversification of Agrobacterium and defined features of ancestral ecological niches, which consistently hint at a strong selective role of host plant rhizospheres.

Rapidly evolving changes and gene loss associated with host switching in Corynebacterium pseudotuberculosis

Phylogenomics and genome scale positive selection analyses were performed on 29 Corynebacterium pseudotuberculosis genomes that were isolated from different hosts, including representatives of the Ovis and Equi biovars. A total of 27 genes were identified as undergoing adaptive changes. An analysis of the clades within this species and these biovars, the genes specific to each branch, and the genes responding to selective pressure show clear differences, indicating that adaptation and specialization is occurring in different clades. These changes are often correlated with the isolation host but could indicate responses to some undetermined factor in the respective niches. The fact that some of these more-rapidly evolving genes have homology to known virulence factors, antimicrobial resistance genes and drug targets shows that this type of analysis could be used to identify novel targets, and that these could be used as a way to control this pathogen.

Niche Construction and Exploitation by Agrobacterium: How to Survive and Face Competition in Soil and Plant Habitats

Agrobacterium populations live in different habitats (bare soil, rhizosphere, host plants), and hence face different environmental constraints. They have evolved the capacity to exploit diverse resources and to escape plant defense and competition from other microbiota. By modifying the genome of their host, Agrobacterium populations exhibit the remarkable ability to construct and exploit the ecological niche of the plant tumors that they incite. This niche is characterized by the accumulation of specific, low molecular weight compounds termed opines that play a critical role in Agrobacterium 's lifestyle. We present and discuss the functions, advantages, and costs associated with this niche construction and exploitation.

Delineation of Taxonomic Species within Complex of Species: Aeromonas media and Related Species as a Test Case

Aeromonas media is an opportunistic pathogen for human and animals mainly found in aquatic habitats and which has been noted for significant genomic and phenotypic heterogeneities. We aimed to better understand the population structure and diversity of strains currently affiliated to A. media and the related species A. rivipollensis. Forty-one strains were included in a population study integrating, multilocus genetics, phylogenetics, comparative genomics, as well as phenotypics, lifestyle, and evolutionary features. Sixteen gene-based multilocus phylogeny delineated three clades. Clades corresponded to different genomic groups or genomospecies defined by phylogenomic metrics ANI (average nucleotide identity) and isDDH (in silico DNA-DNA hybridization) on 14 whole genome sequences. DL-lactate utilization, cefoxitin susceptibility, nucleotide signatures, ribosomal multi-operon diversity, and differences in relative effect of recombination and mutation (i.e., in evolution mode) distinguished the two species Aeromonas media and Aeromonas rivipollensis. The description of these two species was emended accordingly. The genome metrics and comparative genomics suggested that a third clade is a distinct genomospecies. Beside the species delineation, genetic and genomic data analysis provided a more comprehensive knowledge of the cladogenesis determinants at the root and inside A. media species complex among aeromonads. Particular lifestyles and phenotypes as well as major differences in evolution modes may represent putative factors associated with lineage emergence and speciation within the A. media complex. Finally, the integrative and populational approach presented in this study is considered broadly in order to conciliate the delineation of taxonomic species and the population structure in bacterial genera organized in species complexes.

Evolutionary study of Yersinia genomes deciphers emergence of human pathogenic species

On record, there are 17 species in the Yersinia genus, of which three are known to be pathogenic to human. While the chromosomal and pYV (or pCD1) plasmid-borne virulence genes as well as pathogenesis of these three species are well studied, their genomic evolution is poorly understood. Our study aims to predict the key evolutionary events that led to the emergence of pathogenic Yersinia species by analyzing gene gain-and-loss, virulence genes, and “Clustered regularly-interspaced short palindromic repeats”. Our results suggest that the most recent ancestor shared by the human pathogenic Yersinia was most probably an environmental species that had adapted to the human body. This might have led to ecological specialization that diverged Yersinia into ecotypes and distinct lineages based on differential gene gain-and-loss in different niches. Our data also suggest that Y. pseudotuberculosis group might be the donor of the ail virulence gene to Y. enterocolitica. Hence, we postulate that evolution of human pathogenic Yersinia might not be totally in parallel, but instead, there were lateral gene transfer events. Furthermore, the presence of virulence genes seems to be important for the positive selection of virulence plasmid. Our studies provide better insights into the evolutionary biology of these bacteria.

Pangenome Evolution in the Marine Bacterium Alteromonas

We have examined a collection of the free-living marine bacterium Alteromonas genomes with cores diverging in average nucleotide identities ranging from 99.98% to 73.35%, i.e., from microbes that can be considered members of a natural clone (like in a clinical epidemiological outbreak) to borderline genus level. The genomes were largely syntenic allowing a precise delimitation of the core and flexible regions in each. The core was 1.4 Mb (ca. 30% of the typical strain genome size). Recombination rates along the core were high among strains belonging to the same species (37.7-83.7% of all nucleotide polymorphisms) but they decreased sharply between species (18.9-5.1%). Regarding the flexible genome, its main expansion occurred within the boundaries of the species, i.e., strains of the same species already have a large and diverse flexible genome. Flexible regions occupy mostly fixed genomic locations. Four large genomic islands are involved in the synthesis of strain-specific glycosydic receptors that we have called glycotypes. These genomic regions are exchanged by homologous recombination within and between species and there is evidence for their import from distant taxonomic units (other genera within the family). In addition, several hotspots for integration of gene cassettes by illegitimate recombination are distributed throughout the genome. They code for features that give each clone specific properties to interact with their ecological niche and must flow fast throughout the whole genus as they are found, with nearly identical sequences, in different species. Models for the generation of this genomic diversity involving phage predation are discussed.

What Is Speciation?

Concepts and definitions of species have been debated by generations of biologists and remain controversial. Microbes pose a particular challenge because of their genetic diversity, asexual reproduction, and often promiscuous horizontal gene transfer (HGT). However, microbes also present an opportunity to study and understand speciation because of their rapid evolution, both in nature and in the lab, and small, easily sequenced genomes. Here, we review how microbial population genomics has enabled us to catch speciation "in the act" and how the results have challenged and enriched our concepts of species, with implications for all domains of life. We describe how recombination (including HGT and introgression) has shaped the genomes of nascent microbial, animal, and plant species and argue for a prominent role of natural selection in initiating and maintaining speciation. We ask how universal is the process of speciation across the tree of life, and what lessons can be drawn from microbes? Comparative genomics showing the extent of HGT in natural populations certainly jeopardizes the relevance of vertical descent (i.e., the species tree) in speciation. Nevertheless, we conclude that species do indeed exist as clusters of genetic and ecological similarity and that speciation is driven primarily by natural selection, regardless of the balance between horizontal and vertical descent.