Recombination and the nature of bacterial speciation

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Recombination and horizontal transfer of nodulation and ACC deaminase (acdS) genes within Alpha- and Betaproteobacteria nodulating legumes of the Cape Fynbos biome

The goal of this work is to study the evolution and the degree of horizontal gene transfer (HGT) within rhizobial genera of both Alphaproteobacteria (Mesorhizobium, Rhizobium) and Betaproteobacteria (Burkholderia), originating from South African Fynbos legumes. By using a phylogenetic approach and comparing multiple chromosomal and symbiosis genes, we revealed conclusive evidence of high degrees of horizontal transfer of nodulation genes among closely related species of both groups of rhizobia, but also among species with distant genetic backgrounds (Rhizobium and Mesorhizobium), underscoring the importance of lateral transfer of symbiosis traits as an important evolutionary force among rhizobia of the Cape Fynbos biome. The extensive exchange of symbiosis genes in the Fynbos is in contrast with a lack of significant events of HGT among Burkholderia symbionts from the South American Cerrado and Caatinga biome. Furthermore, homologous recombination among selected housekeeping genes had a substantial impact on sequence evolution within Burkholderia and Mesorhizobium. Finally, phylogenetic analyses of the non-symbiosis acdS gene in Mesorhizobium, a gene often located on symbiosis islands, revealed distinct relationships compared to the chromosomal and symbiosis genes, suggesting a different evolutionary history and independent events of gene transfer. The observed events of HGT and incongruence between different genes necessitate caution in interpreting topologies from individual data types.

Positive selection is the main driving force for evolution of citrus canker-causing Xanthomonas

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA(w)) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA(w) and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in 'carbohydrate transport and metabolism' and 'DNA replication, recombination and repair' genes (P<0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.

Predicting microbial traits with phylogenies

Phylogeny reflects genetic and phenotypic traits in Bacteria and Archaea. The phylogenetic conservatism of microbial traits has prompted the application of phylogeny-based algorithms to predict unknown trait values of extant taxa based on the traits of their evolutionary relatives to estimate, for instance, rRNA gene copy numbers, gene contents or tolerance to abiotic conditions. Unlike the 'macrobial' world, microbial ecologists face scenarios potentially compromising the accuracy of trait reconstruction methods, as, for example, extremely large phylogenies and limited information on the traits of interest. We review 990 bacterial and archaeal traits from the literature and support that phylogenetic trait conservatism is widespread through the tree of life, while revealing that it is generally weak for ecologically relevant phenotypic traits and high for genetically complex traits. We then perform a simulation exercise to assess the accuracy of phylogeny-based trait predictions in common scenarios faced by microbial ecologists. Our simulations show that ca. 60% of the variation in phylogeny-based trait predictions depends on the magnitude of the trait conservatism, the number of species in the tree, the proportion of species with unknown trait values and the mean distance in the tree to the nearest neighbour with a known trait value. Results are similar for both binary and continuous traits. We discuss these results under the light of the reviewed traits and provide recommendations for the use of phylogeny-based trait predictions for microbial ecologists.

Microbial Speciation

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

Bacterial genospecies that are not ecologically coherent: population genomics of Rhizobium leguminosarum

Biological species may remain distinct because of genetic isolation or ecological adaptation, but these two aspects do not always coincide. To establish the nature of the species boundary within a local bacterial population, we characterized a sympatric population of the bacterium Rhizobium leguminosarum by genomic sequencing of 72 isolates. Although all strains have 16S rRNA typical of R. leguminosarum, they fall into five genospecies by the criterion of average nucleotide identity (ANI). Many genes, on plasmids as well as the chromosome, support this division: recombination of core genes has been largely within genospecies. Nevertheless, variation in ecological properties, including symbiotic host range and carbon-source utilization, cuts across these genospecies, so that none of these phenotypes is diagnostic of genospecies. This phenotypic variation is conferred by mobile genes. The genospecies meet the Mayr criteria for biological species in respect of their core genes, but do not correspond to coherent ecological groups, so periodic selection may not be effective in purging variation within them. The population structure is incompatible with traditional ‘polyphasic taxonomy′ that requires bacterial species to have both phylogenetic coherence and distinctive phenotypes. More generally, genomics has revealed that many bacterial species share adaptive modules by horizontal gene transfer, and we envisage a more consistent taxonomic framework that explicitly recognizes this. Significant phenotypes should be recognized as ‘biovars' within species that are defined by core gene phylogeny.

Inferring Speciation Processes from Patterns of Natural Variation in Microbial Genomes

Microbial species concepts have long been the focus of contentious debate, fueled by technological limitations to the genetic resolution of species, by the daunting task of investigating phenotypic variation among individual microscopic organisms, and by a lack of understanding of gene flow in reproductively asexual organisms that are prone to promiscuous horizontal gene transfer. Population genomics, the emerging approach of analyzing the complete genomes of a multitude of closely related organisms, is poised to overcome these limitations by providing a window into patterns of genome variation revealing the evolutionary processes through which species diverge. This new approach is more than just an extension of previous multilocus sequencing technologies, in that it provides a comprehensive view of interacting evolutionary processes. Here we argue that the application of population genomic tools in a rigorous population genetic framework will help to identify the processes of microbial speciation and ultimately lead to a general species concept based on the unique biology and ecology of microorganisms.

Competitive strategies differentiate closely related species of marine actinobacteria

Although competition, niche partitioning, and spatial isolation have been used to describe the ecology and evolution of macro-organisms, it is less clear to what extent these principles account for the extraordinary levels of bacterial diversity observed in nature. Ecological interactions among bacteria are particularly challenging to address due to methodological limitations and uncertainties over how to recognize fundamental units of diversity and link them to the functional traits and evolutionary processes that led to their divergence. Here we show that two closely related marine actinomycete species can be differentiated based on competitive strategies. Using a direct challenge assay to investigate inhibitory interactions with members of the bacterial community, we observed a temporal difference in the onset of inhibition. The majority of inhibitory activity exhibited by Salinispora arenicola occurred early in its growth cycle and was linked to antibiotic production. In contrast, most inhibition by Salinispora tropica occurred later in the growth cycle and was more commonly linked to nutrient depletion or other sources. Comparative genomics support these differences, with S. arenicola containing nearly twice the number of secondary metabolite biosynthetic gene clusters as S. tropica, indicating a greater potential for secondary metabolite production. In contrast, S. tropica is enriched in gene clusters associated with the acquisition of growth-limiting nutrients such as iron. Coupled with differences in growth rates, the results reveal that S. arenicola uses interference competition at the expense of growth, whereas S. tropica preferentially employs a strategy of exploitation competition. The results support the ecological divergence of two co-occurring and closely related species of marine bacteria by providing evidence they have evolved fundamentally different strategies to compete in marine sediments.

Recombinant transfer in the basic genome of Escherichia coli

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

Microbial diversity. Fine-scale diversity and extensive recombination in a quasisexual bacterial population occupying a broad niche

Extensive fine-scale genetic diversity is found in many microbial species across varied environments, but for most, the evolutionary scenarios that generate the observed variation remain unclear. Deep sequencing of a thermophilic cyanobacterial population and analysis of the statistics of synonymous single-nucleotide polymorphisms revealed a high rate of homologous recombination and departures from neutral drift consistent with the effects of genetic hitchhiking. A sequenced isolate genome resembled an unlinked random mixture of the allelic diversity at the sampled loci. These observations suggested a quasisexual microbial population that occupies a broad ecological niche, with selection driving frequencies of alleles rather than whole genomes.

Soil microbial community structure is unaltered by plant invasion, vegetation clipping, and nitrogen fertilization in experimental semi-arid grasslands

Global and regional environmental changes often co-occur, creating complex gradients of disturbance on the landscape. Soil microbial communities are an important component of ecosystem response to environmental change, yet little is known about how microbial structure and function respond to multiple disturbances, or whether multiple environmental changes lead to unanticipated interactive effects. Our study used experimental semi-arid grassland plots in a Mediterranean-climate to determine how soil microbial communities in a seasonally variable ecosystem respond to one, two, or three simultaneous environmental changes: exotic plant invasion, plant invasion + vegetation clipping (to simulate common management practices like mowing or livestock grazing), plant invasion + nitrogen (N) fertilization, and plant invasion + clipping + N fertilization. We examined microbial community structure 5-6 years after plot establishment via sequencing of >1 million 16S rRNA genes. Abiotic soil properties (soil moisture, temperature, pH, and inorganic N) and microbial functioning (nitrification and denitrification potentials) were also measured and showed treatment-induced shifts, including altered NO(-) 3 availability, temperature, and nitrification potential. Despite these changes, bacterial and archaeal communities showed little variation in composition and diversity across treatments. Even communities in plots exposed to three interacting environmental changes were similar to those in restored native grassland plots. Historical exposure to large seasonal and inter-annual variations in key soil properties, in addition to prior site cultivation, may select for a functionally plastic or largely dormant microbial community, resulting in a microbial community that is structurally robust to single and multiple environmental changes.

Involvement of aph(3')-IIa in the formation of mosaic aminoglycoside resistance genes in natural environments

Intragenic recombination leading to mosaic gene formation is known to alter resistance profiles for particular genes and bacterial species. Few studies have examined to what extent aminoglycoside resistance genes undergo intragenic recombination. We screened the GenBank database for mosaic gene formation in homologs of the aph(3')-IIa (nptII) gene. APH(3')-IIa inactivates important aminoglycoside antibiotics. The gene is widely used as a selectable marker in biotechnology and enters the environment via laboratory discharges and the release of transgenic organisms. Such releases may provide opportunities for recombination in competent environmental bacteria. The retrieved GenBank sequences were grouped in three datasets comprising river water samples, duck pathogens and full-length variants from various bacterial genomes and plasmids. Analysis for recombination in these datasets was performed with the Recombination Detection Program (RDP4), and the Genetic Algorithm for Recombination Detection (GARD). From a total of 89 homologous sequences, 83% showed 99-100% sequence identity with aph(3')-IIa originally described as part of transposon Tn5. Fifty one were unique sequence variants eligible for recombination analysis. Only a single recombination event was identified with high confidence and indicated the involvement of aph(3')-IIa in the formation of a mosaic gene located on a plasmid of environmental origin in the multi-resistant isolate Pseudomonas aeruginosa PA96. The available data suggest that aph(3')-IIa is not an archetypical mosaic gene as the divergence between the described sequence variants and the number of detectable recombination events is low. This is in contrast to the numerous mosaic alleles reported for certain penicillin or tetracycline resistance determinants.

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.

Do Species Exist in Asexuals? Theory and Evidence from Bdelloid Rotifers

The possibility for independently evolving entities to form and persist in the absence of sexual recombination in eukaryotes has been questioned; nevertheless, there are organisms that are known to be asexual and that have apparently diversified into multiple species as recognized by taxonomists. These organisms have therefore been identified as an evolutionary paradox. We explore three alternative hypotheses attempting to solve the apparent paradox, focusing on bdelloid rotifers, the most studied group of organisms in which all species are considered asexual: (1) they may have some hidden form of sex; (2) species do not represent biological entities but simply convenient names; and (3) sex may not be a necessary requirement for speciation. We provide ample evidence against the first two hypotheses, reporting several studies supporting (1) bdelloids asexuality from different approaches, and (2) the existence of species from genetics, jaw morphology, ecology, and physiology. Thus, we (3) explore the role of sex in speciation comparing bdelloid and monogonont rotifers, and conclude with some caveats that could still change our understanding of bdelloid species.

Genetic diversity of the Pneumococcal CbpA: Implications for next-generation vaccine development

Pneumococci are capable of vaccine escape by genetic recombination at the targeted capsular locus, significantly reducing long-term vaccine effectiveness. Recently, efforts have been redirected to understanding pneumococcal biology related to potential next-generation vaccine candidates. A variety of serotype-independent protein antigens capable of inducing protective immune responses in tissue culture and animal models of infection have been identified. However, ideal vaccine candidates that are conserved across all genotypes, provide broad population coverage, and induce T-cell dependent immune responses are still under investigation. We examined whether immune responses due to the highly polymorphic CbpA antigen are due to a conserved domain capable of evoking specific immune "memory" across all genotypes of pneumococci. We defined the genotypes in a global dataset of 213 pneumococcal isolates. This isolate collection was genotypically diverse and ideal for establishing the presence of conserved CbpA epitopes as potential protein vaccine candidates. Examination of the CbpA locus sequence was highly polymorphic at both the nucleic acid and amino acid level. Despite this high polymorphism some domains are broadly conserved and consist of different amino acid residues with the same physicochemical properties, and therefore have similar tertiary structures. The two most common domains identified in the CbpA gene are modular teichoic acid phosphorylcholine esterase Pce (2bib:A), and R2 domain (1w9r:A). These conserved domains are immunogenic, therefore capable of inducing long-term host immune responses; moreover they are extracellularly located and thus accessible. We proposed their evaluation as suitable next-generation CbpA-fusion protein vaccine candidates.

Gene-specific selective sweeps in bacteria and archaea caused by negative frequency-dependent selection

BACKGROUND

Fixation of beneficial genes in bacteria and archaea (collectively, prokaryotes) is often believed to erase pre-existing genomic diversity through the hitchhiking effect, a phenomenon known as genome-wide selective sweep. Recent studies, however, indicate that beneficial genes spread through a prokaryotic population via recombination without causing genome-wide selective sweeps. These gene-specific selective sweeps seem to be at odds with the existing estimates of recombination rates in prokaryotes, which appear far too low to explain such phenomena.

RESULTS

We use mathematical modeling to investigate potential solutions to this apparent paradox. Most microbes in nature evolve in heterogeneous, dynamic communities, in which ecological interactions can substantially impact evolution. Here, we focus on the effect of negative frequency-dependent selection (NFDS) such as caused by viral predation (kill-the-winner dynamics). The NFDS maintains multiple genotypes within a population, so that a gene beneficial to every individual would have to spread via recombination, hence a gene-specific selective sweep. However, gene loci affected by NFDS often are located in variable regions of microbial genomes that contain genes involved in the mobility of selfish genetic elements, such as integrases or transposases. Thus, the NFDS-affected loci are likely to experience elevated rates of recombination compared with the other loci. Consequently, these loci might be effectively unlinked from the rest of the genome, so that NFDS would be unable to prevent genome-wide selective sweeps. To address this problem, we analyzed population genetic models of selective sweeps in prokaryotes under NFDS. The results indicate that NFDS can cause gene-specific selective sweeps despite the effect of locally elevated recombination rates, provided NFDS affects more than one locus and the basal rate of recombination is sufficiently low. Although these conditions might seem to contradict the intuition that gene-specific selective sweeps require high recombination rates, they actually decrease the effective rate of recombination at loci affected by NFDS relative to the per-locus basal level, so that NFDS can cause gene-specific selective sweeps.

CONCLUSION

Because many free-living prokaryotes are likely to evolve under NFDS caused by ubiquitous viruses, gene-specific selective sweeps driven by NFDS are expected to be a major, general phenomenon in prokaryotic populations.

Ecological Overlap and Horizontal Gene Transfer in Staphylococcus aureus and Staphylococcus epidermidis

The opportunistic pathogens Staphylococcus aureus and Staphylococcus epidermidis represent major causes of severe nosocomial infection, and are associated with high levels of mortality and morbidity worldwide. These species are both common commensals on the human skin and in the nasal pharynx, but are genetically distinct, differing at 24% average nucleotide divergence in 1,478 core genes. To better understand the genome dynamics of these ecologically similar staphylococcal species, we carried out a comparative analysis of 324 S. aureus and S. epidermidis genomes, including 83 novel S. epidermidis sequences. A reference pan-genome approach and whole genome multilocus-sequence typing revealed that around half of the genome was shared between the species. Based on a BratNextGen analysis, homologous recombination was found to have impacted on 40% of the core genes in S. epidermidis, but on only 24% of the core genes in S. aureus. Homologous recombination between the species is rare, with a maximum of nine gene alleles shared between any two S. epidermidis and S. aureus isolates. In contrast, there was considerable interspecies admixture of mobile elements, in particular genes associated with the SaPIn1 pathogenicity island, metal detoxification, and the methicillin-resistance island SCCmec. Our data and analysis provide a context for considering the nature of recombinational boundaries between S. aureus and S. epidermidis and, the selective forces that influence realized recombination between these species.

Mechanisms and impact of genetic recombination in the evolution of Streptococcus pneumoniae

Streptococcus pneumoniae (the pneumococcus) is a highly recombinogenic bacterium responsible for a high burden of human disease globally. Genetic recombination, a process in which exogenous DNA is acquired and incorporated into its genome, is a key evolutionary mechanism employed by the pneumococcus to rapidly adapt to selective pressures. The rate at which the pneumococcus acquires genetic variation through recombination is much higher than the rate at which the organism acquires variation through spontaneous mutations. This higher rate of variation allows the pneumococcus to circumvent the host innate and adaptive immune responses, escape clinical interventions, including antibiotic therapy and vaccine introduction. The rapid influx of whole genome sequence (WGS) data and the advent of novel analysis methods and powerful computational tools for population genetics and evolution studies has transformed our understanding of how genetic recombination drives pneumococcal adaptation and evolution. Here we discuss how genetic recombination has impacted upon the evolution of the pneumococcus.

Phylogenetic relatedness determined between antibiotic resistance and 16S rRNA genes in actinobacteria

BACKGROUND

Distribution and evolutionary history of resistance genes in environmental actinobacteria provide information on intensity of antibiosis and evolution of specific secondary metabolic pathways at a given site. To this day, actinobacteria producing biologically active compounds were isolated mostly from soil but only a limited range of soil environments were commonly sampled. Consequently, soil remains an unexplored environment in search for novel producers and related evolutionary questions.

RESULTS

Ninety actinobacteria strains isolated at contrasting soil sites were characterized phylogenetically by 16S rRNA gene, for presence of erm and ABC transporter resistance genes and antibiotic production. An analogous analysis was performed in silico with 246 and 31 strains from Integrated Microbial Genomes (JGI_IMG) database selected by the presence of ABC transporter genes and erm genes, respectively. In the isolates, distances of erm gene sequences were significantly correlated to phylogenetic distances based on 16S rRNA genes, while ABC transporter gene distances were not. The phylogenetic distance of isolates was significantly correlated to soil pH and organic matter content of isolation sites. In the analysis of JGI_IMG datasets the correlation between phylogeny of resistance genes and the strain phylogeny based on 16S rRNA genes or five housekeeping genes was observed for both the erm genes and ABC transporter genes in both actinobacteria and streptomycetes. However, in the analysis of sequences from genomes where both resistance genes occurred together the correlation was observed for both ABC transporter and erm genes in actinobacteria but in streptomycetes only in the erm gene.

CONCLUSIONS

The type of erm resistance gene sequences was influenced by linkage to 16S rRNA gene sequences and site characteristics. The phylogeny of ABC transporter gene was correlated to 16S rRNA genes mainly above the genus level. The results support the concept of new specific secondary metabolite scaffolds occurring more likely in taxonomically distant producers but suggest that the antibiotic selection of gene pools is also influenced by site conditions.

Population structure of Neisseria gonorrhoeae based on whole genome data and its relationship with antibiotic resistance

Neisseria gonorrhoeae is the causative agent of gonorrhea, a sexually transmitted infection (STI) of major importance. As a result of antibiotic resistance, there are now limited options for treating patients. We collected draft genome sequence data and associated metadata data on 76 N. gonorrhoeae strains from around the globe and searched for known determinants of antibiotics resistance within the strains. The population structure and evolutionary forces within the pathogen population were analyzed. Our results indicated a cosmopolitan gonoccocal population mainly made up of five subgroups. The estimated ratio of recombination to mutation (r/m = 2.2) from our data set indicates an appreciable level of recombination occurring in the population. Strains with resistance phenotypes to more recent antibiotics (azithromycin and cefixime) were mostly found in two of the five population subgroups.

DNA sequence alignment by microhomology sampling during homologous recombination

Homologous recombination (HR) mediates the exchange of genetic information between sister or homologous chromatids. During HR, members of the RecA/Rad51 family of recombinases must somehow search through vast quantities of DNA sequence to align and pair single-strand DNA (ssDNA) with a homologous double-strand DNA (dsDNA) template. Here, we use single-molecule imaging to visualize Rad51 as it aligns and pairs homologous DNA sequences in real time. We show that Rad51 uses a length-based recognition mechanism while interrogating dsDNA, enabling robust kinetic selection of 8-nucleotide (nt) tracts of microhomology, which kinetically confines the search to sites with a high probability of being a homologous target. Successful pairing with a ninth nucleotide coincides with an additional reduction in binding free energy, and subsequent strand exchange occurs in precise 3-nt steps, reflecting the base triplet organization of the presynaptic complex. These findings provide crucial new insights into the physical and evolutionary underpinnings of DNA recombination.

ClonalFrameML: efficient inference of recombination in whole bacterial genomes

Recombination is an important evolutionary force in bacteria, but it remains challenging to reconstruct the imports that occurred in the ancestry of a genomic sample. Here we present ClonalFrameML, which uses maximum likelihood inference to simultaneously detect recombination in bacterial genomes and account for it in phylogenetic reconstruction. ClonalFrameML can analyse hundreds of genomes in a matter of hours, and we demonstrate its usefulness on simulated and real datasets. We find evidence for recombination hotspots associated with mobile elements in Clostridium difficile ST6 and a previously undescribed 310kb chromosomal replacement in Staphylococcus aureus ST582. ClonalFrameML is freely available at http://clonalframeml.googlecode.com/.

Microbial taxonomy in the post-genomic era: rebuilding from scratch?

Microbial taxonomy should provide adequate descriptions of bacterial, archaeal, and eukaryotic microbial diversity in ecological, clinical, and industrial environments. Its cornerstone, the prokaryote species has been re-evaluated twice. It is time to revisit polyphasic taxonomy, its principles, and its practice, including its underlying pragmatic species concept. Ultimately, we will be able to realize an old dream of our predecessor taxonomists and build a genomic-based microbial taxonomy, using standardized and automated curation of high-quality complete genome sequences as the new gold standard.

Contrasting inter- and intraspecies recombination patterns in the "Harveyi clade" vibrio collected over large spatial and temporal scales

Recombination plays an important role in the divergence of bacteria, but the frequency of interspecies and intraspecies recombination events remains poorly understood. We investigated recombination events that occurred within core genomes of 35 Vibrio strains (family Vibrionaceae, Gammaproteobacteria), from six closely related species in the so-called “Harveyi clade.” The strains were selected from a collection of strains isolated in the last 90 years, from various environments worldwide. We found a close relationship between the number of interspecies recombination events within core genomes of the 35 strains and the overall genomic identity, as inferred from calculations of the average nucleotide identity. The relationship between the overall nucleotide identity and the number of detected interspecies recombination events was comparable when analyzing strains isolated over 80 years apart, from different hemispheres, or from different ecologies, as well as in strains isolated from the same geographic location within a short time frame. We further applied the same method of detecting recombination events to analyze 11 strains of Vibrio campbellii, and identified disproportionally high number of intraspecies recombination events within the core genomes of some, but not all, strains. The high number of recombination events was detected between V. campbellii strains that have significant temporal (over 18 years) and geographical (over 10,000 km) differences in their origins of isolation. Results of this study reveal a remarkable stability of Harveyi clade species, and give clues about the origins and persistence of species in the clade.

Modularized evolution in archaeal methanogens phylogenetic forest

Methanogens are methane-producing archaea that plays a key role in the global carbon cycle. To date, the evolutionary history of methanogens and closely related nonmethanogen species remains unresolved among studies conducted upon different genetic markers, attributing to horizontal gene transfers (HGTs). With an effort to decipher both congruent and conflicting evolutionary events, reconstruction of coevolved gene clusters and hierarchical structure in the archaeal methanogen phylogenetic forest, comprehensive evolution, and network analyses were performed upon 3,694 gene families from 41 methanogens and 33 closely related archaea. Our results show that 1) greater than 50% of genes are in topological dissonance with others; 2) the prevalent interorder HGTs, even for core genes, in methanogen genomes led to their scrambled phylogenetic relationships; 3) most methanogenesis-related genes have experienced at least one HGT; 4) greater than 20% of the genes in methanogen genomes were transferred horizontally from other archaea, with genes involved in cell-wall synthesis and defense system having been transferred most frequently; 5) the coevolution network contains seven statistically robust modules, wherein the central module has the highest average node strength and comprises a majority of the core genes; 6) different coevolutionary module genes boomed in different time and evolutionary lineage, constructing diversified pan-genome structures; 7) the modularized evolution is also closely related to the vertical evolution signals and the HGT rate of the genes. Overall, this study presented a modularized phylogenetic forest that describes a combination of complicated vertical and nonvertical evolutionary processes for methanogenic archaeal species.

Population genetic analysis of Streptomyces albidoflavus reveals habitat barriers to homologous recombination in the diversification of streptomycetes

Examining the population structure and the influence of recombination and ecology on microbial populations makes great sense for understanding microbial evolution and speciation. Streptomycetes are a diverse group of bacteria that are widely distributed in nature and a rich source of useful bioactive compounds; however, they are rarely subjected to population genetic investigations. In this study, we applied a five-gene-based multilocus sequence analysis (MLSA) scheme to 41 strains of Streptomyces albidoflavus derived from diverse sources, mainly insects, sea, and soil. Frequent recombination was detected in S. albidoflavus, supported by multiple lines of evidence from the pairwise homoplasy index (Φw) test, phylogenetic discordance, the Shimodaira-Hasegawa (SH) test, and network analysis, underpinning the predominance of homologous recombination within Streptomyces species. A strong habitat signal was also observed in both phylogenetic and Structure 2.3.3 analyses, indicating the importance of ecological difference in shaping the population structure. Moreover, all three habitat-associated groups, particularly the entomic group, demonstrated significantly reduced levels of gene flow with one another, generally revealing habitat barriers to recombination. Therefore, a combined effect of homologous recombination and ecology is inferred for S. albidoflavus, where dynamic evolution is at least partly balanced by the extent that differential distributions of strains among habitats limit genetic exchange. Our study stresses the significance of ecology in microbial speciation and reveals the coexistence of homologous recombination and ecological divergence in the evolution of streptomycetes.

Diversification of bacterial genome content through distinct mechanisms over different timescales

Bacterial populations often consist of multiple co-circulating lineages. Determining how such population structures arise requires understanding what drives bacterial diversification. Using 616 systematically sampled genomes, we show that Streptococcus pneumoniae lineages are typically characterized by combinations of infrequently transferred stable genomic islands: those moving primarily through transformation, along with integrative and conjugative elements and phage-related chromosomal islands. The only lineage containing extensive unique sequence corresponds to a set of atypical unencapsulated isolates that may represent a distinct species. However, prophage content is highly variable even within lineages, suggesting frequent horizontal transmission that would necessitate rapidly diversifying anti-phage mechanisms to prevent these viruses sweeping through populations. Correspondingly, two loci encoding Type I restriction-modification systems able to change their specificity over short timescales through intragenomic recombination are ubiquitous across the collection. Hence short-term pneumococcal variation is characterized by movement of phage and intragenomic rearrangements, with the slower transfer of stable loci distinguishing lineages.

Populations of the pathogenic bacterium Streptococcus pneumoniae consist of distinct co-circulating lineages. Here, the authors show lineages are characterized by particular combinations of stable genomic islands, whereas prophage and restriction-modification systems vary over short timescales.

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.

Distribution, virulence-associated genes and antimicrobial resistance of Aeromonas isolates from diarrheal patients and water, China

OBJECTIVES

To determine the prevalence of Aeromonas infections in diarrheal patients, the distribution of virulence-associated genes and antibiotic resistance among different Aeromonas species in China.

METHODS

We conducted continual active surveillance aimed on Aeromonas from diarrheal patients and aquatic samples. Aeromonas strains were identified by biochemical tests, further confirmed to species level by a multilocus phylogenetic analysis. Potential virulence genes were detected by PCR. Antibiotics susceptibility testing was carried based on the minimal inhibitory concentration.

RESULTS

From 5069 samples (stool specimens, n = 4529; water samples, n = 540) in China, 257 Aeromonas isolates [stools, n = 193 (4.3%); water, n = 64 (11.9%)] were identified by biochemical tests. The most common species from stools and water were Aeromonas veronii (42.5%) and Aeromonas caviae (37.5%), respectively. Distribution of five potential genes were significantly different between stool and water samples, two genes (ast and alt) were higher in stool than in water samples (P < 0.01). Meanwhile, three species (A. veronii, A. caviae and Aeromonas aquariorum) account for the six most prevalent combination patterns of potential genes. Furthermore, strains resistant to nine antibiotics was markedly higher in strains isolated from water than those from stools (P ≤ 0.003); in contrast, resistance to only two antibiotics was higher in strains isolated from stools compared to those from water. In addition, strains containing multiple antibiotic resistance (MAR) from stools (8.6%; 16/187) and water (30.2%; 19/63) were resistant to ten or more antibiotics.

CONCLUSION

Our study highlights the multiple factors involved in the pathogenesis of Aeromonas and reveals that environmental Aeromonas has acquired a wide range of MAR compared to those from clinical sources.

Extensive genomic variation within clonal bacterial groups resulted from homologous recombination

Due to divergence, genetic variation is generally believed to be high among distantly related strains, low among closely related ones and little or none within the same classified clonal groups. Several recent genome-wide studies, however, revealed that significant genetic variation resides in a considerable number of genes among strains with identical MLST (Multilocus sequence typing) types and much of the variation was introduced by homologous recombination. Recognizing and understanding genomic variation within clonal bacterial groups could shed new light on the evolutionary path of infectious agents and the emergence of particularly pathogenic or virulent variants. This commentary presents our recent contributions to this line of work.

Synechococcus: 3 billion years of global dominance

Cyanobacteria are among the most important primary producers on the Earth. However, the evolutionary forces driving cyanobacterial species diversity remain largely enigmatic due to both their distinction from macro-organisms and an undersampling of sequenced genomes. Thus, we present a new genome of a Synechococcus-like cyanobacterium from a novel evolutionary lineage. Further, we analyse all existing 16S rRNA sequences and genomes of Synechococcus-like cyanobacteria. Chronograms showed extremely polyphyletic relationships in Synechococcus, which has not been observed in any other cyanobacteria. Moreover, most Synechococcus lineages bifurcated after the Great Oxidation Event, including the most abundant marine picoplankton lineage. Quantification of horizontal gene transfer among 70 cyanobacterial genomes revealed significant differences among studied genomes. Horizontal gene transfer levels were not correlated with ecology, genome size or phenotype, but were correlated with the age of divergence. All findings were synthetized into a novel model of cyanobacterial evolution, characterized by serial convergence of the features, that is multicellularity and ecology.

Processing of abasic site damaged lesions by APE1 enzyme on DNA adsorbed over normal and organomodified clay

The efficiency of the apurinic/apyrimidinic endonuclease (APE1) DNA repair enzyme in the processing of abasic site DNA damage lesions at precise location in DNA oligomer duplexes that are adsorbed on clay surfaces was evaluated. Three different forms of clay namely montmorillonite, quaternary ammonium salt modified montmorillonite and its boiled counterpart i.e. partially devoid of organic moiety were used for a comparative study of adsorption, desorption and DNA repair efficiency on their surfaces. The interaction between the DNA and the clay was analysed by X-ray diffraction, Atomic force microscopy, UV-Vis spectroscopy and Infrared spectroscopy. The abasic site cleavage efficiency of APE1 enzyme was quantitatively evaluated by polyacrylamide gel electrophoresis. Apart from the difference in the DNA adsorption or desorption capacity of the various forms of clay, substantial variation in the repair efficiency of abasic sites initiated by the APE1 enzyme on the clay surfaces was observed. The incision efficiency of APE1 enzyme at abasic sites was found to be greatly diminished, when the DNA was adsorbed over organomodified montmorillonite. The reduced repair activity indicates an important role of the pendant surfactant groups on the clay surfaces in directing APE1 mediated cleavage of abasic site DNA damage lesions.

Horizontal transfer of short and degraded DNA has evolutionary implications for microbes and eukaryotic sexual reproduction

Horizontal gene transfer in the form of long DNA fragments has changed our view of bacterial evolution. Recently, we discovered that such processes may also occur with the massive amounts of short and damaged DNA in the environment, and even with truly ancient DNA. Although it presently remains unclear how often it takes place in nature, horizontal gene transfer of short and damaged DNA opens up the possibility for genetic exchange across distinct species in both time and space. In this essay, we speculate on the potential evolutionary consequences of this phenomenon. We argue that it may challenge basic assumptions in evolutionary theory; that it may have distant origins in life's history; and that horizontal gene transfer should be viewed as an evolutionary strategy not only preceding but causally underpinning the evolution of sexual reproduction.

An In-depth Analysis of a Multilocus Phylogeny Identifies leuS As a Reliable Phylogenetic Marker for the Genus Pantoea

Partial sequences of six core genes (fusA, gyrB, leuS, pyrG, rlpB, and rpoB) of 37 strains of Pantoea species were analyzed in order to obtain a comprehensive view regarding the phylogenetic relationships within the Pantoea genus and compare tree topologies to identify gene(s) for reliable species and subspecies differentiation. All genes used in this study were effective at species-level delineation, but the internal nodes represented conflicting common ancestors in fusA- and pyrG-based phylogenies. Concatenated gene phylogeny gave the expected DNA relatedness, underscoring the significance of a multilocus sequence analysis. Pairwise comparison of topological distances and percent similarities indicated a significant differential influence of individual genes on the concatenated tree topology. leuS- and fusA-inferred phylogenies exhibited, respectively, the lowest (4) and highest (52) topological distances to the concatenated tree. These correlated well with high (96.3%) and low (64.4%) percent similarities of leuS- and fusA-inferred tree topologies to the concatenated tree, respectively. We conclude that the concatenated tree topology is strongly influenced by the gene with the highest number of polymorphic and non-synonymous sites in the absence of significant recombination events.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes

Background

The bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes.

Description

We have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers.

Conclusions

With a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.

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.

Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus

Horizontal gene transfer is an important driver of bacterial evolution, but genetic exchange in the core genome of clonal species, including the major pathogen Staphylococcus aureus, is incompletely understood. Here we reveal widespread homologous recombination in S. aureus at the species level, in contrast to its near-complete absence between closely related strains. We discover a patchwork of hotspots and coldspots at fine scales falling against a backdrop of broad-scale trends in rate variation. Over megabases, homoplasy rates fluctuate 1.9-fold, peaking towards the origin-of-replication. Over kilobases, we find core recombination hotspots of up to 2.5-fold enrichment situated near fault lines in the genome associated with mobile elements. The strongest hotspots include regions flanking conjugative transposon ICE6013, the staphylococcal cassette chromosome (SCC) and genomic island νSaα. Mobile element-driven core genome transfer represents an opportunity for adaptation and challenges our understanding of the recombination landscape in predominantly clonal pathogens, with important implications for genotype–phenotype mapping.

Horizontal gene transfer occurs in most bacteria, yet it is unclear whether it happens in clonal species. Here, Everitt et al. show widespread within-species recombination, driven by mobile elements, in the genome of the pathogen Staphylococcus aureus, but no recombination between closely related strains.

Cryptic ecology among host generalist Campylobacter jejuni in domestic animals

Homologous recombination between bacterial strains is theoretically capable of preventing the separation of daughter clusters, and producing cohesive clouds of genotypes in sequence space. However, numerous barriers to recombination are known. Barriers may be essential such as adaptive incompatibility, or ecological, which is associated with the opportunities for recombination in the natural habitat. Campylobacter jejuni is a gut colonizer of numerous animal species and a major human enteric pathogen. We demonstrate that the two major generalist lineages of C. jejuni do not show evidence of recombination with each other in nature, despite having a high degree of host niche overlap and recombining extensively with specialist lineages. However, transformation experiments show that the generalist lineages readily recombine with one another in vitro. This suggests ecological rather than essential barriers to recombination, caused by a cryptic niche structure within the hosts.

Evolutionary genomics of Borrelia burgdorferi sensu lato: findings, hypotheses, and the rise of hybrids

Borrelia burgdorferi sensu lato (B. burgdorferi s.l.), the group of bacterial species represented by Lyme disease pathogens, has one of the most complex and variable genomic architectures among prokaryotes. Showing frequent recombination within and limited gene flow among geographic populations, the B. burgdorferi s.l. genomes provide an excellent window into the processes of bacterial evolution at both within- and between-population levels. Comparative analyses of B. burgdorferi s.l. genomes revealed a highly dynamic plasmid composition but a conservative gene repertoire. Gene duplication and loss as well as sequence variations at loci encoding surface-localized lipoproteins (e.g., the PF54 genes) are strongly associated with adaptive differences between species. There are a great many conserved intergenic spacer sequences that are candidates for cis-regulatory elements and non-coding RNAs. Recombination among coexisting strains occurs at a rate approximately three times the mutation rate. The coexistence of a large number of genomic groups within local B. burgdorferi s.l. populations may be driven by immune-mediated diversifying selection targeting major antigen loci as well as by adaptation to multiple host species. Questions remain regarding the ecological causes (e.g., climate change, host movements, or new adaptations) of the ongoing range expansion of B. burgdorferi s.l. and on the genomic variations associated with its ecological and clinical variability. Anticipating an explosive growth of the number of B. burgdorferi s.l. genomes sampled from both within and among species, we propose genome-based methods to test adaptive mechanisms and to identify molecular bases of phenotypic variations. Genome sequencing is also necessary for monitoring a likely increase of genetic admixture of previously isolated species and populations in North America and elsewhere.

Phylogeny poorly predicts the utility of a challenging horizontally transferred gene in Methylobacterium strains

Horizontal gene transfer plays a crucial role in microbial evolution. While much is known about the mechanisms that determine whether physical DNA can be transferred into a new host, the factors determining the utility of the transferred genes are less clear. We have explored this issue using dichloromethane consumption in Methylobacterium strains. Methylobacterium extorquens DM4 expresses a dichloromethane dehalogenase (DcmA) that has been acquired through horizontal gene transfer and allows the strain to grow on dichloromethane as the sole carbon and energy source. We transferred the dcmA gene into six Methylobacterium strains that include both close and distant evolutionary relatives. The transconjugants varied in their ability to grow on dichloromethane, but their fitness on dichloromethane did not correlate with the phylogeny of the parental strains or with any single tested physiological factor. This work highlights an important limiting factor in horizontal gene transfer, namely, the capacity of the recipient strain to accommodate the stress and metabolic disruption resulting from the acquisition of a new enzyme or pathway. Understanding these limitations may help to rationalize historical examples of horizontal transfer and aid deliberate genetic transfers in biotechnology for metabolic engineering.

A reference pan-genome approach to comparative bacterial genomics: identification of novel epidemiological markers in pathogenic Campylobacter

The increasing availability of hundreds of whole bacterial genomes provides opportunities for enhanced understanding of the genes and alleles responsible for clinically important phenotypes and how they evolved. However, it is a significant challenge to develop easy-to-use and scalable methods for characterizing these large and complex data and relating it to disease epidemiology. Existing approaches typically focus on either homologous sequence variation in genes that are shared by all isolates, or non-homologous sequence variation--focusing on genes that are differentially present in the population. Here we present a comparative genomics approach that simultaneously approximates core and accessory genome variation in pathogen populations and apply it to pathogenic species in the genus Campylobacter. A total of 7 published Campylobacter jejuni and Campylobacter coli genomes were selected to represent diversity across these species, and a list of all loci that were present at least once was compiled. After filtering duplicates a 7-isolate reference pan-genome, of 3,933 loci, was defined. A core genome of 1,035 genes was ubiquitous in the sample accounting for 59% of the genes in each isolate (average genome size of 1.68 Mb). The accessory genome contained 2,792 genes. A Campylobacter population sample of 192 genomes was screened for the presence of reference pan-genome loci with gene presence defined as a BLAST match of ≥ 70% identity over ≥ 50% of the locus length--aligned using MUSCLE on a gene-by-gene basis. A total of 21 genes were present only in C. coli and 27 only in C. jejuni, providing information about functional differences associated with species and novel epidemiological markers for population genomic analyses. Homologs of these genes were found in several of the genomes used to define the pan-genome and, therefore, would not have been identified using a single reference strain approach.

Ordering microbial diversity into ecologically and genetically cohesive units

We propose that microbial diversity must be viewed in light of gene flow and selection, which define units of genetic similarity, and of phenotype and ecological function, respectively. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and on evidence from some of the first microbial populations studied with genomics. These show that if recombination is frequent and selection moderate, ecologically adaptive mutations or genes can spread within populations independently of their original genomic background (gene-specific sweeps). Alternatively, if the effect of recombination is smaller than selection, genome-wide selective sweeps should occur. In both cases, however, distinct units of overlapping ecological and genotypic similarity will form if microgeographic separation, likely involving ecological tradeoffs, induces barriers to gene flow. These predictions are supported by (meta)genomic data, which suggest that a 'reverse ecology' approach, in which genomic and gene flow information is used to make predictions about the nature of ecological units, is a powerful approach to ordering microbial diversity.

Explaining microbial genomic diversity in light of evolutionary ecology

Comparisons of closely related microorganisms have shown that individual genomes can be highly diverse in terms of gene content. In this Review, we discuss several studies showing that much of this variation is associated with social and ecological interactions, which have an important role in the population biology of wild populations of bacteria and archaea. These interactions create frequency-dependent selective pressures that can either stabilize gene frequencies at intermediate levels in populations or promote fast gene turnover, which presents as low gene frequencies in genome surveys. Thus, interpretation of gene-content diversity requires the delineation of populations according to cohesive gene flow and ecology, as micro-evolutionary changes arise in response to local selection pressures and population dynamics.

Molecular phylogenetics and temporal diversification in the genus Aeromonas based on the sequences of five housekeeping genes

Several approaches have been developed to estimate both the relative and absolute rates of speciation and extinction within clades based on molecular phylogenetic reconstructions of evolutionary relationships, according to an underlying model of diversification. However, the macroevolutionary models established for eukaryotes have scarcely been used with prokaryotes. We have investigated the rate and pattern of cladogenesis in the genus Aeromonas (γ-Proteobacteria, Proteobacteria, Bacteria) using the sequences of five housekeeping genes and an uncorrelated relaxed-clock approach. To our knowledge, until now this analysis has never been applied to all the species described in a bacterial genus and thus opens up the possibility of establishing models of speciation from sequence data commonly used in phylogenetic studies of prokaryotes. Our results suggest that the genus Aeromonas began to diverge between 248 and 266 million years ago, exhibiting a constant divergence rate through the Phanerozoic, which could be described as a pure birth process.

Insertion sequence distribution bias in Archaea

Insertion sequences (IS) are common transposable elements in Archaea. Intergenic IS elements are usually less harmful than intragenic ISs, simply because they are less likely to disrupt host gene function. However, because regulatory sequences are intergenic and upstream of genes, we hypothesized that not all intergenic regions are selectively equivalent for IS insertion. We tested this hypothesis by analyzing the distributions of intergenic IS elements within 155 fully sequenced archaeal genomes. Of the 22 genomes with enough IS elements for statistical analysis, five have significantly fewer ISs between divergently oriented neighboring genes than expected by chance, and seven have significantly more ISs between convergently oriented genes. Furthermore, of the 85 genomes with at least one expected IS within each of the three possible neighboring gene orientations (i.e., divergent, convergent, and tandem), 73 genomes have fewer ISs between divergently oriented genes than expected, and 60 have more ISs between convergently oriented genes than expected (both values deviate significantly from binomial probabilities of random distribution). We suspect that these non-random IS distributions are molded by natural selection resulting from differential disruption of neighboring gene regulation, and that this selective pressure has affected transposable element distributions in prokaryotes for billions of years.

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.

Harnessing population genomics to understand how bacterial pathogens emerge, adapt to crop hosts, and disseminate

Crop diseases emerge without warning. In many cases, diseases cross borders, or even oceans, before plant pathologists have time to identify and characterize the causative agents. Genome sequencing, in combination with intensive sampling of pathogen populations and application of population genetic tools, is now providing the means to unravel how bacterial crop pathogens emerge from environmental reservoirs, how they evolve and adapt to crops, and what international and intercontinental routes they follow during dissemination. Here, we introduce the field of population genomics and review the population genomics research of bacterial plant pathogens over the past 10 years. We highlight the potential of population genomics for investigating plant pathogens, using examples of population genomics studies of human pathogens. We also describe the complementary nature of the fields of population genomics and molecular plant-microbe interactions and propose how to translate new insights into improved disease prevention and control.

Frequent recombination shapes the epidemic population structure of Planktothrix (Cyanoprokaryota) in Italian subalpine lakes

The planktonic genus Planktothrix, as other cyanobacteria, shows signals of both homologous and nonhomologous recombination. However, the frequency of recombination and its effect on Planktothrix population structuring is unknown. We isolated 290 Planktothrix strains from seven neighboring lakes in the subalpine Italian region and analyzed these using multilocus sequence typing. Four of six loci analyzed were polymorphic, resulting in 20 distinct multilocus genotypes. Association indices among alleles at different loci were suggestive of an "epidemic population structure," resulting from an explosive (and temporary) dominance of one genotype against a panmictic background. ClonalFrame analyses supported this view by detecting: (i) three major clades affected by three distinct recombination events, (ii) a recombination rate about equal to the mutation rate, and (iii) the fact that recombination had an impact on introducing molecular diversity more than double the mutation rate. Furthermore, analysis of molecular variance over an annual cycle in three of seven lakes revealed that both local clonal expansion and recombination processes affected among-lake diversity. Our observations suggest that recombination affects microevolution of Planktothrix and that an epidemic structure can emerge in populations of this genus.