Gradual evolution in bacteria: evidence from Bacillus systematics

Environment Shapes the Intra-species Diversity of Bacillus subtilis Isolates

Cosmopolitan bacteria are those that are found practically everywhere in the world. One of them is Bacillus subtilis, which can travel around the world through dust storms rising from various deserts. Upon landing, bacterial survival is determined by the ability to adjust to the heterogonous environments and bacteria isolated from extremely different environments, such as desert and riverbank soil, are expected to be less related due to the environmental pressure of each region. However, little is known about the influence of soil and habitat on B. subtilis evolution. Here, we show that desert and riverbank B. subtilis strains differ in genetic relatedness and physiological traits, such as biofilm morphology and utilisation of carbon sources. Desert strains showed more diversity at the genetic level and were able to utilise more carbon sources than riverbank strains which were highly genetically conserved. Biofilm morphologies of desert and riverbank strains generally segregated and both groups formed different morphology clusters despite the astonishing diversity observed among riverbank strains. We also show that relatedness of B. subtilis strains does not decrease with distance inside the same habitat, which, together with diversity data implies that the difference in environmental selection pressures plays a fundamental role in the evolution of this species.

Rokubacteria: Genomic Giants among the Uncultured Bacterial Phyla

Recent advances in single-cell genomic and metagenomic techniques have facilitated the discovery of numerous previously unknown, deep branches of the tree of life that lack cultured representatives. Many of these candidate phyla are composed of microorganisms with minimalistic, streamlined genomes lacking some core metabolic pathways, which may contribute to their resistance to growth in pure culture. Here we analyzed single-cell genomes and metagenome bins to show that the "Candidate phylum Rokubacteria," formerly known as SPAM, represents an interesting exception, by having large genomes (6-8 Mbps), high GC content (66-71%), and the potential for a versatile, mixotrophic metabolism. We also observed an unusually high genomic heterogeneity among individual Rokubacteria cells in the studied samples. These features may have contributed to the limited recovery of sequences of this candidate phylum in prior cultivation and metagenomic studies. Our analyses suggest that Rokubacteria are distributed globally in diverse terrestrial ecosystems, including soils, the rhizosphere, volcanic mud, oil wells, aquifers, and the deep subsurface, with no reports from marine environments to date.

An Inducible Operon Is Involved in Inulin Utilization in Lactobacillus plantarum Strains, as Revealed by Comparative Proteogenomics and Metabolic Profiling

The draft genomes of Lactobacillus plantarum strains isolated from Asian fermented foods, infant feces, and shrimp intestines were sequenced and compared to those of well-studied strains. Among 28 strains of L. plantarum, variations in the genomic features involved in ecological adaptation were elucidated. The genome sizes ranged from approximately 3.1 to 3.5 Mb, of which about 2,932 to 3,345 protein-coding sequences (CDS) were predicted. The food-derived isolates contained a higher number of carbohydrate metabolism-associated genes than those from infant feces. This observation correlated to their phenotypic carbohydrate metabolic profile, indicating their ability to metabolize the largest range of sugars. Surprisingly, two strains (P14 and P76) isolated from fermented fish utilized inulin. β-Fructosidase, the inulin-degrading enzyme, was detected in the supernatants and cell wall extracts of both strains. No activity was observed in the cytoplasmic fraction, indicating that this key enzyme was either membrane-bound or extracellularly secreted. From genomic mining analysis, a predicted inulin operon of fosRABCDXE, which encodes β-fructosidase and many fructose transporting proteins, was found within the genomes of strains P14 and P76. Moreover, pts1BCA genes, encoding sucrose-specific IIBCA components involved in sucrose transport, were also identified. The proteomic analysis revealed the mechanism and functional characteristic of the fosRABCDXE operon involved in the inulin utilization of L. plantarum The expression levels of the fos operon and pst genes were upregulated at mid-log phase. FosE and the LPXTG-motif cell wall anchored β-fructosidase were induced to a high abundance when inulin was present as a carbon source.

IMPORTANCE

Inulin is a long-chain carbohydrate that may act as a prebiotic, which provides many health benefits to the host by selectively stimulating the growth and activity of beneficial bacteria in the colon. While certain lactobacilli can catabolize inulin, this has not yet been described for Lactobacillus plantarum, and an associated putative inulin operon has not been reported in this species. By using comparative and functional genomics, we showed that two L. plantarum strains utilized inulin and identified functional inulin operons in their genomes. The proteogenomic data revealed that inulin degradation and uptake routes, which related to the fosRABCDXE operon and pstBCA genes, were widely expressed among L. plantarum strains. The present work provides a novel understanding of gene regulation and mechanisms of inulin utilization in probiotic L. plantarum generating opportunities for synbiotic product development.

Analysis of gene gain and loss in the evolution of predatory bacteria

Predatory bacteria are ubiquitously distributed in nature in including in aquatic environments, sewage, intestinal tracts of animals and humans, rhizophere and, soils. However, our understanding of their evolutionary history is limited. Results of recent studies have shown that acquiring novel genes is a major force driving bacterial evolution. Therefore, to gain a better understanding of the impact of gene gain and loss in the evolution of bacterial predators, this study employed comparative genomic approaches to identify core-set gene families and species-specific gene families, and model gene gain and loss events among 11 genomes that represented diverse lineages. In total, 1977 gene families were classified. Of these 509 (pattern 11111111111) were present all of the 11 species. Among the non-core set gene families, 52 were present only in saltwater bacteria predators and had no ortholog in the other genomes. Similarly 109 and 44 were present only in the genomes of Micavibrio spp. and Bdellovibrio spp., respectively. In this study, the gain loss mapping engine GLOOME was selected to analyze and estimate the expectations and probabilities of both gain and loss events in the predatory bacteria. In total, 354 gene families were involved in significant gene gain events, and 407 gene families were classified into gene loss events with high supported value. Moreover, 18 families from the core set gene family were identified as putative genes under positive selection. The results of this study suggest that acquisition of particular genes that encode functional proteins in metabolism and cellular processes and signaling, especially ABC systems, may help bacterial predators adapt to surrounding environmental changes and present different predation strategies for survival in their habitats.

Gene content dissimilarity for subclassification of highly similar microbial strains

Background

Identification and classification of highly similar microbial strains is a challenging issue in microbiology, ecology and evolutionary biology. Among various available approaches, gene content analysis is also at the core of microbial taxonomy. However, no threshold has been determined for grouping microorgnisms to different taxonomic levels, and it is still not clear that to what extent genomic fluidity should occur to form a microbial taxonomic group.

Results

By taking advantage of the eggNOG database for orthologous groups, we calculated gene content dissimilarity among different microbial strains based on the orthologous gene profiles and tested the possibility of applying gene content dissimilarity as a quantitative index in classifying microbial taxonomic groups, as well as its potential application in subclassification of highly similar microbial strains. Evaluation of gene content dissimilarity to completed microbial genomes at different taxonomic levels suggested that cutoffs of 0.2 and 0.4 can be respectively used for species and family delineation, and that 0.2 gene content dissimilarity cutoff approximately corresponded to 98 % 16S rRNA gene identity and 94 % ANI for microbial species delineation. Furthermore, application of gene content dissimilarity to highly similar microbial strains suggested it as an effective approach in classifying closely related microorganisms into subgroups.

Conclusions

This approach is especially useful in identifying pathogens from commensals in clinical microbiology. It also provides novel insights into how genomic fluidity is linked with microbial taxonomy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2991-9) contains supplementary material, which is available to authorized users.

Variability of rRNA Operon Copy Number and Growth Rate Dynamics of Bacillus Isolated from an Extremely Oligotrophic Aquatic Ecosystem

The ribosomal RNA (rrn) operon is a key suite of genes related to the production of protein synthesis machinery and thus to bacterial growth physiology. Experimental evidence has suggested an intrinsic relationship between the number of copies of this operon and environmental resource availability, especially the availability of phosphorus (P), because bacteria that live in oligotrophic ecosystems usually have few rrn operons and a slow growth rate. The Cuatro Ciénegas Basin (CCB) is a complex aquatic ecosystem that contains an unusually high microbial diversity that is able to persist under highly oligotrophic conditions. These environmental conditions impose a variety of strong selective pressures that shape the genome dynamics of their inhabitants. The genus Bacillus is one of the most abundant cultivable bacterial groups in the CCB and usually possesses a relatively large number of rrn operon copies (6–15 copies). The main goal of this study was to analyze the variation in the number of rrn operon copies of Bacillus in the CCB and to assess their growth-related properties as well as their stoichiometric balance (N and P content). We defined 18 phylogenetic groups within the Bacilli clade and documented a range of from six to 14 copies of the rrn operon. The growth dynamic of these Bacilli was heterogeneous and did not show a direct relation to the number of operon copies. Physiologically, our results were not consistent with the Growth Rate Hypothesis, since the copies of the rrn operon were decoupled from growth rate. However, we speculate that the diversity of the growth properties of these Bacilli as well as the low P content of their cells in an ample range of rrn copy number is an adaptive response to oligotrophy of the CCB and could represent an ecological mechanism that allows these taxa to coexist. These findings increase the knowledge of the variability in the number of copies of the rrn operon in the genus Bacillus and give insights about the physiology of this bacterial group under extreme oligotrophic conditions.

Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis

Closely related bacterial genomes usually differ in gene content, suggesting that nearly every strain in nature may be ecologically unique. We have tested this hypothesis by sequencing the genomes of extremely close relatives within a recognized taxon and analyzing the genomes for evidence of ecological distinctness. We compared the genomes of four Death Valley isolates plus the laboratory strain W23, all previously classified as Bacillus subtilis subsp. spizizenii and hypothesized through multilocus analysis to be members of the same ecotype (an ecologically homogeneous population), named putative ecotype 15 (PE15). These strains showed a history of positive selection on amino acid sequences in 38 genes. Each of the strains was under a different regimen of positive selection, suggesting that each strain is ecologically unique and represents a distinct ecological speciation event. The rate of speciation appears to be much faster than can be resolved with multilocus sequencing. Each PE15 strain contained unique genes known to confer a function for bacteria. Remarkably, no unique gene conferred a metabolic system or subsystem function that was not already present in all the PE15 strains sampled. Thus, the origin of ecotypes within this clade shows no evidence of qualitative divergence in the set of resources utilized. Ecotype formation within this clade is consistent with the nanoniche model of bacterial speciation, in which ecotypes use the same set of resources but in different proportions, and genetic cohesion extends beyond a single ecotype to the set of ecotypes utilizing the same resources.

Marine sediment-derived Streptomyces bacteria from British Columbia, Canada are a promising microbiota resource for the discovery of antimicrobial natural products

Representatives of the genus Streptomyces from terrestrial sources have been the focus of intensive research for the last four decades because of their prolific production of chemically diverse and biologically important compounds. However, metabolite research from this ecological niche had declined significantly in the past years because of the rediscovery of the same bioactive compounds and redundancy of the sample strains. More recently, a new picture has begun to emerge in which marine-derived Streptomyces bacteria have become the latest hot spot as new source for unique and biologically active compounds. Here, we investigated the marine sediments collected in the temperate cold waters from British Columbia, Canada as a valuable source for new groups of marine-derived Streptomyces with antimicrobial activities. We performed culture dependent isolation from 49 marine sediments samples and obtained 186 Streptomyces isolates, 47 of which exhibited antimicrobial activities. Phylogenetic analyses of the active isolates resulted in the identification of four different clusters of bioactive Streptomyces including a cluster with isolates that appear to represent novel species. Moreover, we explored whether these marine-derived Streptomyces produce new secondary metabolites with antimicrobial properties. Chemical analyses revealed structurally diverse secondary metabolites, including four new antibacterial novobiocin analogues. We conducted structure-activity relationships (SAR) studies of these novobiocin analogues against methicillin-resistant Staphylococcus aureus (MRSA). In this study, we revealed the importance of carbamoyl and OMe moieties at positions 3” and 4” of novobiose as well as the hydrogen substituent at position 5 of hydroxybenzoate ring for the anti-MRSA activity. Changes in the substituents at these positions dramatically impede or completely eliminate the inhibitory activity of novobiocins against MRSA.

Analysis of the Lactobacillus casei supragenome and its influence in species evolution and lifestyle adaptation

Background

The broad ecological distribution of L. casei makes it an insightful subject for research on genome evolution and lifestyle adaptation. To explore evolutionary mechanisms that determine genomic diversity of L. casei, we performed comparative analysis of 17 L. casei genomes representing strains collected from dairy, plant, and human sources.

Results

Differences in L. casei genome inventory revealed an open pan-genome comprised of 1,715 core and 4,220 accessory genes. Extrapolation of pan-genome data indicates L. casei has a supragenome approximately 3.2 times larger than the average genome of individual strains. Evidence suggests horizontal gene transfer from other bacterial species, particularly lactobacilli, has been important in adaptation of L. casei to new habitats and lifestyles, but evolution of dairy niche specialists also appears to involve gene decay.

Conclusions

Genome diversity in L. casei has evolved through gene acquisition and decay. Acquisition of foreign genomic islands likely confers a fitness benefit in specific habitats, notably plant-associated niches. Loss of unnecessary ancestral traits in strains collected from bacterial-ripened cheeses supports the hypothesis that gene decay contributes to enhanced fitness in that niche. This study gives the first evidence for a L. casei supragenome and provides valuable insights into mechanisms for genome evolution and lifestyle adaptation of this ecologically flexible and industrially important lactic acid bacterium. Additionally, our data confirm the Distributed Genome Hypothesis extends to non-pathogenic, ecologically flexible species like L. casei.

Comparative genomic analyses of Streptococcus mutans provide insights into chromosomal shuffling and species-specific content

Background

Streptococcus mutans is the major pathogen of dental caries, and it occasionally causes infective endocarditis. While the pathogenicity of this species is distinct from other human pathogenic streptococci, the species-specific evolution of the genus Streptococcus and its genomic diversity are poorly understood.

Results

We have sequenced the complete genome of S. mutans serotype c strain NN2025, and compared it with the genome of UA159. The NN2025 genome is composed of 2,013,587 bp, and the two strains show highly conserved core-genome. However, comparison of the two S. mutans strains showed a large genomic inversion across the replication axis producing an X-shaped symmetrical DNA dot plot. This phenomenon was also observed between other streptococcal species, indicating that streptococcal genetic rearrangements across the replication axis play an important role in Streptococcus genetic shuffling. We further confirmed the genomic diversity among 95 clinical isolates using long-PCR analysis. Genomic diversity in S. mutans appears to occur frequently between insertion sequence (IS) elements and transposons, and these diversity regions consist of restriction/modification systems, antimicrobial peptide synthesis systems, and transporters. S. mutans may preferentially reject the phage infection by clustered regularly interspaced short palindromic repeats (CRISPRs). In particular, the CRISPR-2 region, which is highly divergent between strains, in NN2025 has long repeated spacer sequences corresponding to the streptococcal phage genome.

Conclusion

These observations suggest that S. mutans strains evolve through chromosomal shuffling and that phage infection is not needed for gene acquisition. In contrast, S. pyogenes tolerates phage infection for acquisition of virulence determinants for niche adaptation.

Genome sequence and comparative genome analysis of Lactobacillus casei: insights into their niche-associated evolution

Lactobacillus casei is remarkably adaptable to diverse habitats and widely used in the food industry. To reveal the genomic features that contribute to its broad ecological adaptability and examine the evolution of the species, the genome sequence of L. casei ATCC 334 is analyzed and compared with other sequenced lactobacilli. This analysis reveals that ATCC 334 contains a high number of coding sequences involved in carbohydrate utilization and transcriptional regulation, reflecting its requirement for dealing with diverse environmental conditions. A comparison of the genome sequences of ATCC 334 to L. casei BL23 reveals 12 and 19 genomic islands, respectively. For a broader assessment of the genetic variability within L. casei, gene content of 21 L. casei strains isolated from various habitats (cheeses, n = 7; plant materials, n = 8; and human sources, n = 6) was examined by comparative genome hybridization with an ATCC 334-based microarray. This analysis resulted in identification of 25 hypervariable regions. One of these regions contains an overrepresentation of genes involved in carbohydrate utilization and transcriptional regulation and was thus proposed as a lifestyle adaptation island. Differences in L. casei genome inventory reveal both gene gain and gene decay. Gene gain, via acquisition of genomic islands, likely confers a fitness benefit in specific habitats. Gene decay, that is, loss of unnecessary ancestral traits, is observed in the cheese isolates and likely results in enhanced fitness in the dairy niche. This study gives the first picture of the stable versus variable regions in L. casei and provides valuable insights into evolution, lifestyle adaptation, and metabolic diversity of L. casei.

Looking for Darwin's footprints in the microbial world

As we observe the 200th anniversary of Charles Darwin's birth, microbiologists interested in the application of Darwin's ideas to the microscopic world have a lot to celebrate: an emerging picture of the (mostly microbial) Tree of Life at ever-increasing resolution, an understanding of horizontal gene transfer as a driving force in the evolution of microbes, and thousands of complete genome sequences to help formulate and refine our theories. At the same time, quantitative models of the microevolutionary processes shaping microbial populations remain just out of reach, a point that is perhaps most dramatically illustrated by the lack of consensus on how (or even whether) to define bacterial species. Here, we summarize progress and prospects in bacterial population genetics, with an emphasis on detecting the footprint of positive Darwinian selection in microbial genomes.

Phylogenetic Relationships and Functional Genes: Distribution of a Gene (mnxG) encoding a putative manganese-oxidizing enzyme in Bacillus species

Several Bacillus and Paenibacillus species were isolated from Fe and Mn oxide minerals precipitating at a deep subsurface oxic-anoxic interface at Henderson Molybdenum Mine, Empire, CO. The isolates were investigated for their Mn(II)-oxidizing potential and interrogated for possession of the mnxG gene, a gene that codes for a putative Mn(II)-oxidizing enzyme in Bacillus species. Seven of eight Bacillus species were capable of Mn(II) oxidation; however, the mnxG gene was detected in only one isolate. Using sequences of known Bacillus species both with and without amplifiable mnxG genes and Henderson Mine isolates, the 16S rRNA and mnxG gene phylogenies were compared to determine if 16S rRNA sequences could be used to predict the presence or absence of an amplifiable mnxG gene within the genomes of the isolates. We discovered a strong correspondence between 16S rRNA sequence similarity and the presence/absence of an amplifiable mnxG gene in the isolates. The data revealed a complex phylogenetic distribution of the mnxG gene in which vertical inheritance and gene loss influence the distribution of the gene among the Bacillus species included in this study. Comparisons of 16S rRNA and functional gene phylogenies can be used as a tool to aid in unraveling the history and dispersal of the mnxG gene within the Bacillus clade.

Gene amelioration demonstrated: the journey of nascent genes in bacteria

Gene amelioration is the hypothesis that genes acquired via lateral gene transfer will, over time, acquire the molecular characteristics of the host genome. Species for which multiple strains have been sequenced permit a demonstration that this hypothesis is correct. We use 7 sequenced genomes of Streptococcus pyogenes and 6 sequenced genomes of Staphylococcus aureus to illustrate the action of amelioration on these genomes.

The rise and fall of deleterious mutation

It is well established that selection is less efficient in small populations than in large ones. Here we review the impact of this effect by considering the gradual selective purging of deleterious mutation over time. We outline an approach to explore the dynamics of this process, and highlight its profound implications.

A Systematics for Discovering the Fundamental Units of Bacterial Diversity

Bacterial systematists face unique challenges when trying to identify ecologically meaningful units of biological diversity. Whereas plant and animal systematists are guided by a theory-based concept of species, microbiologists have yet to agree upon a set of ecological and evolutionary properties that will serve to define a bacterial species. Advances in molecular techniques have given us a glimpse of the tremendous diversity present within the microbial world, but significant work remains to be done in order to understand the ecological and evolutionary dynamics that can account for the origin, maintenance, and distribution of that diversity. We have developed a conceptual framework that uses ecological and evolutionary theory to identify the DNA sequence clusters most likely corresponding to the fundamental units of bacterial diversity. Taking into account diverse models of bacterial evolution, we argue that bacterial systematics should seek to identify ecologically distinct groups with evidence of a history of coexistence, as based on interpretation of sequence clusters. This would establish a theory-based species unit that holds the dynamic properties broadly attributed to species outside of microbiology.

The role of laterally transferred genes in adaptive evolution

Background

Bacterial genomes develop new mechanisms to tide them over the imposing conditions they encounter during the course of their evolution. Acquisition of new genes by lateral gene transfer may be one of the dominant ways of adaptation in bacterial genome evolution. Lateral gene transfer provides the bacterial genome with a new set of genes that help it to explore and adapt to new ecological niches.

Methods

A maximum likelihood analysis was done on the five sequenced corynebacterial genomes to model the rates of gene insertions/deletions at various depths of the phylogeny.

Results

The study shows that most of the laterally acquired genes are transient and the inferred rates of gene movement are higher on the external branches of the phylogeny and decrease as the phylogenetic depth increases. The newly acquired genes are under relaxed selection and evolve faster than their older counterparts. Analysis of some of the functionally characterised LGTs in each species has indicated that they may have a possible adaptive role.

Conclusion

The five Corynebacterial genomes sequenced to date have evolved by acquiring between 8 – 14% of their genomes by LGT and some of these genes may have a role in adaptation.

Towards a conceptual and operational union of bacterial systematics, ecology, and evolution

To completely understand the ecology of a bacterial community, we need to identify its ecologically distinct populations (ecotypes). The greatest promise for enumerating a community's constituent ecotypes is held by molecular approaches that identify bacterial ecotypes as DNA sequence clusters. These approaches succeed when ecotypes correspond with sequence clusters, but some models of bacterial speciation predict a one-to-many and others a many-to-one relationship between ecotypes and sequence clusters. A further challenge is that sequence-based phylogenies often contain a hierarchy of clusters and subclusters within clusters, and there is no widely accepted theory to guide systematists and ecologists to the size of cluster most likely to correspond to ecotypes. While present systematics attempts to use universal thresholds of sequence divergence to help demarcate species, the recently developed 'community phylogeny' approach assumes no universal thresholds, but demarcates ecotypes based on the analysis of a lineage's evolutionary dynamics. Theory-based approaches like this one can give a conceptual framework as well as operational criteria for hypothesizing the identity and membership of ecotypes from sequence data; ecology-based approaches can then confirm that the putative ecotypes are actually ecologically distinct. Bacterial ecotypes that are demonstrated to have a history of coexistence as ecologically distinct lineages (based on sequence analysis) and as a prognosis of future coexistence (based on ecological differences), are the fundamental units of bacterial ecology and evolution, and should be recognized by bacterial systematics.

Culture isolation and culture-independent clone libraries reveal new marine Synechococcus ecotypes with distinctive light and N physiologies

Marine microbial communities often contain multiple closely related phylogenetic clades, but in many cases, it is still unclear what physiological traits differentiate these putative ecotypes. The numerically abundant marine cyanobacterium Synechococcus can be divided into at least 14 clades. In order to better understand ecotype differentiation in this genus, we assessed the diversity of a Synechococcus community from a well-mixed water column in the Sargasso Sea during March 2002, a time of year when this genus typically reaches its annual peak in abundance. Diversity was estimated from water sampled at three depths (approximately 5, 70, and 170 m) using both culture isolation and construction of cyanobacterial 16S-23S rRNA internal transcribed sequence clone libraries. Clonal isolates were obtained by enrichment with ammonium, nitrite, or nitrate as the sole N source, followed by pour plating. Each method sampled the in situ diversity differently. The combined methods revealed a total of seven Synechococcus phylotypes including two new putative ecotypes, labeled XV and XVI. Although most other isolates grow on nitrate, clade XV exhibited a reduced efficiency in nitrate utilization, and both clade XV and XVI are capable of chromatic adaptation, demonstrating that this trait is more widely distributed among Synechococcus strains than previously known. Thus, as in its sister genus Prochlorococcus, light and nitrogen utilization are important factors in ecotype differentiation in the marine Synechococcus lineage.

A cryptic Bacillus isolate exhibited narrow 16S rRNA gene sequence divergence with Bacillus thuringiensis and showed low maintenance requirements in hyper-osmotic complex substrate cultivations

A cryptic Bacillus (K90) isolate obtained from soil samples from the Kuwait desert exhibited lower maintenance requirements in complex substrate cultivations than Bacillus thuringiensis. A mathematical model was used to estimate apparent maintenance coefficients (m(c)) and these were found to be 0.336 and 0.041/h for B. thuringiensis and K90, respectively. The results also showed that the values of apparent maintenance coefficients were inversely related to the specific growth rates. Furthermore, 16S rRNA gene sequencing showed that K90 exhibited 99.81% sequence similarity to that of B. mojavensis and 92.9% with B. thuringiensis. It is evident from the dendrogram that the evolution of B. mojavensis (K90) (B. subtilis group), which may have originated after B. licheniformis could have been influenced by prolonged hyper-osmotic conditions, while B. thuringiensis that evolved before B. oleronius exhibited greater sensitivity as implied by the higher maintenance coefficient obtained for the hyper-osmotic cultures. As K90 exhibited low maintenance requirements in hyperosmotic cultures, close phylogenetic relationship with B. thuringiensis, along with the reported property of encapsulation of insecticidal crystal proteins (Cry) in Bacillus strains and endophytic nature of B. mojavensis, strongly suggest that K90 could be a promising surrogate host for the transgenic delivery of "Cry" proteins.

Comparisons of dN/dS are time dependent for closely related bacterial genomes

The ratio of non-synonymous (dN) to synonymous (dS) changes between taxa is frequently computed to assay the strength and direction of selection. Here we note that for comparisons between closely related strains and/or species a second parameter needs to be considered, namely the time since divergence of the two sequences under scrutiny. We demonstrate that a simple time lag model provides a general, parsimonious explanation of the extensive variation in the dN/dS ratio seen when comparing closely related bacterial genomes. We explore this model through simulation and comparative genomics, and suggest a role for hitch-hiking in the accumulation of non-synonymous mutations. We also note taxon-specific differences in the change of dN/dS over time, which may indicate variation in selection, or in population genetics parameters such as population size or the rate of recombination. The effect of comparing intra-species polymorphism and inter-species substitution, and the problems associated with these concepts for asexual prokaryotes, are also discussed. We conclude that, because of the critical effect of time since divergence, inter-taxa comparisons are only possible by comparing trajectories of dN/dS over time and it is not valid to compare taxa on the basis of single time points.

Opinion: Re-evaluating prokaryotic species

There is no widely accepted concept of species for prokaryotes, and assignment of isolates to species is based on measures of phenotypic or genome similarity. The current methods for defining prokaryotic species are inadequate and incapable of keeping pace with the levels of diversity that are being uncovered in nature. Prokaryotic taxonomy is being influenced by advances in microbial population genetics, ecology and genomics, and by the ease with which sequence data can be obtained. Here, we review the classical approaches to prokaryotic species definition and discuss the current and future impact of multilocus nucleotide-sequence-based approaches to prokaryotic systematics. We also consider the potential, and difficulties, of assigning species status to biologically or ecologically meaningful sequence clusters.

The evolution of IncP catabolic plasmids

The recent adoption of whole plasmid genome sequencing as a routine analytical technique has provided the basis for cataloging the historical events through which plasmids are assembled from the available families of modular plasmid components. Horizontal gene transfer mediated by plasmids plays an important role in the adaptation of bacteria to the presence of specific metabolizable compounds, including man-made chemicals, in the surrounding environment. Bacterial plasmid genome sequence comparisons indicate that plasmids have complex genetic histories resulting from transposition, homologous recombination, and illegitimate recombinational events. Evidence from IncP plasmid genome sequences indicates that cryptic plasmid backbones acquire diverse catabolic pathways through gene capture and horizontal gene transfer.

The ecology and genetics of microbial diversity

Natural communities of microbes are often diverse, a fact that is difficult to reconcile with the action of natural selection in simple, uniform environments. We suggest that this apparent paradox may be resolved by considering the origin and fate of diversity in an explicitly ecological context. Here, we review insights into the ecological and genetic causes of diversity that stem from experiments with microbial populations evolving in the defined conditions of the laboratory environment. These studies highlight the importance of environmental structure in governing the fate of diversity and shed light on the genetic mechanisms generating diversity. We conclude by emphasizing the importance of placing detailed molecular-level studies within the context of a sound ecological and evolutionary framework.

Microarray-based resequencing of multiple Bacillus anthracis isolates

Custom-designed resequencing arrays were used to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication and by comparison to independently generated shotgun sequence

We used custom-designed resequencing arrays to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication (discrepancy rate of 7.4 × 10^-7) and by comparison to independently generated shotgun sequence (discrepancy rate < 2.5 × 10^-6). Population genomics studies of microbial pathogens using rapid resequencing technologies such as resequencing arrays are critical for recognizing newly emerging or genetically engineered strains.