FINE-SCALE GENETIC AND PHENOTYPIC STRUCTURE IN NATURAL POPULATIONS OF BACILLUS SUBTILIS AND BACILLUS LICHENIFORMIS: IMPLICATIONS FOR BACTERIAL EVOLUTION AND SPECIATION

Forty Years without Family: Three Novel Bacteriophages with High Similarity to SPP1 Reveal Decades of Evolutionary Stasis since the Isolation of Their Famous Relative

SPP1, an extensively studied bacteriophage of the Gram-positive Bacillus subtilis, is a model system for the study of phage-host interactions. Despite progress in the isolation and characterization of Bacillus phages, no previously fully sequenced phages have shared more than passing genetic similarity to SPP1. Here, we describe three virulent phages very similar to SPP1; SPP1 has greater than 80% nucleotide sequence identity and shares more that 85% of its protein coding genes with these phages. This is remarkable, given more than 40 years between the isolation of SPP1 and these phages. All three phages have somewhat larger genomes and more genes than SPP1. We identified a new putative gene in SPP1 based on a conserved sequence found in all phages. Gene conservation connotes purifying selection and is observed in structural genes and genes involved in DNA metabolism, but also in genes of unknown function, suggesting an important role in phage survival independent of the environment. Patterns of divergence point to genes or gene domains likely involved in adaptation to diverse hosts or different environments. Ultimately, comparative genomics of related phages provides insight into the long-term selective pressures that affect phage-bacteria interactions and alter phage genome content.

Complete Genome Sequences of Two Temperate Bacillus subtilis Phages Isolated at Tumamoc Hill Desert Laboratory

Bacteriophages are important in structuring bacterial communities, including desert soils dominated by Bacillus species. Here, we describe two genetically similar temperate phages isolated on a Bacillus subtilis strain from soil in Tucson, Arizona. Their double-stranded DNA (dsDNA) genomes contain 98 and 102 genes, with a set of 4 genes being found in only one phage.

Whole genome genetic variation and linkage disequilibrium in a diverse collection of Listeria monocytogenes isolates

We performed whole-genome multi-locus sequence typing for 2554 genes in a large and heterogenous panel of 180 Listeria monocytogenes strains having diverse geographical and temporal origins. The subtyping data was used for characterizing genetic variation and evaluating patterns of linkage disequilibrium in the pan-genome of L. monocytogenes. Our analysis revealed the presence of strong linkage disequilibrium in L. monocytogenes, with ~99% of genes showing significant non-random associations with a large majority of other genes in the genome. Twenty-seven loci having lower levels of association with other genes were considered to be potential “hot spots” for horizontal gene transfer (i.e., recombination via conjugation, transduction, and/or transformation). The patterns of linkage disequilibrium in L. monocytogenes suggest limited exchange of foreign genetic material in the genome and can be used as a tool for identifying new recombinant strains. This can help understand processes contributing to the diversification and evolution of this pathogenic bacteria, thereby facilitating development of effective control measures.

Genotyping of B. licheniformis based on a novel multi-locus sequence typing (MLST) scheme

BACKGROUND

Bacillus licheniformis has for many years been used in the industrial production of enzymes, antibiotics and detergents. However, as a producer of dormant heat-resistant endospores B. licheniformis might contaminate semi-preserved foods. The aim of this study was to establish a robust and novel genotyping scheme for B. licheniformis in order to reveal the evolutionary history of 53 strains of this species. Furthermore, the genotyping scheme was also investigated for its use to detect food-contaminating strains.

RESULTS

A multi-locus sequence typing (MLST) scheme, based on the sequence of six house-keeping genes (adk, ccpA, recF, rpoB, spo0A and sucC) of 53 B. licheniformis strains from different sources was established. The result of the MLST analysis supported previous findings of two different subgroups (lineages) within this species, named "A" and "B" Statistical analysis of the MLST data indicated a higher rate of recombination within group "A". Food isolates were widely dispersed in the MLST tree and could not be distinguished from the other strains. However, the food contaminating strain B. licheniformis NVH1032, represented by a unique sequence type (ST8), was distantly related to all other strains.

CONCLUSIONS

In this study, a novel and robust genotyping scheme for B. licheniformis was established, separating the species into two subgroups. This scheme could be used for further studies of evolution and population genetics in B. licheniformis.

Phylogeny in aid of the present and novel microbial lineages: diversity in Bacillus

Bacillus represents microbes of high economic, medical and biodefense importance. Bacillus strain identification based on 16S rRNA sequence analyses is invariably limited to species level. Secondly, certain discrepancies exist in the segregation of Bacillus subtilis strains. In the RDP/NCBI databases, out of a total of 2611 individual 16S rDNA sequences belonging to the 175 different species of the genus Bacillus, only 1586 have been identified up to species level. 16S rRNA sequences of Bacillus anthracis (153 strains), B. cereus (211 strains), B. thuringiensis (108 strains), B. subtilis (271 strains), B. licheniformis (131 strains), B. pumilus (83 strains), B. megaterium (47 strains), B. sphaericus (42 strains), B. clausii (39 strains) and B. halodurans (36 strains) were considered for generating species-specific framework and probes as tools for their rapid identification. Phylogenetic segregation of 1121, 16S rDNA sequences of 10 different Bacillus species in to 89 clusters enabled us to develop a phylogenetic frame work of 34 representative sequences. Using this phylogenetic framework, 305 out of 1025, 16S rDNA sequences presently classified as Bacillus sp. could be identified up to species level. This identification was supported by 20 to 30 nucleotides long signature sequences and in silico restriction enzyme analysis specific to the 10 Bacillus species. This integrated approach resulted in identifying around 30% of Bacillus sp. up to species level and revealed that B. subtilis strains can be segregated into two phylogenetically distinct groups, such that one of them may be renamed.

Importance of 3-hydroxy fatty acid composition of lipopeptides for biosurfactant activity

Biosurfactant production may be an economic approach to improving oil recovery. To obtain candidates most suitable for oil recovery, 207 strains, mostly belonging to the genus Bacillus, were tested for growth and biosurfactant production in medium with 5% NaCl under aerobic and anaerobic conditions. All strains grew aerobically with 5% NaCl, and 147 strains produced a biosurfactant. Thirty-five strains grew anaerobically with 5% NaCl, and two produced a biosurfactant. In order to relate structural differences to activity, eight lipopeptide biosurfactants with different specific activities produced by various Bacillus species were purified by a new protocol. The amino acid compositions of the eight lipopeptides were the same (Glu/Gln:Asp/Asn:Val:Leu, 1:1:1:4), but the fatty acid compositions differed. Multiple regression analysis showed that the specific biosurfactant activity depended on the ratios of both iso to normal even-numbered fatty acids and anteiso to iso odd-numbered fatty acids. A multiple regression model accurately predicted the specific biosurfactant activities of four newly purified biosurfactants (r2= 0.91). The fatty acid composition of the biosurfactant produced by Bacillus subtilis subsp. subtilis strain T89-42 was altered by the addition of branched-chain amino acids to the growth medium. The specific activities of biosurfactants produced in cultures with different amino acid additions were accurately predicted by the multiple regression model derived from the fatty acid compositions (r2= 0.95). Our work shows that many strains of Bacillus mojavensis and Bacillus subtilis produce biosurfactants and that the fatty acid composition is important for biosurfactant activity.

Adaptation and incipient sympatric speciation of Bacillus simplex under microclimatic contrast at "Evolution Canyons" I and II, Israel

The microevolutionary dynamics of prokaryotes in natural habitats, such as soil, is poorly understood in contrast to our increasing knowledge on their immense diversity. We performed microevolutionary analyses on 945 soil isolates of Bacillus simplex from "Evolution Canyons" I (Carmel, Israel) and II (Galilee, Israel). These canyons represent similar ecological replicates, separated by 40 km, with highly contrasting interslope abiotic and biotic conditions in each (within a distance of only 100-400 m). Strains representing genetic groups were identical in their 16S sequences, suggesting high genetic similarity and monophyletic origin. Parallel and nested phylogenetic structures correlated with ecological contrasts rather than geographical distance. Additionally, slope-specific populations differed substantially in their diversity. The levels of DNA repair (determined by UV sensitivity) and spontaneous mutation rate (resistance to rifampicin) relate to ecological stress and phylogeny. Altogether, the results suggest adaptive radiation at a microscale. We discuss the observed adaptive population structures in the context of incipient sympatric speciation in soil bacteria. We conclude that, despite different biology, prokaryotes, like sexually reproducing eukaryotes, may consist of true species and parallel ecological speciation in eukaryotes.

Anaerobic growth of Bacillus mojavensis and Bacillus subtilis requires deoxyribonucleosides or DNA

Bacillus mojavensis strains JF-2 (ATCC 39307), ROB2, and ABO21191(T) and Bacillus subtilis strains 168 (ATCC 23857) and ATCC 12332 required four deoxyribonucleosides or DNA for growth under strict anaerobic conditions. Bacillus licheniformis strains L89-11 and L87-11, Bacillus sonorensis strain TG8-8, and Bacillus cereus (ATCC 14579) did not require DNA for anaerobic growth. The requirement for the deoxyribonucleosides or DNA did not occur under aerobic growth conditions. The addition of a mixture of five nucleic acid bases, four ribonucleotides, or four ribonucleosides to the basal medium did not replace the requirement of B. mojavensis JF-2 for the four deoxyribonucleosides. However, the addition of salmon sperm DNA, herring sperm DNA, Escherichia coli DNA, or synthetic DNA (single or double stranded) to the basal medium supported anaerobic growth. The addition of four deoxyribonucleosides to the basal medium allowed aerobic growth of B. mojavensis JF-2 in the presence of hydroxyurea. B. mojavensis did not grow in DNA-supplemented basal medium that lacked sucrose as the energy source. These data provide strong evidence that externally supplied deoxyribonucleosides can be used to maintain a balanced deoxyribonucleotide pool for DNA synthesis and suggest that ribonucleotide reductases may not be essential to the bacterial cell cycle nor are they necessarily part of a minimal bacterial genome.

Genetic diversity, recombination and cryptic clades in Pseudomonas viridiflava infecting natural populations of Arabidopsis thaliana

Species-level genetic diversity and recombination in bacterial pathogens of wild plant populations have been nearly unexplored. Pseudomonas viridiflava is a common natural bacterial pathogen of Arabidopsis thaliana, for which pathogen defense genes and mechanisms are becoming increasing well known. The genetic variation contained within a worldwide sample of P. viridiflava collected from wild populations of A. thaliana was investigated using five genomic sequence fragments totaling 2.3 kb. Two distinct and deeply diverged clades were found within the P. viridiflava sample and in close proximity in multiple populations, each genetically diverse with synonymous variation as high as 9.3% in one of these clades. Within clades, there is evidence of frequent recombination within and between each sequenced locus and little geographic differentiation. Isolates from both clades were also found in a small sample of other herbaceous species in Midwest populations, indicating a possibly broad host range for P. viridiflava. The high levels of genetic variation and recombination together with a lack of geographic differentiation in this pathogen distinguish it from other bacterial plant pathogens for which intraspecific variation has been examined.

Genetic differentiation between sympatric populations of Bacillus cereus and Bacillus thuringiensis

Little is known about genetic exchanges in natural populations of bacteria of the spore-forming Bacillus cereus group, because no population genetics studies have been performed with local sympatric populations. We isolated strains of Bacillus thuringiensis and B. cereus from small samples of soil collected at the same time from two separate geographical sites, one within the forest and the other at the edge of the forest. A total of 100 B. cereus and 98 B. thuringiensis strains were isolated and characterized by electrophoresis to determine allelic composition at nine enzymatic loci. We observed genetic differentiation between populations of B. cereus and B. thuringiensis. Populations of a given Bacillus species--B. thuringiensis or B. cereus--were genetically more similar to each other than to populations of the other Bacillus species. Hemolytic activity provided further evidence of this genetic divergence, which remained evident even if putative clones were removed from the data set. Our results suggest that the rate of gene flow was higher between strains of the same species, but that exchanges between B. cereus and B. thuringiensis were nonetheless possible. Linkage disequilibrium analysis revealed sufficient recombination for B. cereus populations to be considered panmictic units. In B. thuringiensis, the balance between clonal proliferation and recombination seemed to depend on location. Overall, our data indicate that it is not important for risk assessment purposes to determine whether B. cereus and B. thuringiensis belong to a single or two species. Assessment of the biosafety of pest control based on B. thuringiensis requires evaluation of the extent of genetic exchange between strains in realistic natural conditions.

Genetic structure of natural populations of Escherichia coli in wild hosts on different continents

Current knowledge of genotypic and phenotypic diversity in the species Escherichia coli is based almost entirely on strains recovered from humans or zoo animals. In this study, we analyzed a collection of 202 strains obtained from 81 mammalian species representing 39 families and 14 orders in Australia and the Americas, as well as several reference strains; we also included a strain from a reptile and 10 from different families of birds collected in Mexico. The strains were characterized genotypically by multilocus enzyme electrophoresis (MLEE) and phenotypically by patterns of sugar utilization, antibiotic resistance, and plasmid profile. MLEE analysis yielded an estimated genetic diversity (H) of 0.682 for 11 loci. The observed genetic diversity in this sample is the greatest yet reported for E. coli. However, this genetic diversity is not randomly distributed; geographic effects and host taxonomic group accounted for most of the genetic differentiation. The genetic relationship among the strains showed that they are more associated by origin and host order than is expected by chance. In a dendrogram, the ancestral cluster includes primarily strains from Australia and ECOR strains from groups B and C. The most differentiated E. coli in our analysis are strains from Mexican carnivores and strains from humans, including those in the ECOR group A. The kinds and numbers of sugars utilized by the strains varied by host taxonomic group and country of origin. Strains isolated from bats were found to exploit the greatest range of sugars, while those from primates utilized the fewest. Toxins are more frequent in strains from rodents from both continents than in any other taxonomic group. Strains from Mexican wild mammals were, on average, as resistant to antibiotics as strains from humans in cities. On average, the Australian strains presented a lower antibiotic resistance than the Mexican strains. However, strains recovered from hosts in cities carried significantly more plasmids than did strains isolated from wild mammals. Previous studies have shown that natural populations of E. coli harbor an extensive genetic diversity that is organized in a limited number of clones. However, knowledge of this worldwide bacterium has been limited. Here, we suggest that the strains from a wide range of wild hosts from different regions of the world are organized in an ecotypic structure where adaptation to the host plays an important role in the population structure.