Temporal fragmentation of speciation in bacteria

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.

Identification of new SdiA regulon members of Escherichia coli, Enterobacter cloacae, and Salmonella enterica serovars Typhimurium and Typhi

Bacteria can coordinate behavior in response to population density through the production, release, and detection of small molecules, a phenomenon known as quorum sensing. Salmonella enterica is among a group of Enterobacteriaceae that can detect signaling molecules of the N-acyl homoserine lactone (AHL) type but lack the ability to produce them. The AHLs are detected by the LuxR-type transcription factor, SdiA. This enables a behavior known as eavesdropping, where organisms can sense the signaling molecules of other species of bacteria. The role of SdiA remains largely unknown. Here, we use RNA-seq to more completely identify the sdiA regulons of two clinically significant serovars of Salmonella enterica: Typhimurium and Typhi. We find that their sdiA regulons are largely conserved despite the significant differences in pathogenic strategy and host range of these two serovars. Previous studies identified sdiA-regulated genes in Escherichia coli and Enterobacter cloacae, but there is surprisingly no overlap in regulon membership between the different species. This led us to individually test orthologs of each regulon member in the other species and determine that there is indeed some overlap. Unfortunately, the functions of most sdiA-regulated genes are unknown, with the overall function of eavesdropping in these organisms remaining unclear.

IMPORTANCE

Many bacterial species detect their own population density through the production, release, and detection of small molecules (quorum sensing). Salmonella and other Enterobacteriaceae have a modified system that detects the N-acyl-homoserine lactones of other bacteria through the solo quorum sensing receptor SdiA, a behavior known as eavesdropping. The roles of sdiA-dependent eavesdropping in the lifecycles of these bacteria are unknown. In this study, we identify sdiA-dependent transcriptional responses in two clinically relevant serovars of Salmonella, Typhimurium and Typhi, and note that their responses are partially conserved. We also demonstrate for the first time that sdiA-dependent regulation of genes is partially conserved in Enterobacter cloacae and Escherichia coli as well, indicating a degree of commonality in eavesdropping among the Enterobacteriaceae.

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.

Using the pan-genomic framework for the discovery of genomic islands in the haloarchaeon Halorubrum ezzemoulense

In this study, we use pan-genomics to characterize the genomic variability of the widely dispersed halophilic archaeal species Halorubrum ezzemoulense (Hez). We include a multi-regional sampling of newly sequenced, high-quality draft genomes. The pan-genome graph of the species reveals 50 genomic islands that represent rare accessory genetic capabilities available to members. Most notably, we observe rearrangements that have led to the insertion/recombination/replacement of mutually exclusive genomic islands in equivalent genome positions ("homeocassettes"). These conflicting islands encode for similar functions, but homologs from islands located between the same core genes exhibit high divergence on the amino acid level, while the neighboring core genes are nearly identical. Both islands of a homeocassette often coexist in the same geographic location, suggesting that either island may be beyond the reach of selective sweeps and that these loci of divergence between Hez members are maintained and persist long term. This implies that subsections of the population have different niche preferences and rare metabolic capabilities. After an evaluation of the gene content in the homeocassettes, we speculate that these islands may play a role in the speciation, niche adaptability, and group selection dynamics in Hez. Though homeocassettes are first described in this study, similar replacements and divergence of genes on genomic islands have been previously reported in other Haloarchaea and distantly related Archaea, suggesting that homeocassettes may be a feature in a wide range of organisms outside of Hez.IMPORTANCEThis study catalogs the rare genes discovered in strains of the species Halorubrum ezzemoulense (Hez), an obligate halophilic archaeon, through the perspective of its pan-genome. These rare genes are often found to be arranged on islands that confer metabolic and transport functions and contain genes that have eluded previous studies. The discovery of divergent, but homologous islands occupying equivalent genome positions ("homeocassettes") in different genomes, reveals significant new information on genome evolution in Hez. Homeocassette pairs encode for similar functions, but their dissimilarity and distribution imply high rates of recombination, different specializations, and niche preferences in Hez. The coexistence of both islands of a homeocassette pair in multiple environments demonstrates that both islands are beyond the reach of selective sweeps and that these genome content differences between strains persist long term. The switch between islands through recombination under different environmental conditions may lead to a greater range of niche adaptability in Hez.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Complete genome sequences and characteristics of mycobacteriophages Diminimus, Dulcita, Glaske16, and Koreni

Complete genome sequences of four novel mycobacteriophages, Diminimus, Dulcita, Glaske16, and Koreni, isolated from soil are presented. All these bacteriophages belong to subcluster M1, except Koreni that belongs to subcluster A4. Moreover, all have siphovirus morphologies, with genome sizes ranging from 51,055 to 81,156 bp.

Impact of lead (Pb2+) on the growth and biological activity of Serratia marcescens selected for wastewater treatment and identification of its zntR gene-a metal efflux regulator

Microorganisms isolated from contaminated areas play an important role in bioremediation processes. They promote heavy metal removal from the environment by adsorbing ions onto the cell wall surface, accumulating them inside the cells, or reducing, complexing, or precipitating these substances in the environment. Microorganism-based bioremediation processes can be highly efficient, low-cost and have low environmental impact. Thus, the present study aimed to select Pb²⁺-resistant bacteria and evaluate the growth rate, biological activity, and the presence of genes associated with metal resistance. Serratia marcescens CCMA 1010, that was previously isolated from coffee processing wastewater, was selected since was able to growth in Pb²⁺ concentrations of up to 4.0 mM. The growth rate and generation time did not differ from those of the control (without Pb²⁺), although biological activity decreased in the first hour of exposure to these ions and stabilized after this period. The presence of the zntR, zntA and pbrA genes was analysed, and only zntR was detected. The zntR gene encodes a protein responsible for regulating the production of ZntA, a transmembrane protein that facilitates Pb²⁺ extrusion out of the cell. S. marcescens CCMA 1010 demonstrated a potential for use as bioindicator that has potential to be used in bioremediation processes due to its resistance to high concentrations of Pb²⁺, ability to grow until 24 h of exposure, and possession of a gene that indicates the existence of mechanisms associated with resistance to lead (Pb²⁺).

The role of hand size in body representation: a developmental investigation

Knowledge of one's own body size is a crucial facet of body representation, both for acting on the environment and perhaps also for constraining body ownership. However, representations of body size may be somewhat plastic, particularly to allow for physical growth in childhood. Here we report a developmental investigation into the role of hand size in body representation (the sense of body ownership, perception of hand position, and perception of own-hand size). Using the rubber hand illusion paradigm, this study used different fake hand sizes (60%, 80%, 100%, 120% or 140% of typical size) in three age groups (6- to 7-year-olds, 12- to 13-year-olds, and adults; N = 229). We found no evidence that hand size constrains ownership or position: participants embodied hands which were both larger and smaller than their own, and indeed judged their own hands to have changed size following the illusion. Children and adolescents embodied the fake hands more than adults, with a greater tendency to feel their own hand had changed size. Adolescents were particularly sensitive to multisensory information. In sum, we found substantial plasticity in the representation of own-body size, with partial support for the hypothesis that children have looser representations than adults.

Diversity and Evolution of Pigment Types in Marine Synechococcus Cyanobacteria

Synechococcus cyanobacteria are ubiquitous and abundant in the marine environment and contribute to an estimated 16% of the ocean net primary productivity. Their light-harvesting complexes, called phycobilisomes (PBS), are composed of a conserved allophycocyanin core, from which radiates six to eight rods with variable phycobiliprotein and chromophore content. This variability allows Synechococcus cells to optimally exploit the wide variety of spectral niches existing in marine ecosystems. Seven distinct pigment types or subtypes have been identified so far in this taxon based on the phycobiliprotein composition and/or the proportion of the different chromophores in PBS rods. Most genes involved in their biosynthesis and regulation are located in a dedicated genomic region called the PBS rod region. Here, we examine the variability of gene content and organization of this genomic region in a large set of sequenced isolates and natural populations of Synechococcus representative of all known pigment types. All regions start with a tRNA-PheGAA and some possess mobile elements for DNA integration and site-specific recombination, suggesting that their genomic variability relies in part on a "tycheposon"-like mechanism. Comparison of the phylogenies obtained for PBS and core genes revealed that the evolutionary history of PBS rod genes differs from the core genome and is characterized by the co-existence of different alleles and frequent allelic exchange. We propose a scenario for the evolution of the different pigment types and highlight the importance of incomplete lineage sorting in maintaining a wide diversity of pigment types in different Synechococcus lineages despite multiple speciation events.

Niche Preference of Escherichia coli in a Peri-Urban Pond Ecosystem

Escherichia coli comprises diverse strains with a large accessory genome, indicating functional diversity and the ability to adapt to a range of niches. Specific strains would display greatest fitness in niches matching their combination of phenotypic traits. Given this hypothesis, we sought to determine whether E. coli in a peri-urban pond and associated cattle pasture display niche preference. Samples were collected from water, sediment, aquatic plants, water snails associated with the pond, as well as bovine feces from cattle in an adjacent pasture. Isolates (120) were obtained after plating on Membrane Lactose Glucuronide Agar (MLGA). We used the uidA and mutS sequences for all isolates to determine phylogeny by maximum likelihood, and population structure through gene flow analysis. PCR was used to allocate isolates to phylogroups and to determine the presence of pathogenicity/virulence genes (stxI, stxII, eaeA, hlyA, ST, and LT). Antimicrobial resistance was determined using a disk diffusion assay for Tetracycline, Gentamicin, Ciprofloxacin, Meropenem, Ceftriaxone, and Azithromycin. Our results showed that isolates from water, sediment, and water plants were similar by phylogroup distribution, virulence gene distribution, and antibiotic resistance while both snail and feces populations were significantly different. Few of the feces isolates were significantly similar to aquatic ones, and most of the snail isolates were also different. Population structure analysis indicated three genetic backgrounds associated with bovine, snail, and aquatic environments. Collectively these data support niche preference of E. coli isolates occurring in this ecosystem.

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

From sequence to information

Today massive amounts of sequenced metagenomic and metatranscriptomic data from different ecological niches and environmental locations are available. Scientific progress depends critically on methods that allow extracting useful information from the various types of sequence data. Here, we will first discuss types of information contained in the various flavours of biological sequence data, and how this information can be interpreted to increase our scientific knowledge and understanding. We argue that a mechanistic understanding of biological systems analysed from different perspectives is required to consistently interpret experimental observations, and that this understanding is greatly facilitated by the generation and analysis of dynamic mathematical models. We conclude that, in order to construct mathematical models and to test mechanistic hypotheses, time-series data are of critical importance. We review diverse techniques to analyse time-series data and discuss various approaches by which time-series of biological sequence data have been successfully used to derive and test mechanistic hypotheses. Analysing the bottlenecks of current strategies in the extraction of knowledge and understanding from data, we conclude that combined experimental and theoretical efforts should be implemented as early as possible during the planning phase of individual experiments and scientific research projects.

This article is part of the theme issue ‘Integrative research perspectives on marine conservation’.

Diversity within species: interpreting strains in microbiomes

Studying within-species variation has traditionally been limited to culturable bacterial isolates and low-resolution microbial community fingerprinting. Metagenomic sequencing and technical advances have enabled culture-free, high-resolution strain and subspecies analyses at high throughput and in complex environments. This holds great scientific promise but has also led to an overwhelming number of methods and terms to describe infraspecific variation. This Review aims to clarify these advances by focusing on the diversity within bacterial and archaeal species in the context of microbiomics. We cover foundational microevolutionary concepts relevant to population genetics and summarize how within-species variation can be studied and stratified directly within microbial communities with a focus on metagenomics. Finally, we describe how common applications of within-species variation can be achieved using metagenomic data. We aim to guide the selection of appropriate terms and analytical approaches to facilitate researchers in benefiting from the increasing availability of large, high-resolution microbiome genetic sequencing data.

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.

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.

Ecological Divergence with Gene Flow in a Thermophilic Cyanobacterium

How ecological diversity is maintained and distributed within populations is a longstanding question in microbial ecology. In the thermophilic cyanobacterium Synechococcus B', high observed levels of recombination are predicted to maintain ecological variation despite the simultaneous action of diverse selective pressures on different regions of the genome. To investigate ecological diversity in these bacteria, we directly isolated laboratory strains of Synechococcus B' from samples collected along the thermal gradients of two geothermal environments in Yellowstone National Park. Extensive recombination was evident for a multi-locus sequence data set, and, consequently, our sample did not exhibit the sequence clustering expected for distinct ecotypes evolving by periodic clonal selection. Evidence for local selective sweeps at specific loci suggests that sweeps may be common but that recombination is effective for maintaining diversity of unlinked genomic regions. Thermal performance for strain growth was positively associated with the temperature of the environment, indicating that Synechococcus B' populations consist of locally adapted ecological specialists that occupy specific thermal niches. Because this ecological differentiation is observed despite the absence of dispersal barriers among sites, we conclude that these bacteria may freely exchange much of the genome but that barriers to gene flow exist for loci under direct temperature selection.

Delineation of ecologically distinct units of marine Bacteroidetes in the Northwestern Mediterranean Sea

Bacteroidetes is one of the dominant phyla of ocean bacterioplankton, yet its diversity and population structure is poorly understood. To advance in the delineation of ecologically meaningful units within this group, we constructed near full-length 16S rRNA gene clone libraries from contrasting marine environments in the NW Mediterranean. Based on phylogeny and the associated ecological variables (depth and season), 24 different Bacteroidetes clades were delineated. By considering their relative abundance (from iTag amplicon sequencing studies), we described the distribution patterns of each of these clades, delimiting them as Ecologically Significant Taxonomic Units (ESTUs). Spatially, there was almost no overlap among ESTUs at different depths. In deep waters there was predominance of Owenweeksia, Leeuwenhoekiella, Muricauda-related genera, and some depth-associated ESTUs within the NS5 and NS2b marine clades. Seasonally, multi-annual dynamics of recurring ESTUs were present with dominance of some ESTUs within the NS4, NS5 and NS2b marine clades along most of the year, but with variable relative frequencies between months. A drastic change towards the predominance of Formosa-related ESTUs and one ESTU from the NS5 marine clade was typically present after the spring bloom. Even though there are no isolates available for these ESTUs to determine their physiology, correlation analyses identified the environmental preference of some of them. Overall, our results suggest that there is a high degree of niche specialisation within these closely related clades. This work constitutes a step forward in disentangling the ecology of marine Bacteroidetes, which are essential players in organic matter processing in the oceans.

Marker genes as predictors of shared genomic function

Background

Although high-throughput marker gene studies provide valuable insight into the diversity and relative abundance of taxa in microbial communities, they do not provide direct measures of their functional capacity. Recently, scientists have shown a general desire to predict functional profiles of microbial communities based on phylogenetic identification inferred from marker genes, and recent tools have been developed to link the two. However, to date, no large-scale examination has quantified the correlation between the marker gene based taxonomic identity and protein coding gene conservation. Here we utilize 4872 representative prokaryotic genomes from NCBI to investigate the relationship between marker gene identity and shared protein coding gene content.

Results

Even at 99–100% marker gene identity, genomes share on average less than 75% of their protein coding gene content. This occurs regardless of the marker gene(s) used: V4 region of the 16S rRNA, complete 16S rRNA, or single copy orthologs through a multi-locus sequence analysis. An important aspect related to this observation is the intra-organism variation of 16S copies from a single genome. Although the majority of 16S copies were found to have high sequence similarity (> 99%), several genomes contained copies that were highly diverged (< 97% identity).

Conclusions

This is the largest comparison between marker gene similarity and shared protein coding gene content to date. The study highlights the limitations of inferring a microbial community’s functions based on marker gene phylogeny. The data presented expands upon the results of previous studies that examined one or few bacterial species and supports the hypothesis that 16S rRNA and other marker genes cannot be directly used to fully predict the functional potential of a bacterial community.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5641-1) contains supplementary material, which is available to authorized users.

The role of recombination, niche-specific gene pools and flexible genomes in the ecological speciation of bacteria

Bacteria diversify into genetic clusters analogous to those observed in sexual eukaryotes, but the definition of bacterial species is an ongoing problem. Recent work has focused on adaptation to distinct ecological niches as the main driver of clustering, but there remains debate about the role of recombination in that process. One view is that homologous recombination occurs too rarely for gene flow to constrain divergent selection. Another view is that homologous recombination is frequent enough in many bacterial populations that barriers to gene flow are needed to permit divergence. Niche-specific gene pools have been proposed as a general mechanism to limit gene flow. We use theoretical models to evaluate additional hypotheses that evolving genetic architecture, specifically the effect sizes of genes and gene gain and loss, can limit gene flow between diverging populations. Our model predicts that (a) in the presence of gene flow and recombination, ecological divergence is concentrated in few loci of large effect and (b) high rates of gene flow plus recombination promote gene loss and favor the evolution of niche-specific genes. The results show that changing genetic architecture and gene loss can facilitate ecological divergence, even without niche-specific gene pools. We discuss these results in the context of recent studies of sympatric divergence in microbes.

Introduction to Genome Biology and Diversity

Organisms display astonishing levels of cell and molecular diversity, including genome size, shape, and architecture. In this chapter, we review how the genome can be viewed as both a structural and an informational unit of biological diversity and explicitly define our intended meaning of genetic information. A brief overview of the characteristic features of bacterial, archaeal, and eukaryotic cell types and viruses sets the stage for a review of the differences in organization, size, and packaging strategies of their genomes. We include a detailed review of genetic elements found outside the primary chromosomal structures, as these provide insights into how genomes are sometimes viewed as incomplete informational entities. Lastly, we reassess the definition of the genome in light of recent advancements in our understanding of the diversity of genomic structures and the mechanisms by which genetic information is expressed within the cell. Collectively, these topics comprise a good introduction to genome biology for the newcomer to the field and provide a valuable reference for those developing new statistical or computation methods in genomics. This review also prepares the reader for anticipated transformations in thinking as the field of genome biology progresses.

Homologous Recombination between Genetically Divergent Campylobacter fetus Lineages Supports Host-Associated Speciation

Homologous recombination is a major driver of bacterial speciation. Genetic divergence and host association are important factors influencing homologous recombination. Here, we study these factors for Campylobacter fetus, which shows a distinct intraspecific host dichotomy. Campylobacter fetus subspecies fetus (Cff) and venerealis are associated with mammals, whereas C. fetus subsp. testudinum (Cft) is associated with reptiles. Recombination between these genetically divergent C. fetus lineages is extremely rare. Previously it was impossible to show whether this barrier to recombination was determined by the differential host preferences, by the genetic divergence between both lineages or by other factors influencing recombination, such as restriction-modification, CRISPR/Cas, and transformation systems. Fortuitously, a distinct C. fetus lineage (ST69) was found, which was highly related to mammal-associated C. fetus, yet isolated from a chelonian. The whole genome sequences of two C. fetus ST69 isolates were compared with those of mammal- and reptile-associated C. fetus strains for phylogenetic and recombination analysis. In total, 5.1–5.5% of the core genome of both ST69 isolates showed signs of recombination. Of the predicted recombination regions, 80.4% were most closely related to Cft, 14.3% to Cff, and 5.6% to C. iguaniorum. Recombination from C. fetus ST69 to Cft was also detected, but to a lesser extent and only in chelonian-associated Cft strains. This study shows that despite substantial genetic divergence no absolute barrier to homologous recombination exists between two distinct C. fetus lineages when occurring in the same host type, which provides valuable insights in bacterial speciation and evolution.

Genomic analyses of the ancestral Manila family of Mycobacterium tuberculosis

With its airborne transmission and prolonged latency period, Mycobacterium tuberculosis spreads worldwide as one of the most successful bacterial pathogens and continues to kill millions of people every year. M. tuberculosis lineage 1 is inferred to originate ancestrally based on the presence of the 52-bp TbD1 sequence and analysis of single nucleotide polymorphisms. Previously, we briefly reported the complete genome sequencing of M. tuberculosis strains 96121 and 96075, which belong to the ancient Manila family and modern Beijing family respectively. Here we present the comprehensive genomic analyses of the Manila family in lineage 1 compared to complete genomes in lineages 2-4. Principal component analysis of the presence and absence of CRISPR spacers suggests that Manila isolate 96121 is distinctly distant from lineages 2-4. We further identify a truncated whiB5 gene and a putative operon consisting of genes encoding a putative serine/threonine kinase PknH and a putative ABC transporter, which are only found in the genomes of Manila family isolates. Six single nucleotide polymorphisms are uniquely conserved in 38 Manila strains. Moreover, when compared to M. tuberculosis H37Rv, 59 proteins are under positive selection in Manila family isolate 96121 but not in Beijing family isolate 96075. The unique features further serve as biomarkers for Manila strains and may shed light on the limited transmission of this ancestral lineage outside of its Filipino host population.

Evolutionary Roots and Diversification of the Genus Aeromonas

Despite the importance of diversification rates in the study of prokaryote evolution, they have not been quantitatively assessed for the majority of microorganism taxa. The investigation of evolutionary patterns in prokaryotes constitutes a challenge due to a very scarce fossil record, limited morphological differentiation and frequently complex taxonomic relationships, which make even species recognition difficult. Although the speciation models and speciation rates in eukaryotes have traditionally been established by analyzing the fossil record data, this is frequently incomplete, and not always available. More recently, several methods based on molecular sequence data have been developed to estimate speciation and extinction rates from phylogenies reconstructed from contemporary taxa. In this work, we determined the divergence time and temporal diversification of the genus Aeromonas by applying these methods widely used with eukaryotic taxa. Our analysis involved 150 Aeromonas strains using the concatenated sequences of two housekeeping genes (approximately 2,000 bp). Dating and diversification model analyses were performed using two different approaches: obtaining the consensus sequence from the concatenated sequences corresponding to all the strains belonging to the same species, or generating the species tree from multiple alignments of each gene. We used BEAST to perform a Bayesian analysis to estimate both the phylogeny and the divergence times. A global molecular clock cannot be assumed for any gene. From the chronograms obtained, we carried out a diversification analysis using several approaches. The results suggest that the genus Aeromonas began to diverge approximately 250 millions of years (Ma) ago. All methods used to determine Aeromonas diversification gave similar results, suggesting that the speciation process in this bacterial genus followed a rate-constant (Yule) diversification model, although there is a small probability that a slight deceleration occurred in recent times. We also determined the constant of diversification (λ) values, which in all cases were very similar, about 0.01 species/Ma, a value clearly lower than those described for different eukaryotes.

Evolution, genomics and epidemiology of Pseudomonas syringae: Challenges in Bacterial Molecular Plant Pathology

A remarkable shift in our understanding of plant-pathogenic bacteria is underway. Until recently, nearly all research on phytopathogenic bacteria was focused on a small number of model strains, which provided a deep, but narrow, perspective on plant-microbe interactions. Advances in genome sequencing technologies have changed this by enabling the incorporation of much greater diversity into comparative and functional research. We are now moving beyond a typological understanding of a select collection of strains to a more generalized appreciation of the breadth and scope of plant-microbe interactions. The study of natural populations and evolution has particularly benefited from the expansion of genomic data. We are beginning to have a much deeper understanding of the natural genetic diversity, niche breadth, ecological constraints and defining characteristics of phytopathogenic species. Given this expanding genomic and ecological knowledge, we believe the time is ripe to evaluate what we know about the evolutionary dynamics of plant pathogens.

The O-antigen mediates differential survival of Salmonella against communities of natural predators

Antigenically distinct members of bacterial species can be differentially distributed in the environment. Predators known to consume antigenically distinct prey with different efficiencies are also differentially distributed. Here we show that antigenically distinct, but otherwise isogenic and physiologically indistinct, strains of Salmonella enterica show differential survival in natural soil, sediment and intestinal environments, where they would face a community of predators. Decline in overall cell numbers is attenuated by factors that inhibit the action of predators, including heat and antiprotozoal and antihelminthic drugs. Moreover, the fitness of strains facing these predators - calculated by comparing survival with and without treatments attenuating predator activity - varies between environments. These results suggest that relative survival in natural environments is arbitrated by communities of natural predators whose feeding preferences, if not species composition, vary between environments. These data support the hypothesis that survival against natural predators may drive the differential distribution of bacteria among microenvironments.

Gene flow in microbial communities could explain unexpected patterns of synonymous variation in the Escherichia coli core genome

Researchers contest the importance of gene flow in bacterial core genomes, as traditionalists view microbes as predominantly clonal, asexually reproducing organisms. Contrary to the traditional perspective, Escherichia coli core genes vary greatly in their levels of synonymous genetic diversity. This observation indicates that the relative importance of evolutionary forces such as mutation, selection, and recombination varies from gene to gene. In this paper, I highlight why the synonymous diversity observation is broadly relevant to researchers interested in the evolutionary dynamics of microbial populations and communities. I explain how a model of evolution called the coalescent relates neutral diversity (i.e. mutations with negligible fitness effects) to mutation rates, evolutionary time, and a parameter called effective population size. I then describe the possible ways in which mutation, selection, and recombination can explain observed patterns of synonymous diversity in E. coli. Finally, I describe a model for E. coli genome evolution in which different loci are subject to varying levels of gene flow among co-occurring microbes and viruses in the environment. Researchers can falsify the gene flow hypothesis by sequencing genes and strains isolated from stable microbiomes or by carrying out evolution experiments that trace gene genealogies in real-time.

The Contribution of Genetic Recombination to CRISPR Array Evolution

CRISPR (clustered regularly interspaced short palindromic repeats) is a microbial immune system against foreign DNA. Recognition sequences (spacers) encoded within the CRISPR array mediate the immune reaction in a sequence-specific manner. The known mechanisms for the evolution of CRISPR arrays include spacer acquisition from foreign DNA elements at the time of invasion and array erosion through spacer deletion. Here, we consider the contribution of genetic recombination between homologous CRISPR arrays to the evolution of spacer repertoire. Acquisition of spacers from exogenic arrays via recombination may confer the recipient with immunity against unencountered antagonists. For this purpose, we develop a novel method for the detection of recombination in CRISPR arrays by modeling the spacer order in arrays from multiple strains from the same species. Because the evolutionary signal of spacer recombination may be similar to that of pervasive spacer deletions or independent spacer acquisition, our method entails a robustness analysis of the recombination inference by a statistical comparison to resampled and perturbed data sets. We analyze CRISPR data sets from four bacterial species: two Gammaproteobacteria species harboring CRISPR type I and two Streptococcus species harboring CRISPR type II loci. We find that CRISPR array evolution in Escherichia coli and Streptococcus agalactiae can be explained solely by vertical inheritance and differential spacer deletion. In Pseudomonas aeruginosa, we find an excess of single spacers potentially incorporated into the CRISPR locus during independent acquisition events. In Streptococcus thermophilus, evidence for spacer acquisition by recombination is present in 5 out of 70 strains. Genetic recombination has been proposed to accelerate adaptation by combining beneficial mutations that arose in independent lineages. However, for most species under study, we find that CRISPR evolution is shaped mainly by spacer acquisition and loss rather than recombination. Since the evolution of spacer content is characterized by a rapid turnover, it is likely that recombination is not beneficial for improving phage resistance in the strains under study, or that it cannot be detected in the resolution of intraspecies comparisons.

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.

GC-Content evolution in bacterial genomes: the biased gene conversion hypothesis expands

The characterization of functional elements in genomes relies on the identification of the footprints of natural selection. In this quest, taking into account neutral evolutionary processes such as mutation and genetic drift is crucial because these forces can generate patterns that may obscure or mimic signatures of selection. In mammals, and probably in many eukaryotes, another such confounding factor called GC-Biased Gene Conversion (gBGC) has been documented. This mechanism generates patterns identical to what is expected under selection for higher GC-content, specifically in highly recombining genomic regions. Recent results have suggested that a mysterious selective force favouring higher GC-content exists in Bacteria but the possibility that it could be gBGC has been excluded. Here, we show that gBGC is probably at work in most if not all bacterial species. First we find a consistent positive relationship between the GC-content of a gene and evidence of intra-genic recombination throughout a broad spectrum of bacterial clades. Second, we show that the evolutionary force responsible for this pattern is acting independently from selection on codon usage, and could potentially interfere with selection in favor of optimal AU-ending codons. A comparison with data from human populations shows that the intensity of gBGC in Bacteria is comparable to what has been reported in mammals. We propose that gBGC is not restricted to sexual Eukaryotes but also widespread among Bacteria and could therefore be an ancestral feature of cellular organisms. We argue that if gBGC occurs in bacteria, it can account for previously unexplained observations, such as the apparent non-equilibrium of base substitution patterns and the heterogeneity of gene composition within bacterial genomes. Because gBGC produces patterns similar to positive selection, it is essential to take this process into account when studying the evolutionary forces at work in bacterial genomes.

Classical population genetics models indicate that the efficiency of selection, and hence adaptation, depends on a number of non-selective factors, such as the size of a population or the intensity of recombination. In the last 10 years, evidence has accumulated that another mechanism called GC-Biased Gene Conversion (gBGC) can interfere with selection and even mimic its effects. This phenomenon, which arises from a particularity of the recombination machinery, was first thought to be restricted to sexual eukaryotic organisms. Here, we show that this mechanism probably exists in Bacteria and has a strong impact on their genome evolution. This discovery not only explains many previously unconnected features of bacterial genome evolution, but also highlights the importance of non-adaptive evolutionary processes in Bacteria.

Origins of major archaeal clades correspond to gene acquisitions from bacteria

The mechanisms that underlie the origin of major prokaryotic groups are poorly understood. In principle, the origin of both species and higher taxa among prokaryotes should entail similar mechanisms — ecological interactions with the environment paired with natural genetic variation involving lineage-specific gene innovations and lineage-specific gene acquisitions^1,2,3,4. To investigate the origin of higher taxa in archaea, we have determined gene distributions and gene phylogenies for the 267,568 protein coding genes of 134 sequenced archaeal genomes in the context of their homologs from 1,847 reference bacterial genomes. Archaea-specific gene families define 13 traditionally recognized archaeal higher taxa in our sample. Here we report that the origins of these 13 groups unexpectedly correspond to 2,264 group-specific gene acquisitions from bacteria. Interdomain gene transfer is highly asymmetric, transfers from bacteria to archaea are more than 5-fold more frequent than vice versa. Gene transfers identified at major evolutionary transitions among prokaryotes specifically implicate gene acquisitions for metabolic functions from bacteria as key innovations in the origin of higher archaeal taxa.

Divergence with gene flow in a population of thermophilic bacteria: a potential role for spatially varying selection

A fundamental goal of evolutionary biology is to understand how ecological diversity arises and is maintained in natural populations. We have investigated the contributions of gene flow and divergent selection to the distribution of genetic variation in an ecologically differentiated population of a thermophilic cyanobacterium (Mastigocladus laminosus) found along the temperature gradient of a nitrogen-limited stream in Yellowstone National Park. For most loci sampled, gene flow appears to be sufficient to prevent substantial genetic divergence. However, one locus (rfbC) exhibited a comparatively low migration rate as well as other signatures expected for a gene experiencing spatially varying selection, including an excess of common variants, an elevated level of polymorphism and extreme genetic differentiation along the gradient. rfbC is part of an expression island involved in the production of the polysaccharide component of the protective envelope of the heterocyst, the specialized nitrogen-fixing cell of these bacteria. SNP genotyping in the vicinity of rfbC revealed a ~5-kbp region including a gene content polymorphism that is tightly associated with environmental temperature and therefore likely contains the target of selection. Two genes have been deleted both in the predominant haplotype found in the downstream region of White Creek and in strains from other Yellowstone populations of M. laminosus, which may result in the production of heterocysts with different envelope properties. This study implicates spatially varying selection in the maintenance of variation related to thermal performance at White Creek despite on-going or recent gene flow.

MALDI-TOF mass spectrometry confirms difficulties in separating species of the Avibacterium genus

In the present study a well-characterized strain collection (n = 33) of Avibacterium species was investigated by matrix-assisted laser desorption ionization-time-of flight mass spectrometry (MALDI-TOF MS). The robustness of the currently available reference database (Bruker Biotyper 3.0) was tested to determine the degree of identification of these strains. Reproducible signal patterns were obtained from all strains. However, identification of most strains was only possible at genus level. Furthermore, two strains could not be identified by this method. Based on their protein spectra profiles, a MALDI main spectra dendrogram was created to determine their relationship. Most strains were closely related-for example, 26 strains formed cluster 1 including the type strains of Avibacterium volantium, Avibacterium gallinarum, Avibacterium endocarditidis and Avibacterium avium-while Avibacterium paragallinarum biovars 1 and 2 formed cluster 2 and, finally, strain 55000 remained on its own. The present MALDI-TOF MS results confirm recent findings that only certain isolates of Av. paragallinarum represent a well-defined species within the genus Avibacterium, making a taxonomic revision essential. To improve identification of Avibacterium at species level by MALDI-TOF MS, relevant reference strains were included in the newly created database and results are presented. In conclusion, Av. paragallinarum can be identified by MALDI/Biotyper and not the other species of the genus.

Detecting rare gene transfer events in bacterial populations

Horizontal gene transfer (HGT) enables bacteria to access, share, and recombine genetic variation, resulting in genetic diversity that cannot be obtained through mutational processes alone. In most cases, the observation of evolutionary successful HGT events relies on the outcome of initially rare events that lead to novel functions in the new host, and that exhibit a positive effect on host fitness. Conversely, the large majority of HGT events occurring in bacterial populations will go undetected due to lack of replication success of transformants. Moreover, other HGT events that would be highly beneficial to new hosts can fail to ensue due to lack of physical proximity to the donor organism, lack of a suitable gene transfer mechanism, genetic compatibility, and stochasticity in tempo-spatial occurrence. Experimental attempts to detect HGT events in bacterial populations have typically focused on the transformed cells or their immediate offspring. However, rare HGT events occurring in large and structured populations are unlikely to reach relative population sizes that will allow their immediate identification; the exception being the unusually strong positive selection conferred by antibiotics. Most HGT events are not expected to alter the likelihood of host survival to such an extreme extent, and will confer only minor changes in host fitness. Due to the large population sizes of bacteria and the time scales involved, the process and outcome of HGT are often not amenable to experimental investigation. Population genetic modeling of the growth dynamics of bacteria with differing HGT rates and resulting fitness changes is therefore necessary to guide sampling design and predict realistic time frames for detection of HGT, as it occurs in laboratory or natural settings. Here we review the key population genetic parameters, consider their complexity and highlight knowledge gaps for further research.

Genetic characterization of atypical Citrobacter freundii

The ability of a bacterial population to survive in different niches, as well as in stressful and rapidly changing environmental conditions, depends greatly on its genetic content. To survive such fluctuating conditions, bacteria have evolved different mechanisms to modulate phenotypic variations and related strategies to produce high levels of genetic diversity. Laboratories working in microbiological diagnosis have shown that Citrobacter freundii is very versatile in its colony morphology, as well as in its biochemical, antigenic and pathogenic behaviours. This phenotypic versatility has made C. freundii difficult to identify and it is frequently confused with both Salmonella enterica and Escherichia coli. In order to determine the genomic events and to explain the mechanisms involved in this plasticity, six C. freundii isolates were selected from a phenotypic variation study. An I-CeuI genomic cleavage map was created and eight housekeeping genes, including 16S rRNA, were sequenced. In general, the results showed a range of both phenotypes and genotypes among the isolates with some revealing a greater similarity to C. freundii and some to S. enterica, while others were identified as phenotypic and genotypic intermediary states between the two species. The occurrence of these events in natural populations may have important implications for genomic diversification in bacterial evolution, especially when considering bacterial species boundaries. In addition, such events may have a profound impact on medical science in terms of treatment, course and outcomes of infectious diseases, evading the immune response, and understanding host-pathogen interactions.

Analysis of genome content evolution in pvc bacterial super-phylum: assessment of candidate genes associated with cellular organization and lifestyle

The Planctomycetes, Verrucomicrobia, Chlamydiae (PVC) super-phylum contains bacteria with either complex cellular organization or simple cell structure; it also includes organisms of different lifestyles (pathogens, mutualists, commensal, and free-living). Genome content evolution of this group has not been studied in a systematic fashion, which would reveal genes underlying the emergence of PVC-specific phenotypes. Here, we analyzed the evolutionary dynamics of 26 PVC genomes and several outgroup species. We inferred HGT, duplications, and losses by reconciliation of 27,123 gene trees with the species phylogeny. We showed that genome expansion and contraction have driven evolution within Planctomycetes and Chlamydiae, respectively, and balanced each other in Verrucomicrobia and Lentisphaerae. We also found that for a large number of genes in PVC genomes the most similar sequences are present in Acidobacteria, suggesting past and/or current ecological interaction between organisms from these groups. We also found evidence of shared ancestry between carbohydrate degradation genes in the mucin-degrading human intestinal commensal Akkermansia muciniphila and sequences from Acidobacteria and Bacteroidetes, suggesting that glycoside hydrolases are transferred laterally between gut microbes and that the process of carbohydrate degradation is crucial for microbial survival within the human digestive system. Further, we identified a highly conserved genetic module preferentially present in compartmentalized PVC species and possibly associated with the complex cell plan in these organisms. This conserved machinery is likely to be membrane targeted and involved in electron transport, although its exact function is unknown. These genes represent good candidates for future functional studies.

Comparative genomic analysis of the genus Nocardiopsis provides new insights into its genetic mechanisms of environmental adaptability

The genus Nocardiopsis, a widespread group in phylum Actinobacteria, has received much attention owing to its ecological versatility, pathogenicity, and ability to produce a rich array of bioactive metabolites. Its high environmental adaptability might be attributable to its genome dynamics, which can be estimated through comparative genomic analysis targeting microorganisms with close phylogenetic relationships but different phenotypes. To shed light on speciation, gene content evolution, and environmental adaptation in these unique actinobacteria, we sequenced draft genomes for 16 representative species of the genus and compared them with that of the type species N. dassonvillei subsp. dassonvillei DSM 43111^T. The core genome of 1,993 orthologous and paralogous gene clusters was identified, and the pan-genomic reservoir was found not only to accommodate more than 22,000 genes, but also to be open. The top ten paralogous genes in terms of copy number could be referred to three functional categories: transcription regulators, transporters, and synthases related to bioactive metabolites. Based on phylogenomic reconstruction, we inferred past evolutionary events, such as gene gains and losses, and identified a list of clade-specific genes implicated in environmental adaptation. These results provided insights into the genetic causes of environmental adaptability in this cosmopolitan actinobacterial group and the contributions made by its inherent features, including genome dynamics and the constituents of core and accessory proteins.

Comparative genomics of Wolbachia and the bacterial species concept

The importance of host-specialization to speciation processes in obligate host-associated bacteria is well known, as is also the ability of recombination to generate cohesion in bacterial populations. However, whether divergent strains of highly recombining intracellular bacteria, such as Wolbachia, can maintain their genetic distinctness when infecting the same host is not known. We first developed a protocol for the genome sequencing of uncultivable endosymbionts. Using this method, we have sequenced the complete genomes of the Wolbachia strains wHa and wNo, which occur as natural double infections in Drosophila simulans populations on the Seychelles and in New Caledonia. Taxonomically, wHa belong to supergroup A and wNo to supergroup B. A comparative genomics study including additional strains supported the supergroup classification scheme and revealed 24 and 33 group-specific genes, putatively involved in host-adaptation processes. Recombination frequencies were high for strains of the same supergroup despite different host-preference patterns, leading to genomic cohesion. The inferred recombination fragments for strains of different supergroups were of short sizes, and the genomes of the co-infecting Wolbachia strains wHa and wNo were not more similar to each other and did not share more genes than other A- and B-group strains that infect different hosts. We conclude that Wolbachia strains of supergroup A and B represent genetically distinct clades, and that strains of different supergroups can co-exist in the same arthropod host without converging into the same species. This suggests that the supergroups are irreversibly separated and that barriers other than host-specialization are able to maintain distinct clades in recombining endosymbiont populations. Acquiring a good knowledge of the barriers to genetic exchange in Wolbachia will advance our understanding of how endosymbiont communities are constructed from vertically and horizontally transmitted genes.

Speciation in sexual organisms is defined as the inability of two populations to get viable offspring. Speciation in asexual, obligate endosymbionts is thought to be an indirect consequence of host-specialization. An important question is if divergent endosymbionts would start blending if the host barrier isolating them were removed. Here, we have studied Wolbachia, an abundant group of bacteria in the insect world. Wolbachia is classified into supergroups based on multi-locus sequence typing. We have sequenced the genomes from the Wolbachia strains wNo and wHa. These are particularly interesting since they belong to different supergroups yet co-occur as a double-infection in natural populations of Drosophila simulans. A comparative genomics study showed that wHa and wNo contain no uniquely shared genes. Instead, each strain shares unique gene functions with members of the same supergroup that infect other hosts. This unexpected finding suggests an alternative means of ecological speciation, indicating that speciation is not restricted to host-specialization but rather that related endosymbionts can coexist as separate species in the same host. Our study sheds light on the genomic divergence between different partners inhabiting the intracellular niche of the same host organism.

Evolutionary Genomics of Salmonella enterica Subspecies

Six subspecies are currently recognized in Salmonella enterica. Subspecies I (subspecies enterica) is responsible for nearly all infections in humans and warm-blooded animals, while five other subspecies are isolated principally from cold-blooded animals. We sequenced 21 phylogenetically diverse strains, including two representatives from each of the previously unsequenced five subspecies and 11 diverse new strains from S. enterica subspecies enterica, to put this species into an evolutionary perspective. The phylogeny of the subspecies was partly obscured by abundant recombination events between lineages and a relatively short period of time within which subspeciation took place. Nevertheless, a variety of different tree-building methods gave congruent evolutionary tree topologies for subspeciation. A total of 285 gene families were identified that were recruited into subspecies enterica, and most of these are of unknown function. At least 2,807 gene families were identified in one or more of the other subspecies that are not found in subspecies I or Salmonella bongori. Among these gene families were 13 new candidate effectors and 7 new candidate fimbrial clusters. A third complete type III secretion system not present in subspecies enterica (I) isolates was found in both strains of subspecies salamae (II). Some gene families had complex taxonomies, such as the type VI secretion systems, which were recruited from four different lineages in five of six subspecies. Analysis of nonsynonymous-to-synonymous substitution rates indicated that the more-recently acquired regions in S. enterica are undergoing faster fixation rates than the rest of the genome. Recently acquired AT-rich regions, which often encode virulence functions, are under ongoing selection to maintain their high AT content.

We have sequenced 21 new genomes which encompass the phylogenetic diversity of Salmonella, including strains of the previously unsequenced subspecies arizonae, diarizonae, houtenae, salamae, and indica as well as new diverse strains of subspecies enterica. We have deduced possible evolutionary paths traversed by this very important zoonotic pathogen and identified novel putative virulence factors that are not found in subspecies I. Gene families gained at the time of the evolution of subspecies enterica are of particular interest because they include mechanisms by which this subspecies adapted to warm-blooded hosts.

Multilocus sequence analysis of Pasteurella multocida demonstrates a type species under development

The aim of the present study was to use multilocus sequence typing (MLST) of a diverse collection of Pasteurella multocida with regard to animal source, place and date of collection, including all available serovars of Carter, Heddleston, Little & Lyon, Namioka, Cornelius and Roberts, to further investigate the evolution of this species with a focus on two lineages, A (P. multocida subsp. multocida and P. multocida subsp. gallicida) and B (P. multocida subsp. septica), previously reported. Isolates of P. multocida (n = 116) including reference strains of major serotyping systems were investigated by MLST based on partial sequences of the genes adk, est, gdh, mdh, pgi, pmi and zwf, and 67 sequence types (STs) were observed. Phylogenetic analysis of these concatenated sequences confirmed the separation of groups A (41 STs, 71 isolates) and B (22 STs, 38 isolates) out of the 67 STs. All Carter serovars, 12 Heddleston serovars, all three Little-Lyon types, six out of seven Namioka serovars, all five Roberts types and all four Cornelius serovars were allocated to the A group, while group B included the remaining four Heddleston serovars, 6, 7, 8 and 13, in addition to Namioka type 8 : A. The overrepresentation of reference strains of serotyping systems in the A group contrasts with the high number of isolates obtained from diseased birds in the B group, the effect of which should be addressed in future vaccine development. Isolates from birds (25) dominated the B group, which also included four isolates from Felidae, whereas group A included isolates from all types of hosts. The evolutionary implications of the lack of capsular type D, pig and bovine isolates in group B, as well as its association with Aves and Felidae that also applied to the whole Rural Industries Research and Development Corporation (RIRDC) MLST database, need further investigation. The combination of rpoB and 16S rRNA gene sequence comparison as well as the developed PCR test assigned isolates to lineage A, represented by the type strain of P. multocida subsp. Multocida, or lineage B represented by the type strain of P. multocida subsp. septica. It was not possible to circumscribe either the A or B lineages with a set of conserved phenotypic characters, calling into question the validity of subspecies within P. multocida. Phylogenetic analysis carried out on individual MLST genes showed deviations as to single or multiple genes for 17 % of group A and 43 % of group B, indicating that lineage A probably developed from lineage B, and that major changes are ongoing. From a genotypical point of view, we conclude that P. multocida subsp. gallicida represents an artificial unit.

Reproducibility of Vibrionaceae population structure in coastal bacterioplankton

How reproducibly microbial populations assemble in the wild remains poorly understood. Here, we assess evidence for ecological specialization and predictability of fine-scale population structure and habitat association in coastal ocean Vibrionaceae across years. We compare Vibrionaceae lifestyles in the bacterioplankton (combinations of free-living, particle, or zooplankton associations) measured using the same sampling scheme in 2006 and 2009 to assess whether the same groups show the same environmental association year after year. This reveals complex dynamics with populations falling primarily into two categories: (i) nearly equally represented in each of the two samplings and (ii) highly skewed, often to an extent that they appear exclusive to one or the other sampling times. Importantly, populations recovered at the same abundance in both samplings occupied highly similar habitats suggesting predictable and robust environmental association while skewed abundances of some populations may be triggered by shifts in ecological conditions. The latter is supported by difference in the composition of large eukaryotic plankton between years, with samples in 2006 being dominated by copepods, and those in 2009 by diatoms. Overall, the comparison supports highly predictable population-habitat linkage but highlights the fact that complex, and often unmeasured, environmental dynamics in habitat occurrence may have strong effects on population dynamics.

Natural transformation facilitates transfer of transposons, integrons and gene cassettes between bacterial species

We have investigated to what extent natural transformation acting on free DNA substrates can facilitate transfer of mobile elements including transposons, integrons and/or gene cassettes between bacterial species. Naturally transformable cells of Acinetobacter baylyi were exposed to DNA from integron-carrying strains of the genera Acinetobacter, Citrobacter, Enterobacter, Escherichia, Pseudomonas, and Salmonella to determine the nature and frequency of transfer. Exposure to the various DNA sources resulted in acquisition of antibiotic resistance traits as well as entire integrons and transposons, over a 24 h exposure period. DNA incorporation was not solely dependent on integrase functions or the genetic relatedness between species. DNA sequence analyses revealed that several mechanisms facilitated stable integration in the recipient genome depending on the nature of the donor DNA; homologous or heterologous recombination and various types of transposition (Tn21-like and IS26-like). Both donor strains and transformed isolates were extensively characterized by antimicrobial susceptibility testing, integron- and cassette-specific PCRs, DNA sequencing, pulsed field gel electrophoreses (PFGE), Southern blot hybridizations, and by re-transformation assays. Two transformant strains were also genome-sequenced. Our data demonstrate that natural transformation facilitates interspecies transfer of genetic elements, suggesting that the transient presence of DNA in the cytoplasm may be sufficient for genomic integration to occur. Our study provides a plausible explanation for why sequence-conserved transposons, IS elements and integrons can be found disseminated among bacterial species. Moreover, natural transformation of integron harboring populations of competent bacteria revealed that interspecies exchange of gene cassettes can be highly efficient, and independent on genetic relatedness between donor and recipient. In conclusion, natural transformation provides a much broader capacity for horizontal acquisitions of genetic elements and hence, resistance traits from divergent species than previously assumed.

Genetic elements, such as transposons and integrons, frequently carry antimicrobial resistance determinants and can be found widely disseminated among pathogenic bacteria. Their distribution pattern suggests dissemination through horizontal gene transfer. The role of natural transformation in horizontal transfer of genetic elements other than those that are self-replicative (plasmids) has remained largely unexplored. We have tested if natural transformation can facilitate transfer of transposons and class 1 integrons between bacterial species. We here provide experimental evidence showing that natural transformation can be a general mechanism for dissemination of genetic elements that by themselves do not encode interspecies transfer functions (e.g. transposons, insertion sequences). We demonstrate that antibiotic resistance determinants present in such genetic elements can spread by natural transformation between species of clinical interest. We show by quantitative data that interspecies exchange of resistance gene cassettes is highly efficient among integron-containing strains and species. Our study also provides a plausible explanation for how sequence-conserved integrons can become distributed among bacterial species.

The ecology of bacterial genes and the survival of the new

Much of the observed variation among closely related bacterial genomes is attributable to gains and losses of genes that are acquired horizontally as well as to gene duplications and larger amplifications. The genomic flexibility that results from these mechanisms certainly contributes to the ability of bacteria to survive and adapt in varying environmental challenges. However, the duplicability and transferability of individual genes imply that natural selection should operate, not only at the organismal level, but also at the level of the gene. Genes can be considered semiautonomous entities that possess specific functional niches and evolutionary dynamics. The evolution of bacterial genes should respond both to selective pressures that favor competition, mostly among orthologs or paralogs that may occupy the same functional niches, and cooperation, with the majority of other genes coexisting in a given genome. The relative importance of either type of selection is likely to vary among different types of genes, based on the functional niches they cover and on the tightness of their association with specific organismal lineages. The frequent availability of new functional niches caused by environmental changes and biotic evolution should enable the constant diversification of gene families and the survival of new lineages of genes.

Ecological adaptation in bacteria: speciation driven by codon selection

In bacteria, physiological change may be effected by a single gene acquisition, producing ecological differentiation without genetic isolation. Natural selection acting on such differences can reduce the frequency of genotypes that arise from recombination at these loci. However, gene acquisition can only account for recombination interference in the fraction of the genome that is tightly linked to the integration site. To identify additional loci that contribute to adaptive differences, we examined orthologous genes in species of Enterobacteriaceae to identify significant differences in the degree of codon selection. Significance was assessed using the Adaptive Codon Enrichment metric, which accounts for the variation in codon usage bias that is expected to arise from mutation and drift; large differences in codon usage bias were identified in more genes than would be expected to arise from stochastic processes alone. Genes in the same operon showed parallel differences in codon usage bias, suggesting that changes in the overall levels of gene expression led to changes in the degree of adaptive codon usage. Most significant differences between orthologous operons were found among those involved with specific environmental adaptations, whereas "housekeeping" genes rarely showed significant changes. When considered together, the loci experiencing significant changes in codon selection outnumber potentially adaptive gene acquisition events. The identity of genes under strong codon selection seems to be influenced by the habitat from which the bacteria were isolated. We propose a two-stage model for how adaptation to different selective regimes can drive bacterial speciation. Initially, gene acquisitions catalyze rapid ecological differentiation, which modifies the utilization of genes, thereby changing the strength of codon selection on them. Alleles develop fitness variation by substitution, producing recombination interference at these loci in addition to those flanking acquired genes, allowing sequences to diverge across the entire genome and establishing genetic isolation (i.e., protection from frequent homologous recombination).